An optimal test with maximum average power while controlling FDR with application to RNA-seq data

The recent RNA-seq technology is an attractive method to study gene expression. One of the most important goals in RNA-seq data analysis is to detect genes differentially expressed across treatments. Although several statistical methods have been published, there are no theoretical justifications for whether these methods are optimal or how to search for the optimal test. Furthermore, most proposed tests are designed for testing whether the mean expression levels are exactly the same or not across treatments, whereas sometimes, biologists are interested in detecting genes with expression changes larger than a certain threshold. Another issue with current methods is that the false discovery rate (FDR) control is not well studied. In this manuscript, we propose a test to address all the above issues. Under model assumptions, we derive an optimal test that achieves the maximum of average power among those that control FDR at the same level. We also provide an approximated version, the approximated most average powerful (AMAP) test, for practical implementation. The proposed method allows for testing null hypotheses that are much more general than the ones most previous studies have considered, and it leads to a natural way of controlling the FDR. Through simulation studies, we show that our test has a higher power than other methods, including the widely-used edgeR, DESeq, and baySeq methods, as well as better FDR control than two other FDR control procedures commonly used in practice. For demonstration, we also apply the proposed method to a real RNA-seq dataset obtained from maize.

Prostate cancer invasion and metastasis: insights from mining genomic data

Prostate cancer (PCa) is the second most commonly diagnosed malignancy in men in the Western world and the second leading cause of cancer-related deaths among men worldwide. Although most cancers have the potential to metastasize under appropriate conditions, PCa favors the skeleton as a primary site of metastasis, suggesting that the bone microenvironment is conducive to its growth. PCa metastasis proceeds through a complex series of molecular events that include angiogenesis at the site of the original tumor, local migration within the primary site, intravasation into the blood stream, survival within the circulation, extravasation of the tumor cells to the target organ and colonization of those cells within the new site. In turn, each one of these steps involves a complicated chain of events that utilize multiple protein-protein interactions, protein signaling cascades and transcriptional changes. Despite the urgent need to improve current biomarkers for diagnosis, prognosis and drug resistance, advances have been slow. Global gene expression methods such as gene microarrays and RNA sequencing enable the study of thousands of genes simultaneously and allow scientists to examine molecular pathways of cancer pathogenesis. In this review, we summarize the current literature that explored high-throughput transcriptome analysis toward the advancement of biomarker discovery for PCa. Novel biomarkers are strongly needed to enable more accurate detection of PCa, improve prediction of tumor aggressiveness and facilitate the discovery of new therapeutic targets for tailored medicine. Promising molecular markers identified from gene expression profiling studies include HPN, CLU1, WT1, WNT5A, AURKA and SPARC.

Impacts of variation in the human genome on gene regulation

Recent advances in fast and inexpensive DNA sequencing have enabled the extensive study of genomic and transcriptomic variation in humans. Human genomic variation is composed of sequence and structural changes including single-nucleotide and multinucleotide variants, short insertions or deletions (indels), larger copy number variants, and similarly sized copy neutral inversions and translocations. It is now well established that any two genomes differ extensively and that structural changes constitute a prominent source of this variation. There have also been major technological advances in RNA sequencing to globally quantify and describe diversity in transcripts. Large consortia such as the 1000 Genomes Project and the ENCODE (ENCyclopedia Of DNA Elements) Project are producing increasingly comprehensive maps outlining the regions of the human genome containing variants and functional elements, respectively. Integration of genetic variation data and extensive annotation of functional genomic elements, along with the ability to measure global transcription, allow the impacts of genetic variants on gene expression to be resolved. There are several well-established models by which genetic variants affect gene regulation depending on the type, nature, and position of the variant with respect to the affected genes. These effects can be manifested in two ways: changes to transcript sequences and isoforms by coding variants, and changes to transcript abundance by dosage or regulatory variants. Here, we review the current state of how genetic variations impact gene regulation locally and globally in the human genome.

MicroRNA discovery by similarity search to a database of RNA-seq profiles

In silico generated search for microRNAs (miRNAs) has been driven by methods compiling structural features of the miRNA precursor hairpin, as well as to some degree combining this with the analysis of RNA-seq profiles for which the miRNA typically leave the drosha/dicer fingerprint of 1-2 ~22 nt blocks of reads corresponding to the mature and star miRNA. In complement to the previous methods, we present a study where we systematically exploit these patterns of read profiles. We created two datasets comprised of 2540 and 4795 read profiles obtained after preprocessing short RNA-seq data from miRBase and ENCODE, respectively. Out of 4795 ENCODE read profiles, 1361 are annotated as non-coding RNAs (ncRNAs) and of which 285 are further annotated as miRNAs. Using deepBlockAlign (dba), we align ncRNA read profiles from ENCODE against the miRBase read profiles (cleaned for "self-matches") and are able to separate ENCODE miRNAs from the other ncRNAs by a Matthews Correlation Coefficient (MCC) of 0.8 and obtain an area under the curve of 0.93. Based on the dba score cut-off of 0.7 at which we observed the maximum MCC of 0.8, we predict 523 novel miRNA candidates. An additional RNA secondary structure analysis reveal that 42 of the candidates overlap with predicted conserved secondary structure. Further analysis reveal that the 523 miRNA candidates are located in genomic regions with MAF block (UCSC) fragmentation and poor sequence conservation, which in part might explain why they have been overlooked in previous efforts. We further analyzed known human and mouse miRNA read profiles and found two distinct classes; the first containing two blocks and the second containing >2 blocks of reads. Also the latter class holds read profiles that have less well defined arrangement of reads in comparison to the former class. On comparison of miRNA read profiles from plants and animals, we observed kingdom specific read profiles that are distinct in terms of both length and distribution of reads within the read profiles to each other. All the data, as well as a server to search miRBase read profiles by uploading a BED file, is available at http://rth.dk/resources/mirdba.

The effects of MucR on expression of type IV secretion system, quorum sensing system and stress responses in Brucella melitensis

MucR is a transcriptional regulator in many bacterial pathogens and is required for virulence in mice and macrophages, resistance to stress responses, and modification of the cell envelope in Brucella spp. To determine why the mucR deleted mutant is attenuated in vivo and in vitro, we performed RNA-seq analysis using Brucella melitensis RNA obtained from B. melitensis 16M and 16MΔmucR grown under the same conditions. We found 442 differentially expressed genes; 310 were over expressed, and 132 were less expressed in 16MΔmucR. Many genes identified are involved in metabolism, cell wall/envelope biogenesis, replication, and translation. Notably, genes involved in type IV secretion system and quorum sensing system were down-regulated in 16MΔmucR. In addition, genes involved in tolerance to acid and iron-limitation were also affected and experimentally verified in this study. The effects of MucR on Brucella survival and persistence in mice and macrophages were related to type IV secretion system, quorum sensing system, and stress tolerance, which also provide added insight to the MucR regulon.

The family that eats together stays together: new p53 family transcriptional targets in autophagy

Autophagy is a biological process that is crucial to maintain cellular homeostasis and is regulated by several metabolic pathways, including the p53 tumor suppressor pathway. In this issue of Genes & Development, Kenzelmann Broz and colleagues (pp. 1016-1031) show how the p53 family as a whole, including p63 and p73, collaborate in controlling autophagy to support tumor suppression.

Global genomic profiling reveals an extensive p53-regulated autophagy program contributing to key p53 responses

The mechanisms by which the p53 tumor suppressor acts remain incompletely understood. To gain new insights into p53 biology, we used high-throughput sequencing to analyze global p53 transcriptional networks in primary mouse embryo fibroblasts in response to DNA damage. Chromatin immunoprecipitation sequencing reveals 4785 p53-bound sites in the genome located near 3193 genes involved in diverse biological processes. RNA sequencing analysis shows that only a subset of p53-bound genes is transcriptionally regulated, yielding a list of 432 p53-bound and regulated genes. Interestingly, we identify a host of autophagy genes as direct p53 target genes. While the autophagy program is regulated predominantly by p53, the p53 family members p63 and p73 contribute to activation of this autophagy gene network. Induction of autophagy genes in response to p53 activation is associated with enhanced autophagy in diverse settings and depends on p53 transcriptional activity. While p53-induced autophagy does not affect cell cycle arrest in response to DNA damage, it is important for both robust p53-dependent apoptosis triggered by DNA damage and transformation suppression by p53. Together, our data highlight an intimate connection between p53 and autophagy through a vast transcriptional network and indicate that autophagy contributes to p53-dependent apoptosis and cancer suppression.

Dynamic binding of RBPJ is determined by Notch signaling status

Notch signaling plays crucial roles in mediating cell fate choices in all metazoans largely by specifying the transcriptional output of one cell in response to a neighboring cell. The DNA-binding protein RBPJ is the principle effector of this pathway in mammals and, together with the transcription factor moiety of Notch (NICD), regulates the expression of target genes. The prevalent view presumes that RBPJ statically occupies consensus binding sites while exchanging repressors for activators in response to NICD. We present the first specific RBPJ chromatin immunoprecipitation and high-throughput sequencing study in mammalian cells. To dissect the mode of transcriptional regulation by RBPJ and identify its direct targets, whole-genome binding profiles were generated for RBPJ; its coactivator, p300; NICD; and the histone H3 modifications H3 Lys 4 trimethylation (H3K4me3), H3 Lys 4 monomethylation (H3K4me1), and histone H3 Lys 27 acetylation (H3K27ac) in myogenic cells under active or inhibitory Notch signaling conditions. Our results demonstrate dynamic binding of RBPJ in response to Notch activation at essentially all sites co-occupied by NICD. Additionally, we identify a distinct set of sites where RBPJ recruits neither NICD nor p300 and binds DNA statically, irrespective of Notch activity. These findings significantly modify our views on how RBPJ and Notch signaling mediate their activities and consequently impact on cell fate decisions.

Human housekeeping genes, revisited

Housekeeping genes are involved in basic cell maintenance and, therefore, are expected to maintain constant expression levels in all cells and conditions. Identification of these genes facilitates exposure of the underlying cellular infrastructure and increases understanding of various structural genomic features. In addition, housekeeping genes are instrumental for calibration in many biotechnological applications and genomic studies. Advances in our ability to measure RNA expression have resulted in a gradual increase in the number of identified housekeeping genes. Here, we describe housekeeping gene detection in the era of massive parallel sequencing and RNA-seq. We emphasize the importance of expression at a constant level and provide a list of 3804 human genes that are expressed uniformly across a panel of tissues. Several exceptionally uniform genes are singled out for future experimental use, such as RT-PCR control genes. Finally, we discuss both ways in which current technology can meet some of past obstacles encountered, and several as yet unmet challenges.

Deep sequencing of the transcriptome reveals inflammatory features of porcine visceral adipose tissue

Functional differences in the different types of adipose tissue and the impact of their dysfunction on metabolism are associated with the regional distribution of adipose depots. Here we show a genome-wide comparison between the transcriptomes of one source of subcutaneous and two sources of visceral adipose tissue in the pig using an RNA-seq approach. We obtained ~32.3 million unique mapped reads which covered ~80.2% of the current annotated transcripts across these three sources of adipose tissue. We identified various genes differentially expressed between subcutaneous and visceral adipose tissue, which are potentially associated with the inflammatory features of visceral adipose tissue. These results are of benefit for understanding the phenotypic, metabolic and functional differences between different types of adipose tissue that are deposited in different body sites.

A quantitative reference transcriptome for Nematostella vectensis early embryonic development: a pipeline for de novo assembly in emerging model systems

BACKGROUND

The de novo assembly of transcriptomes from short shotgun sequences raises challenges due to random and non-random sequencing biases and inherent transcript complexity. We sought to define a pipeline for de novo transcriptome assembly to aid researchers working with emerging model systems where well annotated genome assemblies are not available as a reference. To detail this experimental and computational method, we used early embryos of the sea anemone, Nematostella vectensis, an emerging model system for studies of animal body plan evolution. We performed RNA-seq on embryos up to 24 h of development using Illumina HiSeq technology and evaluated independent de novo assembly methods. The resulting reads were assembled using either the Trinity assembler on all quality controlled reads or both the Velvet and Oases assemblers on reads passing a stringent digital normalization filter. A control set of mRNA standards from the National Institute of Standards and Technology (NIST) was included in our experimental pipeline to invest our transcriptome with quantitative information on absolute transcript levels and to provide additional quality control.

RESULTS

We generated >200 million paired-end reads from directional cDNA libraries representing well over 20 Gb of sequence. The Trinity assembler pipeline, including preliminary quality control steps, resulted in more than 86% of reads aligning with the reference transcriptome thus generated. Nevertheless, digital normalization combined with assembly by Velvet and Oases required far less computing power and decreased processing time while still mapping 82% of reads. We have made the raw sequencing reads and assembled transcriptome publically available.

CONCLUSIONS

Nematostella vectensis was chosen for its strategic position in the tree of life for studies into the origins of the animal body plan, however, the challenge of reference-free transcriptome assembly is relevant to all systems for which well annotated gene models and independently verified genome assembly may not be available. To navigate this new territory, we have constructed a pipeline for library preparation and computational analysis for de novo transcriptome assembly. The gene models defined by this reference transcriptome define the set of genes transcribed in early Nematostella development and will provide a valuable dataset for further gene regulatory network investigations.

A digital compendium of genes mediating the reversible phosphorylation of proteins in fe-deficient Arabidopsis roots

Post-translational modifications of proteins such as reversible phosphorylation provide an important but understudied regulatory network that controls important nodes in the adaptation of plants to environmental conditions. Iron (Fe) is an essential mineral nutrient for plants, but due to its low solubility often a limiting factor for optimal growth. To understand the role of protein phosphorylation in the regulation of cellular Fe homeostasis, we analyzed the expression of protein kinases (PKs) and phosphatases (PPs) in Arabidopsis roots by mining differentially expressed PK and PP genes. Transcriptome analysis using RNA-seq revealed that subsets of 203 PK and 39 PP genes were differentially expressed under Fe-deficient conditions. Functional modules of these PK and PP genes were further generated based on co-expression analysis using the MACCU toolbox on the basis of 300 publicly available root-related microarray data sets. Results revealed networks comprising 87 known or annotated PK and PP genes that could be subdivided into one large and several smaller highly co-expressed gene modules. The largest module was composed of 58 genes, most of which have been assigned to the leucine-rich repeat protein kinase superfamily and associated with the biological processes "hypotonic salinity response," "potassium ion import," and "cellular potassium ion homeostasis." The comprehensive transcriptional information on PK and PP genes in iron-deficient roots provided here sets the stage for follow-up experiments and contributes to our understanding of the post-translational regulation of Fe deficiency and potassium ion homeostasis.

Mapping Novel Pathways in Cardiovascular Disease Using eQTL Data: The Past, Present, and Future of Gene Expression Analysis

Genome-wide association studies (GWAS) have identified genetic variants associated with numerous cardiovascular and metabolic diseases. Newly identified polymorphisms associated with myocardial infarction, dyslipidemia, hypertension, diabetes, and insulin resistance suggest novel mechanistic pathways that underlie these and other complex diseases. Working out the connections between the polymorphisms identified in GWAS and their biological mechanisms has been especially challenging given the number of non-coding variants identified thus far. In this review, we discuss the utility of expression quantitative trait locus (eQTL) databases in the study of non-coding variants with respect to cardiovascular and metabolic phenotypes. Recent successes in using eQTL data to link variants with functional candidate genes will be reviewed, and the shortcomings of this approach will be outlined. Finally, we discuss the emerging next generation of eQTL studies that take advantage of the ability to generate induced pluripotent stem cell lines from population cohorts.

Mapping the RNA-Seq trash bin: unusual transcripts in prokaryotic transcriptome sequencing data

Prokaryotic transcripts constitute almost always uninterrupted intervals when mapped back to the genome. Split reads, i.e., RNA-seq reads consisting of parts that only map to discontiguous loci, are thus disregarded in most analysis pipelines. There are, however, some well-known exceptions, in particular, tRNA splicing and circularized small RNAs in Archaea as well as self-splicing introns. Here, we reanalyze a series of published RNA-seq data sets, screening them specifically for non-contiguously mapping reads. We recover most of the known cases together with several novel archaeal ncRNAs associated with circularized products. In Eubacteria, only a handful of interesting candidates were obtained beyond a few previously described group I and group II introns. Most of the atypically mapping reads do not appear to correspond to well-defined, specifically processed products. Whether this diffuse background is, at least in part, an incidental by-product of prokaryotic RNA processing or whether it consists entirely of technical artifacts of reverse transcription or amplification remains unknown.

eQTL Mapping Using RNA-seq Data

As RNA-seq is replacing gene expression microarrays to assess genome-wide transcription abundance, gene expression Quantitative Trait Locus (eQTL) studies using RNA-seq have emerged. RNA-seq delivers two novel features that are important for eQTL studies. First, it provides information on allele-specific expression (ASE), which is not available from gene expression microarrays. Second, it generates unprecedentedly rich data to study RNA-isoform expression. In this paper, we review current methods for eQTL mapping using ASE and discuss some future directions. We also review existing works that use RNA-seq data to study RNA-isoform expression and we discuss the gaps between these works and isoform-specific eQTL mapping.

Evolution of the eye transcriptome under constant darkness in Sinocyclocheilus cavefish

In adaptating to perpetual darkness, cave species gradually lose eyes and body pigmentation and evolve alternatives for exploring their environments. Although troglodyte features evolved independently many times in cavefish, we do not yet know whether independent evolution of these characters involves common genetic mechanisms. Surface-dwelling and many cave-dwelling species make the freshwater teleost genus Sinocyclocheilus an excellent model for studying the evolution of adaptations to life in constant darkness. We compared the mature retinal histology of surface and cave species in Sinocyclocheilus and found that adult cavefish showed a reduction in the number and length of photoreceptor cells. To identify genes and genetic pathways that evolved in constant darkness, we used RNA-seq to compare eyes of surface and cave species. De novo transcriptome assemblies were developed for both species, and contigs were annotated with gene ontology. Results from cave-dwelling Sinocyclocheilus revealed reduced transcription of phototransduction and other genes important for retinal function. In contrast to the blind Mexican tetra cavefish Astyanax mexicanus, our results on morphologies and gene expression suggest that evolved retinal reduction in cave-dwelling Sinocyclocheilus occurs in a lens-independent fashion by the reduced proliferation and downregulation of transcriptional factors shown to have direct roles in retinal development and maintenance, including cone-rod homeobox (crx) and Wnt pathway members. These results show that the independent evolution of retinal degeneration in cavefish can occur by different developmental genetic mechanisms.

Comparative systems biology between human and animal models based on next-generation sequencing methods

Animal models provide myriad benefits to both experimental and clinical research. Unfortunately, in many situations, they fall short of expected results or provide contradictory results. In part, this can be the result of traditional molecular biological approaches that are relatively inefficient in elucidating underlying molecular mechanism. To improve the efficacy of animal models, a technological breakthrough is required. The growing availability and application of the high-throughput methods make systematic comparisons between human and animal models easier to perform. In the present study, we introduce the concept of the comparative systems biology, which we define as "comparisons of biological systems in different states or species used to achieve an integrated understanding of life forms with all their characteristic complexity of interactions at multiple levels". Furthermore, we discuss the applications of RNA-seq and ChIP-seq technologies to comparative systems biology between human and animal models and assess the potential applications for this approach in the future studies.

Phenotypic and molecular differences between rats selectively bred to voluntarily run high vs. low nightly distances

The purpose of the present study was to partially phenotype male and female rats from generations 8-10 (G8-G10) that had been selectively bred to possess low (LVR) vs. high voluntary running (HVR) behavior. Over the first 6 days with wheels, 34-day-old G8 male and female LVRs ran shorter distances (P < 0.001), spent less time running (P < 0.001), and ran slower (P < 0.001) than their G8 male and female HVR counterparts, respectively. HVR and LVR lines consumed similar amounts of standard chow with or without wheels. No inherent difference existed in PGC-1α mRNA in the plantaris and soleus muscles of LVR and HVR nonrunners, although G8 LVR rats inherently possessed less NADH-positive superficial plantaris fibers compared with G8 HVR rats. While day 28 body mass tended to be greater in both sexes of G9-G10 LVR nonrunners vs. G9-G10 HVR nonrunners (P = 0.06), body fat percentage was similar between lines. G9-G10 HVRs had fat mass loss after 6 days of running compared with their prerunning values, while LVR did not lose or gain fat mass during the 6-day voluntary running period. RNA deep sequencing efforts in the nucleus accumbens showed only eight transcripts to be >1.5-fold differentially expressed between lines in HVR and LVR nonrunners. Interestingly, HVRs presented less Oprd1 mRNA, which ties in to potential differences in dopaminergic signaling between lines. This unique animal model provides further evidence as to how exercise may be mechanistically regulated.

Mycobacterium tuberculosis Transcriptome Profiling in Mice with Genetically Different Susceptibility to Tuberculosis

Whole transcriptome profiling is now almost routinely used in various fields of biology, including microbiology. In vivo transcriptome studies usually provide relevant information about the biological processes in the organism and thus are indispensable for the formulation of hypotheses, testing, and correcting. In this study, we describe the results of genome-wide transcriptional profiling of the major human bacterial pathogen M. tuberculosis during its persistence in lungs. Two mouse strains differing in their susceptibility to tuberculosis were used for experimental infection with M. tuberculosis. Mycobacterial transcriptomes obtained from the infected tissues of the mice at two different time points were analyzed by deep sequencing and compared. It was hypothesized that the changes in the M. tuberculosis transcriptome may attest to the activation of the metabolism of lipids and amino acids, transition to anaerobic respiration, and increased expression of the factors modulating the immune response. A total of 209 genes were determined whose expression increased with disease progression in both host strains (commonly upregulated genes, CUG). Among them, the genes related to the functional categories of lipid metabolism, cell wall, and cell processes are of great interest. It was assumed that the products of these genes are involved in M. tuberculosis adaptation to the host immune system defense, thus being potential targets for drug development.

A new method for stranded whole transcriptome RNA-seq

This report describes an improved protocol to generate stranded, barcoded RNA-seq libraries to capture the whole transcriptome. By optimizing the use of duplex specific nuclease (DSN) to remove ribosomal RNA reads from stranded barcoded libraries, we demonstrate improved efficiency of multiplexed next generation sequencing (NGS). This approach detects expression profiles of all RNA types, including miRNA (microRNA), piRNA (Piwi-interacting RNA), snoRNA (small nucleolar RNA), lincRNA (long non-coding RNA), mtRNA (mitochondrial RNA) and mRNA (messenger RNA) without the use of gel electrophoresis. The improved protocol generates high quality data that can be used to identify differential expression in known and novel coding and non-coding transcripts, splice variants, mitochondrial genes and SNPs (single nucleotide polymorphisms).

Catalyzing plant science research with RNA-seq

Next generation DNA sequencing technologies are driving increasingly rapid, affordable and high resolution analyses of plant transcriptomes through sequencing of their associated cDNA (complementary DNA) populations; an analytical platform commonly referred to as RNA-sequencing (RNA-seq). Since entering the arena of whole genome profiling technologies only a few years ago, RNA-seq has proven itself to be a powerful tool with a remarkably diverse range of applications, from detailed studies of biological processes at the cell type-specific level, to providing insights into fundamental questions in plant biology on an evolutionary time scale. Applications include generating genomic data for heretofore unsequenced species, thus expanding the boundaries of what had been considered "model organisms," elucidating structural and regulatory gene networks, revealing how plants respond to developmental cues and their environment, allowing a better understanding of the relationships between genes and their products, and uniting the "omics" fields of transcriptomics, proteomics, and metabolomics into a now common systems biology paradigm. We provide an overview of the breadth of such studies and summarize the range of RNA-seq protocols that have been developed to address questions spanning cell type-specific-based transcriptomics, transcript secondary structure and gene mapping.

EDGE-pro: Estimated Degree of Gene Expression in Prokaryotic Genomes

Background

The expression levels of bacterial genes can be measured directly using next-generation sequencing (NGS) methods, offering much greater sensitivity and accuracy than earlier, microarray-based methods. Most bioinformatics software for estimating levels of gene expression from NGS data has been designed for eukaryotic genomes, with algorithms focusing particularly on detection of splicing patterns. These methods do not perform well on bacterial genomes.

Results

Here we describe the first software system designed explicitly for quantifying the degree of gene expression in bacteria and other prokaryotes. EDGE-pro (Estimated Degree of Gene Expression in PROkaryotes) processes the raw data from an RNA-seq experiment on a bacterial or archaeal species and produces estimates of the expression levels for each gene in these gene-dense genomes.

Software

The EDGE-pro tool is implemented as a pipeline of C++ and Perl programs and is freely available as open-source code at http://www.genomics.jhu.edu/software/EDGE/index.shtml.

Investigation of Aspergillus fumigatus biofilm formation by various "omics" approaches

In the lung, Aspergillus fumigatus usually forms a dense colony of filaments embedded in a polymeric extracellular matrix called biofilm (BF). This extracellular matrix embeds and glues hyphae together and protects the fungus from an outside hostile environment. This extracellular matrix is absent in fungal colonies grown under classical liquid shake conditions (PL), which were historically used to understand A. fumigatus pathobiology. Recent works have shown that the fungus in this aerial grown BF-like state exhibits reduced susceptibility to antifungal drugs and undergoes major metabolic changes that are thought to be associated to virulence. These differences in pathological and physiological characteristics between BF and liquid shake conditions suggest that the PL condition is a poor in vitro disease model. In the laboratory, A. fumigatus mycelium embedded by the extracellular matrix can be produced in vitro in aerial condition using an agar-based medium. To provide a global and accurate understanding of A. fumigatus in vitro BF growth, we utilized microarray, RNA-sequencing, and proteomic analysis to compare the global gene and protein expression profiles of A. fumigatus grown under BF and PL conditions. In this review, we will present the different signatures obtained with these three “omics” methods. We will discuss the advantages and limitations of each method and their complementarity.

Autophagy in Antarctica: combating dehydration stress in the world's southernmost insect

The midge Belgica antarctica is the only insect endemic to Antarctica and has the southernmost range of any insect. In its natural environment, B. antarctica frequently faces desiccating conditions, as environmental water is frozen for up to 9 months annually. The molecular mechanisms by which B. antarctica tolerates extreme dehydration are poorly understood, but recent work from our laboratory reports genome-wide expression changes in response to extreme dehydration (~40% water loss), the first genome-scale transcriptome reported for an Antarctic animal. Among transcripts differentially regulated during dehydration, there is coordinated upregulation of numerous genes involved in autophagy, including genes responsible for autophagosome synthesis and autophagy-associated transcription factors. Also, several genes and pathways that interact with and regulate autophagy, e.g., sestrins and proteasomal genes, are concurrently upregulated. This suggests that autophagy and related processes are key elements regulating stress tolerance in this extreme environment.

Identification and characterization of small RNAs in Yersinia pestis

Yersinia pestis, the etiologic agent of plague, is closely related to Yersinia pseudotuberculosis evolutionarily but has a very different mode of infection. The RNA-binding regulatory protein, Hfq, mediates regulation by small RNAs (sRNAs) and is required for virulence of both Y. pestis and Y. pseudotuberculosis. Moreover, Hfq is required for growth of Y. pestis, but not Y. pseudotuberculosis, at 37°C. Together, these observations suggest that sRNAs play important roles in the virulence and survival of Y. pestis, and that regulation by sRNAs may account for some of the differences between Y. pestis and Y. pseudotuberculosis. We have used a deep sequencing approach to identify 31 sRNAs in Y. pestis. The majority of these sRNAs are not conserved outside the Yersiniae. Expression of the sRNAs was confirmed by Northern analysis and we developed deep sequencing approaches to map 5ʹ and 3ʹ ends of many sRNAs simultaneously. Expression of the majority of the sRNAs we identified is dependent upon Hfq. We also observed temperature-dependent effects on the expression of many sRNAs, and differences in expression patterns between Y. pestis and Y. pseudotuberculosis. Thus, our data suggest that regulation by sRNAs plays an important role in the lifestyle switch from flea to mammalian host, and that regulation by sRNAs may contribute to the phenotypic differences between Y. pestis and Y. pseudotuberculosis.

Next-generation sequencing technologies and their impact on microbial genomics

Next-generation sequencing technologies have had a dramatic impact in the field of genomic research through the provision of a low cost, high-throughput alternative to traditional capillary sequencers. These new sequencing methods have surpassed their original scope and now provide a range of utility-based applications, which allow for a more comprehensive analysis of the structure and content of microbial genomes than was previously possible. With the commercialization of a third generation of sequencing technologies imminent, we discuss the applications of current next-generation sequencing methods and explore their impact on and contribution to microbial genome research.

Dual functions of TAF7L in adipocyte differentiation

The diverse transcriptional mechanisms governing cellular differentiation and development of mammalian tissue remains poorly understood. Here we report that TAF7L, a paralogue of TFIID subunit TAF7, is enriched in adipocytes and white fat tissue (WAT) in mouse. Depletion of TAF7L reduced adipocyte-specific gene expression, compromised adipocyte differentiation, and WAT development as well. Ectopic expression of TAF7L in myoblasts reprograms these muscle precursors into adipocytes upon induction. Genome-wide mRNA-seq expression profiling and ChIP-seq binding studies confirmed that TAF7L is required for activating adipocyte-specific genes via a dual mechanism wherein it interacts with PPARγ at enhancers and TBP/Pol II at core promoters. In vitro binding studies confirmed that TAF7L forms complexes with both TBP and PPARγ. These findings suggest that TAF7L plays an integral role in adipocyte gene expression by targeting enhancers as a cofactor for PPARγ and promoters as a component of the core transcriptional machinery.

DOI:http://dx.doi.org/10.7554/eLife.00170.001

The development of a single fertilized egg into a highly complex animal is determined by its genome, with a process called differential gene regulation exerting exquisite control over gene expression to ensure that various specialized cells are generated and that many types of tissue are produced. However, the mechanisms responsible for controlling gene expression and, therefore mammalian development, are poorly understood.

Researchers have developed a number of in vitro cell culture models to elucidate the details of differential gene regulation, and this approach has been used to characterize adipocytes—cells that store energy in the form of fat—for close to two decades. The formation of adipocytes, a process known as adipogenesis, has been extensively studied, but there remain major gaps in our knowledge: for example, the identities of many of the transcriptional regulators that are responsible for the differentiation of mesenchymal stem cells into adipocytes remain a mystery. This task is complicated by the fact that some of these regulators are involved in the differentiation of multiple cell lines, and that some of them also have multiple roles in the generation of a single cell type. In addition to being of fundamental interest, improving our knowledge of the properties and behavior of adipocytes is essential for tackling the increasing prevalence of obesity in the developed world.

Zhou et al. now report that TAF7L—a gene that was previously thought to be involved only in the production of sperm cells—has two roles in the differentiation of stem cells to form adipocytes. Using a combination of cellular, biochemical, genetic and genomic techniques, they show that TAF7L interacts with PPARγ, an important adipocyte transcriptional regulator at enhancer sites on the genome to increase the transcription of genes that are involved in adipogenesis. They also show that TAF7L interacts with a general transcription factor called TBP (short for TATA-binding protein) at promoter sequences, again to increase the expression of genes involved in adipogenesis. Moreover, they show that the expression of TAF7L in myoblasts—precursor cells that usually become muscle cells—can induce the formation of fat cells rather than muscle cells. Furthermore, mice lacking TAF7L are lean compared to their normal littermates. A clearer understanding of the underlying causes of fat cell formation could lead to the development of new approaches for the treatment of obesity and associated diseases.

DOI:http://dx.doi.org/10.7554/eLife.00170.002

The fate of duplicated genes in a polyploid plant genome

Polyploidy is generally not tolerated in animals, but is widespread in plant genomes and may result in extensive genetic redundancy. The fate of duplicated genes is poorly understood, both functionally and evolutionarily. Soybean (Glycine max L.) has undergone two separate polyploidy events (13 and 59 million years ago) that have resulted in 75% of its genes being present in multiple copies. It therefore constitutes a good model to study the impact of whole-genome duplication on gene expression. Using RNA-seq, we tested the functional fate of a set of approximately 18 000 duplicated genes. Across seven tissues tested, approximately 50% of paralogs were differentially expressed and thus had undergone expression sub-functionalization. Based on gene ontology and expression data, our analysis also revealed that only a small proportion of the duplicated genes have been neo-functionalized or non-functionalized. In addition, duplicated genes were often found in collinear blocks, and several blocks of duplicated genes were co-regulated, suggesting some type of epigenetic or positional regulation. We also found that transcription factors and ribosomal protein genes were differentially expressed in many tissues, suggesting that the main consequence of polyploidy in soybean may be at the regulatory level.

Widespread Exonization of Transposable Elements in Human Coding Sequences is Associated with Epigenetic Regulation of Transcription

Background

Transposable Elements (TEs) have long been regarded as selfish or junk DNA having little or no role in the regulation or functioning of the human genome. However, over the past several years this view came to be challenged as several studies provided anecdotal as well as global evidence for the contribution of TEs to the regulatory and coding needs of human genes. In this study, we explored the incorporation and epigenetic regulation of coding sequences donated by TEs using gene expression and other ancillary genomics data from two human hematopoietic cell-lines: GM12878 (a lymphoblastoid cell line) and K562 (a Chronic Myelogenous Leukemia cell line). In each cell line, we found several thousand instances of TEs donating coding sequences to human genes. We compared the transcriptome assembly of the RNA sequencing (RNA-Seq) reads with and without the aid of a reference transcriptome and found that the percentage of genes that incorporate TEs in their coding sequences is significantly greater than that obtained from the reference transcriptome assemblies using Refseq and Gencode gene models. We also used histone modifications chromatin immunoprecipitation sequencing (ChIP-Seq) data, Cap Analysis of Gene Expression (CAGE) data and DNAseI Hypersensitivity Site (DHS) data to demonstrate the epigenetic regulation of the TE derived coding sequences. Our results suggest that TEs form a significantly higher percentage of coding sequences than represented in gene annotation databases and these TE derived sequences are epigenetically regulated in accordance with their expression in the two cell types.

The transcriptional response of Arabidopsis leaves to Fe deficiency

Due to its ease to donate or accept electrons, iron (Fe) plays a crucial role in respiration and metabolism, including tetrapyrrole synthesis, in virtually all organisms. In plants, Fe is a component of the photosystems and thus essential for photosynthesis. Fe deficiency compromises chlorophyll (Chl) synthesis, leading to interveinal chlorosis in developing leaves and decreased photosynthetic activity. To gain insights into the responses of photosynthetically active cells to Fe deficiency, we conducted transcriptional profiling experiments on leaves from Fe-sufficient and Fe-deficient plants using the RNA-seq technology. As anticipated, genes associated with photosynthesis and tetrapyrrole metabolism were dramatically down-regulated by Fe deficiency. A sophisticated response comprising the down-regulation of HEMA1 and NYC1, which catalyze the first committed step in tetrapyrrole biosynthesis and the conversion of Chl b to Chl a at the commencement of Chl breakdown, respectively, and the up-regulation of CGLD27, which is conserved in plastid-containing organisms and putatively involved in xanthophyll biosynthesis, indicates a carefully orchestrated balance of potentially toxic tetrapyrrole intermediates and functional end products to avoid photo-oxidative damage. Comparing the responses to Fe deficiency in leaves to that in roots confirmed subgroup 1b bHLH transcription factors and POPEYE/BRUTUS as important regulators of Fe homeostasis in both leaf and root cells, and indicated six novel players with putative roles in Fe homeostasis that were highly expressed in leaves and roots and greatly induced by Fe deficiency. The data further revealed down-regulation of organ-specific subsets of genes encoding ribosomal proteins, which may be indicative of a change in ribosomal composition that could bias translation. It is concluded that Fe deficiency causes a massive reorganization of plastid activity, which is adjusting leaf function to the availability of Fe.

Chromatin is an ancient innovation conserved between Archaea and Eukarya

The eukaryotic nucleosome is the fundamental unit of chromatin, comprising a protein octamer that wraps ∼147 bp of DNA and has essential roles in DNA compaction, replication and gene expression. Nucleosomes and chromatin have historically been considered to be unique to eukaryotes, yet studies of select archaea have identified homologs of histone proteins that assemble into tetrameric nucleosomes. Here we report the first archaeal genome-wide nucleosome occupancy map, as observed in the halophile Haloferax volcanii. Nucleosome occupancy was compared with gene expression by compiling a comprehensive transcriptome of Hfx. volcanii. We found that archaeal transcripts possess hallmarks of eukaryotic chromatin structure: nucleosome-depleted regions at transcriptional start sites and conserved -1 and +1 promoter nucleosomes. Our observations demonstrate that histones and chromatin architecture evolved before the divergence of Archaea and Eukarya, suggesting that the fundamental role of chromatin in the regulation of gene expression is ancient.DOI:http://dx.doi.org/10.7554/eLife.00078.001.

Gene Expression Under the Influence: Transcriptional Profiling of Ethanol in the Brain

Sensitivity to ethanol intoxication, propensity to drink ethanol and vulnerability to develop alcoholism are all influenced by genetic factors. Conversely, exposure to ethanol or subsequent withdrawal produce gene expression changes, which, in combination with environmental variables, may participate in the emergence of compulsive drinking and relapse. The present review offers an integrated perspective on brain gene expression profiling in rodent models of predisposition to differential ethanol sensitivity or consumption, in rats and mice subjected to acute or chronic ethanol exposure, as well as in human alcoholics. The functional categories over-represented among differentially expressed genes suggest that the transcriptional effects of chronic ethanol consumption contribute to the neuroplasticity and neurotoxicity characteristic of alcoholism. Importantly, ethanol produces distinct transcriptional changes within the different brain regions involved in intoxication, reinforcement and addiction. Special emphasis is put on recent profiling studies that have provided some insights into the molecular mechanisms potentially mediating genome-wide regulation of gene expression by ethanol. In particular, current evidence for a role of transcription factors, chromatin remodeling and microRNAs in coordinating the expression of large sets of genes in animals predisposed to excessive ethanol drinking or exposed to protracted abstinence, as well as in human alcoholics, is presented. Finally, studies that have compared ethanol with other drugs of abuse have highlighted common gene expression patterns that may play a central role in drug addiction. The availability of novel technologies and a focus on mechanistic approaches are shaping the future of ethanol transcriptomics.

A new way to start: nanoRNA-mediated priming of transcription initiation

A recent study provides evidence that RNA polymerase uses 2- to ~4-nt RNAs, species termed "nanoRNAs," to prime transcription initiation in Escherichia coli. Priming of transcription initiation with nanoRNAs represents a previously undocumented component of transcription start site selection and gene expression.

Identification of hundreds of novel UPF1 target transcripts by direct determination of whole transcriptome stability

UPF1 eliminates aberrant mRNAs harboring premature termination codons, and regulates the steady-state levels of normal physiological mRNAs. Although genome-wide studies of UPF1 targets performed, previous studies did not distinguish indirect UPF1 targets because they could not determine UPF1-dependent altered RNA stabilities. Here, we measured the decay rates of the whole transcriptome in UPF1-depleted HeLa cells using BRIC-seq, an inhibitor-free method for directly measuring RNA stability. We determined the half-lives and expression levels of 9,229 transcripts. An amount of 785 transcripts were stabilized in UPF1-depleted cells. Among these, the expression levels of 76 transcripts were increased, but those of the other 709 transcripts were not altered. RNA immunoprecipitation showed UPF1 bound to the stabilized transcripts, suggesting that UPF1 directly degrades the 709 transcripts. Many UPF1 targets in this study were newly identified. This study clearly demonstrates that direct determination of RNA stability is a powerful approach for identifying targets of RNA degradation factors.

A next-generation sequencing approach to study the transcriptomic changes during the differentiation of physarum at the single-cell level

Physarum polycephalum is a unicellular eukaryote belonging to the amoebozoa group of organisms. The complex life cycle involves various cell types that differ in morphology, function, and biochemical composition. Sporulation, one step in the life cycle, is a stimulus-controlled differentiation response of macroscopic plasmodial cells that develop into fruiting bodies. Well-established Mendelian genetics and the occurrence of macroscopic cells with a naturally synchronous population of nuclei as source of homogeneous cell material for biochemical analyses make Physarum an attractive model organism for studying the regulatory control of cell differentiation. Here, we develop an approach using RNA-sequencing (RNA-seq), without needing to rely on a genome sequence as a reference, for studying the transcriptomic changes during stimulus-triggered commitment to sporulation in individual plasmodial cells. The approach is validated through the obtained expression patterns and annotations, and particularly the results from up- and downregulated genes, which correlate well with previous studies.

A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis

During the last 3 years, a number of approaches for the normalization of RNA sequencing data have emerged in the literature, differing both in the type of bias adjustment and in the statistical strategy adopted. However, as data continue to accumulate, there has been no clear consensus on the appropriate normalization method to be used or the impact of a chosen method on the downstream analysis. In this work, we focus on a comprehensive comparison of seven recently proposed normalization methods for the differential analysis of RNA-seq data, with an emphasis on the use of varied real and simulated datasets involving different species and experimental designs to represent data characteristics commonly observed in practice. Based on this comparison study, we propose practical recommendations on the appropriate normalization method to be used and its impact on the differential analysis of RNA-seq data.

A quantitative model of transcriptional differentiation driving host-pathogen interactions

Despite our expanding knowledge about the biochemistry of gene regulation involved in host-pathogen interactions, a quantitative understanding of this process at a transcriptional level is still limited. We devise and assess a computational framework that can address this question. This framework is founded on a mixture model-based likelihood, equipped with functionality to cluster genes per dynamic and functional changes of gene expression within an interconnected system composed of the host and pathogen. If genes from the host and pathogen are clustered in the same group due to a similar pattern of dynamic profiles, they are likely to be reciprocally co-evolving. If genes from the two organisms are clustered in different groups, this means that they experience strong host-pathogen interactions. The framework can test the rates of change for individual gene clusters during pathogenic infection and quantify their impacts on host-pathogen interactions. The framework was validated by a pathological study of poplar leaves infected by fungal Marssonina brunnea in which co-evolving and interactive genes that determine poplar-fungus interactions are identified. The new framework should find its wide application to studying host-pathogen interactions for any other interconnected systems.

Gene expression divergence and nucleotide differentiation between males of different color morphs and mating strategies in the ruff

By next generation transcriptome sequencing, it is possible to obtain data on both nucleotide sequence variation and gene expression. We have used this approach (RNA-Seq) to investigate the genetic basis for differences in plumage coloration and mating strategies in a non-model bird species, the ruff (Philomachus pugnax). Ruff males show enormous variation in the coloration of ornamental feathers, used for individual recognition. This polymorphism is linked to reproductive strategies, with dark males (Independents) defending territories on leks against other Independents, whereas white morphs (Satellites) co-occupy Independent's courts without agonistic interactions. Previous work found a strong genetic component for mating strategy, but the genes involved were not identified. We present feather transcriptome data of more than 6,000 de-novo sequenced ruff genes (although with limited coverage for many of them). None of the identified genes showed significant expression divergence between males, but many genetic markers showed nucleotide differentiation between different color morphs and mating strategies. These include several feather keratin genes, splicing factors, and the Xg blood-group gene. Many of the genes with significant genetic structure between mating strategies have not yet been annotated and their functions remain to be elucidated. We also conducted in-depth investigations of 28 pre-identified coloration candidate genes. Two of these (EDNRB and TYR) were specifically expressed in black- and rust-colored males, respectively. We have demonstrated the utility of next generation transcriptome sequencing for identifying and genotyping large number of genetic markers in a non-model species without previous genomic resources, and highlight the potential of this approach for addressing the genetic basis of ecologically important variation.

RNA-Seq Study of Microbially Induced Hemocyte Transcripts from Larval Heliothis virescens (Lepidoptera: Noctuidae)

Larvae of the tobacco budworm are major polyphagous pests throughout the Americas. Development of effective microbial biopesticides for this and related noctuid pests has been stymied by the natural resistance mediated innate immune response. Hemocytes play an early and central role in activating and coordinating immune responses to entomopathogens. To approach this problem we completed RNA-seq expression profiling of hemocytes collected from larvae following an in vivo challenge with bacterial and fungal cell wall components to elicit an immune response. A de novo exome assembly was constructed by combination of sequence tags from all treatments. Sequence tags from each treatment were aligned separately with the assembly to measure expression. The resulting table of differential expression had >22,000 assemblies each with a distinct combination of annotation and expression. Within these assemblies >1,400 were upregulated and >1,500 downregulated by immune activation with bacteria or fungi. Orthologs to innate immune components of other insects were identified including pattern recognition, signal transduction pathways, antimicrobial peptides and enzymes, melanization and coagulation. Additionally orthologs of components regulating hemocytic functions such as autophagy, apoptosis, phagocytosis and nodulation were identified. Associated cellular oxidative defenses and detoxification responses were identified providing a comprehensive snapshot of the early response to elicitation.

Statistical and Computational Methods for High-Throughput Sequencing Data Analysis of Alternative Splicing

The burgeoning field of high-throughput sequencing significantly improves our ability to understand the complexity of transcriptomes. Alternative splicing, as one of the most important driving forces for transcriptome diversity, can now be studied at an unprecedent resolution. Efficient and powerful computational and statistical methods are in urgent need to facilitate the characterization and quantification of alternative splicing events. Here we discuss methods in splice junction read mapping, and methods in exon-centric or isoform-centric quantification of alternative splicing. In addition, we discuss HITS-CLIP and splicing QTL analyses which are novel high-throughput sequencing based approaches in the dissection of splicing regulation.

Quantification of Transcriptome Responses of the Rumen Epithelium to Butyrate Infusion using RNA-seq Technology

Short-chain fatty acids (SCFAs), such as butyrate, produced by gut microorganisms, play a critical role in energy metabolism and physiology of ruminants as well as in human health. In this study, the temporal effect of elevated butyrate concentrations on the transcriptome of the rumen epithelium was quantified via serial biopsy sampling using RNA-seq technology. The mean number of genes transcribed in the rumen epithelial transcriptome was 17,323.63 ± 277.20 (±SD; N = 24) while the core transcriptome consisted of 15,025 genes. Collectively, 80 genes were identified as being significantly impacted by butyrate infusion across all time points sampled. Maximal transcriptional effect of butyrate on the rumen epithelium was observed at the 72-h infusion when the abundance of 58 genes was altered. The initial reaction of the rumen epithelium to elevated exogenous butyrate may represent a stress response as Gene Ontology (GO) terms identified were predominantly related to responses to bacteria and biotic stimuli. An algorithm for the reconstruction of accurate cellular networks (ARACNE) inferred regulatory gene networks with 113,738 direct interactions in the butyrate-epithelium interactome using a combined cutoff of an error tolerance (ɛ = 0.10) and a stringent P-value threshold of mutual information (5.0 × 10⁻¹¹). Several regulatory networks were controlled by transcription factors, such as CREBBP and TTF2, which were regulated by butyrate. Our findings provide insight into the regulation of butyrate transport and metabolism in the rumen epithelium, which will guide our future efforts in exploiting potential beneficial effect of butyrate in animal well-being and human health.

RNA-sequencing analysis of 5' capped RNAs identifies many new differentially expressed genes in acute hepatitis C virus infection

We describe the first report of RNA sequencing of 5' capped (Pol II) RNAs isolated from acutely hepatitis C virus (HCV) infected Huh 7.5 cells that provides a general approach to identifying differentially expressed annotated and unannotated genes that participate in viral-host interactions. We identified 100, 684, and 1,844 significantly differentially expressed annotated genes in acutely infected proliferative Huh 7.5 cells at 6, 48, and 72 hours, respectively (fold change ≥ 1.5 and Bonferroni adjusted p-values < 0.05). Most of the differentially expressed genes (>80%) and biological pathways (such as adipocytokine, Notch, Hedgehog and NOD-like receptor signaling) were not identified by previous gene array studies. These genes are critical components of host immune, inflammatory and oncogenic pathways and provide new information regarding changes that may benefit the virus or mediate HCV induced pathology. RNAi knockdown studies of newly identified highly upregulated FUT1 and KLHDC7B genes provide evidence that their gene products regulate and facilitate HCV replication in hepatocytes. Our approach also identified novel Pol II unannotated transcripts that were upregulated. Results further identify new pathways that regulate HCV replication in hepatocytes and suggest that our approach will have general applications in studying viral-host interactions in model systems and clinical biospecimens.

An Interspecies Regulatory Network Inferred from Simultaneous RNA-seq of Candida albicans Invading Innate Immune Cells

The ability to adapt to diverse micro-environmental challenges encountered within a host is of pivotal importance to the opportunistic fungal pathogen Candida albicans. We have quantified C. albicans and M. musculus gene expression dynamics during phagocytosis by dendritic cells in a genome-wide, time-resolved analysis using simultaneous RNA-seq. A robust network inference map was generated from this dataset using NetGenerator, predicting novel interactions between the host and the pathogen. We experimentally verified predicted interdependent sub-networks comprising Hap3 in C. albicans, and Ptx3 and Mta2 in M. musculus. Remarkably, binding of recombinant Ptx3 to the C. albicans cell wall was found to regulate the expression of fungal Hap3 target genes as predicted by the network inference model. Pre-incubation of C. albicans with recombinant Ptx3 significantly altered the expression of Mta2 target cytokines such as IL-2 and IL-4 in a Hap3-dependent manner, further suggesting a role for Mta2 in host-pathogen interplay as predicted in the network inference model. We propose an integrated model for the functionality of these sub-networks during fungal invasion of immune cells, according to which binding of Ptx3 to the C. albicans cell wall induces remodeling via fungal Hap3 target genes, thereby altering the immune response to the pathogen. We show the applicability of network inference to predict interactions between host-pathogen pairs, demonstrating the usefulness of this systems biology approach to decipher mechanisms of microbial pathogenesis.

Transcriptome complexity and riboregulation in the human pathogen Helicobacter pylori

The Gram-negative Epsilonproteobacterium Helicobacter pylori is considered as one of the major human pathogens and many studies have focused on its virulence mechanisms as well as genomic diversity. In contrast, only very little is known about post-transcriptional regulation and small regulatory RNAs (sRNAs) in this spiral-shaped microaerophilic bacterium. Considering the absence of the common RNA chaperone Hfq, which is a key-player in post-transcriptional regulation in enterobacteria, H. pylori was even regarded as an organism without riboregulation. However, analysis of the H. pylori primary transcriptome using RNA-seq revealed a very complex transcriptional output from its small genome. Furthermore, the identification of a wealth of sRNAs as well as massive antisense transcription indicates that H. pylori uses riboregulation for its gene expression control. The ongoing functional characterization of sRNAs along with the identification of associated RNA binding proteins will help to understand their potential roles in Helicobacter virulence and stress response. Moreover, research on riboregulation in H. pylori will provide new insights into its virulence mechanisms and will also help to shed light on post-transcriptional regulation in other Epsilonproteobacteria, including widespread and emerging pathogens such as Campylobacter.

Strain Variation in the Transcriptome of the Dengue Fever Vector, Aedes aegypti

Studies of transcriptome dynamics provide a basis for understanding functional elements of the genome and the complexity of gene regulation. The dengue vector mosquito, Aedes aegypti, exhibits great adaptability to diverse ecological conditions, is phenotypically polymorphic, and shows variation in vectorial capacity to arboviruses. Previous genome sequencing showed richness in repetitive DNA and transposable elements that can contribute to genome plasticity. Population genetic studies revealed a varying degree of worldwide genetic polymorphism. However, the extent of functional genetic polymorphism across strains is unknown. The transcriptomes of three Ae. aegypti strains, Chetumal (CTM), Rexville D-Puerto Rico (Rex-D) and Liverpool (LVP), were compared. CTM is more susceptible than Rex- D to infection by dengue virus serotype 2. A total of 4188 transcripts exhibit either no or small variation (<2-fold) among sugar-fed samples of the three strains and between sugar- and blood-fed samples within each strain, corresponding most likely to genes encoding products necessary for vital functions. Transcripts enriched in blood-fed mosquitoes encode proteins associated with catalytic activities, molecular transport, metabolism of lipids, carbohydrates and amino acids, and functions related to blood digestion and the progression of the gonotropic cycle. Significant qualitative and quantitative differences were found in individual transcripts among strains including differential representation of paralogous gene products. The majority of immunity-associated transcripts decreased in accumulation after a bloodmeal and the results are discussed in relation to the different susceptibility of CTM and Rex-D mosquitoes to DENV2 infection.

RNA-Seq Assembly - Are We There Yet?

Transcriptomic sequence resources represent invaluable assets for research, in particular for non-model species without a sequenced genome. To date, the Next Generation Sequencing technologies 454/Roche and Illumina have been used to generate transcriptome sequence databases by mRNA-Seq for more than fifty different plant species. While some of the databases were successfully used for downstream applications, such as proteomics, the assembly parameters indicate that the assemblies do not yet accurately reflect the actual plant transcriptomes. Two different assembly strategies have been used, overlap consensus based assemblers for long reads and Eulerian path/de Bruijn graph assembler for short reads. In this review, we discuss the challenges and solutions to the transcriptome assembly problem. A list of quality control parameters and the necessary scripts to produce them are provided.

Next generation quantitative genetics in plants

Most characteristics in living organisms show continuous variation, which suggests that they are controlled by multiple genes. Quantitative trait loci (QTL) analysis can identify the genes underlying continuous traits by establishing associations between genetic markers and observed phenotypic variation in a segregating population. The new high-throughput sequencing (HTS) technologies greatly facilitate QTL analysis by providing genetic markers at genome-wide resolution in any species without previous knowledge of its genome. In addition HTS serves to quantify molecular phenotypes, which aids to identify the loci responsible for QTLs and to understand the mechanisms underlying diversity. The constant improvements in price, experimental protocols, computational pipelines, and statistical frameworks are making feasible the use of HTS for any research group interested in quantitative genetics. In this review I discuss the application of HTS for molecular marker discovery, population genotyping, and expression profiling in QTL analysis.

The decade of the epigenomes?

The beginning of this century was not only marked by the publication of the first draft of the human genome but also set off a decade of intense research on epigenetic phenomena. Apart from DNA methylation, it became clear that many other factors including a wide range of histone modifications, different shades of chromatin accessibility, and a vast suite of noncoding RNAs comprise the epigenome. With the recent advances in sequencing technologies, it has now become possible to analyze many of these features in depth, allowing for the first time the establishment of complete epigenomic profiles for basically every cell type of interest. Here, we will discuss the recent advances that allow comprehensive epigenetic mapping, highlight several projects that set out to better understand the epigenome, and discuss the impact that epigenomic mapping can have on our understanding of both healthy and diseased cells.

Baculovirus induced transcripts in hemocytes from the larvae of Heliothis virescens

Using RNA-seq digital difference expression profiling methods, we have assessed the gene expression profiles of hemocytes harvested from Heliothis virescens that were challenged with Helicoverpa zea single nucleopolyhedrovirus (HzSNPV). A reference transcriptome of hemocyte-expressed transcripts was assembled from 202 million 42-base tags by combining the sequence data of all samples, and the assembled sequences were then subject to BLASTx analysis to determine gene identities. We used the fully sequenced HzSNPV reference genome to align 477,264 Illumina sequence tags from infected hemocytes in order to document expression of HzSNPV genes at early points during infection. A comparison of expression profiles of control insects to those lethally infected with HzSNPV revealed differential expression of key cellular stress response genes and genes involved in lipid metabolism. Transcriptional regulation of specific insect hormones in baculovirus-infected insects was also altered. A number of transcripts bearing homology to retroviral elements that were detected add to a growing body of evidence for extensive invasion of errantiviruses into the insect genome. Using this method, we completed the first and most comprehensive gene expression survey of both baculoviral infection and host immune defense in lepidopteran larvae.

DEB: A web interface for RNA-seq digital gene expression analysis

Digital expression (DE) is an important application of RNA-seq technology to quantify the transcriptome. The number of mapped reads to each transcript or gene varies under different conditions and replicates. Currently, three different statistical algorithms (edgeR, DESeq and bayseq) are available as R packages, to compare the reads to identify significantly expressed transcripts or genes. So far, users have to manually install and run each R package separately. It is also of users' interest to compare the results of different approaches. Here, we present a pipeline DEB which automates all the steps in file preparation, computation and result comparison.

AVAILABILITY

The database is available for free at http://www.ijbcb.org/DEB/php/onlinetool.php.