Search

Background Ribosomal RNA (rRNA) comprises at least 90% of total RNA extracted from mammalian tissue or cell line samples. Informative transcriptional profiling using massively parallel sequencing technologies requires either enrichment of mature poly-adenylated transcripts or targeted depletion of the rRNA fraction. The latter method is of particular interest because it is compatible with degraded samples such as those extracted from FFPE and also captures transcripts that are not poly-adenylated such as some non-coding RNAs. Here we provide a cross-site study that evaluates the performance of ribosomal RNA removal kits from Illumina, Takara/Clontech, Kapa Biosystems, Lexogen, New England Biolabs and Qiagen on intact and degraded RNA samples. Results We find that all of the kits are capable of performing significant ribosomal depletion, though there are differences in their ease of use. All kits were able to remove ribosomal RNA to below 20% with intact RNA and identify ~ 14,000 protein coding genes from the Universal Human Reference RNA sample at >1FPKM. Analysis of differentially detected genes between kits suggests that transcript length may be a key factor in library production efficiency. Conclusions These results provide a roadmap for labs on the strengths of each of these methods and how best to utilize them. Keywords: RNAseqr; RNA depletion; Illumina; NGS; ABRF; Transcriptomics, National Cancer Institute (U.S.) (Grant P30-CA14051), National Institute of Environmental Health Sciences (Grant P30-ES002109)

Efficient lignin depolymerization is unique to the wood decay basidiomycetes, collectively referred to as white rot fungi. Phanerochaete chrysosporium simultaneously degrades lignin and cellulose, whereas the closely related species, Ceriporiopsis subvermispora, also depolymerizes lignin but may do so with relatively little cellulose degradation. To investigate the basis for selective ligninolysis, we conducted comparative genome analysis of C. subvermispora and P. chrysosporium . Genes encoding manganese peroxidase numbered 13 and five in C. subvermispora and P. chrysosporium , respectively. In addition, the C. subvermispora genome contains at least seven genes predicted to encode laccases, whereas the P. chrysosporium genome contains none. We also observed expansion of the number of C. subvermispora desaturase-encoding genes putatively involved in lipid metabolism. Microarray-based transcriptome analysis showed substantial up-regulation of several desaturase and MnP genes in wood-containing medium. MS identified MnP proteins in C. subvermispora culture filtrates, but none in P. chrysosporium cultures. These results support the importance of MnP and a lignin degradation mechanism whereby cleavage of the dominant nonphenolic structures is mediated by lipid peroxidation products. Two C. subvermispora genes were predicted to encode peroxidases structurally similar to P. chrysosporium lignin peroxidase and, following heterologous expression in Escherichia coli , the enzymes were shown to oxidize high redox potential substrates, but not Mn 2+ . Apart from oxidative lignin degradation, we also examined cellulolytic and hemicellulolytic systems in both fungi. In summary, the C. subvermispora genetic inventory and expression patterns exhibit increased oxidoreductase potential and diminished cellulolytic capability relative to P. chrysosporium .

Epstein-Barr virus (EBV) reactivation involves the ordered induction of approximately 90 viral genes that participate in the generation of infectious virions. Using strand-specific RNA-seq to assess the EBV transcriptome during reactivation, we found extensive bidirectional transcription extending across nearly the entire genome. In contrast, only 4% of the EBV genome is currently bidirectionally annotated. Most of the newly identified transcribed regions show little evidence of coding potential, supporting noncoding roles for most of these RNAs. Based on previous cellular long noncoding RNA size calculations, we estimate that there are likely hundreds more EBV genes expressed during reactivation than was previously known. Limited 5′ and 3′ rapid amplification of cDNA ends (RACE) experiments and findings of novel splicing events by RNA-seq suggest that the complexity of the viral genome during reactivation may be even greater. Further analysis of antisense transcripts at some of the EBV latency gene loci showed that they are “late” genes, they are nuclear, and they tend to localize in areas of the nucleus where others find newly synthesized viral genomes. This raises the possibility that these transcripts perform functions such as new genome processing, stabilization, organization, etc. The finding of a significantly more complex EBV transcriptome during reactivation changes our view of the viral production process from one that is facilitated and regulated almost entirely by previously identified viral proteins to a process that also involves the contribution of a wide array of virus encoded noncoding RNAs. IMPORTANCE Epstein-Barr virus (EBV) is a herpesvirus that infects the majority of the world's population, in rare cases causing serious disease such as lymphoma and gastric carcinoma. Using strand-specific RNA-seq, we have studied viral gene expression during EBV reactivation and have discovered hundreds more viral transcripts than were previously known. The finding of alternative splicing and the prevalence of overlapping transcripts indicate additional complexity. Most newly identified transcribed regions do not encode proteins but instead likely function as noncoding RNA molecules which could participate in regulating gene expression, gene splicing or even activities such as viral genome processing. These findings broaden the scope of what we need to consider to understand the viral manufacturing process. As more detailed studies are undertaken they will likely change the way we view this process as a whole.

Certain wood decay basidiomycetes, collectively referred to as brown rot fungi, rapidly depolymerize cellulose while leaving behind the bulk of cell wall lignin as a modified residue. The mechanism(s) employed is unclear, but considerable evidence implicates the involvement of diffusible oxidants generated via Fenton-like chemistry. Toward a better understanding of this process, we have examined the transcriptome and secretome of Wolfiporia cocos when cultivated on media containing glucose, purified crystalline cellulose, aspen ( Populus grandidentata ), or lodgepole pine ( Pinus contorta ) as the sole carbon source. Compared to the results obtained with glucose, 30, 183, and 207 genes exhibited 4-fold increases in transcript levels in cellulose, aspen, and lodgepole pine, respectively. Mass spectrometry identified peptides corresponding to 64 glycoside hydrolase (GH) proteins, and of these, 17 corresponded to transcripts upregulated on one or both woody substrates. Most of these genes were broadly categorized as hemicellulases or chitinases. Consistent with an important role for hydroxyl radical in cellulose depolymerization, high transcript levels and upregulation were observed for genes involved in iron homeostasis, iron reduction, and extracellular peroxide generation. These patterns of regulation differ markedly from those of the closely related brown rot fungus Postia placenta and expand the number of enzymes potentially involved in the oxidative depolymerization of cellulose. IMPORTANCE The decomposition of wood is an essential component of nutrient cycling in forest ecosystems. Few microbes have the capacity to efficiently degrade woody substrates, and the mechanism(s) is poorly understood. Toward a better understanding of these processes, we show that when grown on wood as a sole carbon source the brown rot fungus W. cocos expresses a unique repertoire of genes involved in oxidative and hydrolytic conversions of cell walls.

Identification of specific genes and enzymes involved in conversion of lignocellulosics from an expanding number of potential feedstocks is of growing interest to bioenergy process development. The basidiomycetous wood decay fungi Phanerochaete chrysosporium and Postia placenta are promising in this regard because they are able to utilize a wide range of simple and complex carbon compounds. However, systematic comparative studies with different woody substrates have not been reported. To address this issue, we examined gene expression of these fungi colonizing aspen ( Populus grandidentata ) and pine ( Pinus strobus ). Transcript levels of genes encoding extracellular glycoside hydrolases, thought to be important for hydrolytic cleavage of hemicelluloses and cellulose, showed little difference for P. placenta colonizing pine versus aspen as the sole carbon source. However, 164 genes exhibited significant differences in transcript accumulation for these substrates. Among these, 15 cytochrome P450s were upregulated in pine relative to aspen. Of 72 P. placenta extracellular proteins identified unambiguously by mass spectrometry, 52 were detected while colonizing both substrates and 10 were identified in pine but not aspen cultures. Most of the 178 P. chrysosporium glycoside hydrolase genes showed similar transcript levels on both substrates, but 13 accumulated >2-fold higher levels on aspen than on pine. Of 118 confidently identified proteins, 31 were identified in both substrates and 57 were identified in pine but not aspen cultures. Thus, P. placenta and P. chrysosporium gene expression patterns are influenced substantially by wood species. Such adaptations to the carbon source may also reflect fundamental differences in the mechanisms by which these fungi attack plant cell walls.

We have developed and optimized the necessary laboratory materials to make DNA microarray technology accessible to all high school students at a fraction of both cost and data size. The primary component is a DNA chip/array that students "print" by hand and then analyze using research tools that have been adapted for classroom use. The primary adaptation is the use of a simulated cDNA target. The low density DNA array we discuss here was used to demonstrate differential expression of several Arabidopsis thaliana genes related to photosynthesis and photomorphogenesis. The methods we present here can be used with any biological organism whose sequence is known. Furthermore, these methods can be adapted to exhibit a variety of differential gene expression patterns under different experimental conditions. The materials and tools we discuss have been applied in classrooms at West High School in Madison, WI. We have also shared these materials with high school teachers attending professional development courses at the University of Wisconsin-Madison.

20 p.-3 tab.-9 fig., Collectively classified as white-rot fungi, certain basidiomycetes efficiently degrade the major structural polymers of wood cell walls. A small subset of these Agaricomycetes, exemplified by Phlebiopsis gigantea, is capable of colonizing freshly exposed conifer sapwood despite its high content of extractives, which retards the establishment of other fungal species. The mechanism(s) by which P. gigantea tolerates and metabolizes resinous compounds have not been explored. Here, we report the annotated P. gigantea genome and compare profiles of its transcriptome and secretome when cultured on freshcut versus solvent-extracted loblolly pine wood. The P. gigantea genome contains a conventional repertoire of hydrolase genes involved in cellulose/hemicellulose degradation, whose patterns of expression were relatively unperturbed by the absence of extractives. The expression of genes typically ascribed to lignin degradation was also largely unaffected. In contrast, genes likely involved in the transformation and detoxification of wood extractives were highly induced in its presence. Their products included an ABC transporter, lipases, cytochrome P450s, glutathione S-transferase and aldehyde dehydrogenase. Other regulated genes of unknown function and several constitutively expressed genes are also likely involved in P. gigantea’s extractives metabolism. These results contribute to our fundamental understanding of pioneer colonization of conifer wood and provide insight into the diverse chemistries employed by fungi in carbon cycling processes., The major portions of this work were performed under US Department of Agriculture Cooperative State, Research, Education, and Extension Service Grant 2007-35504-18257 (to DC and RAB). The US Department of Energy Joint Genome Institute is supported by the Office of Science of the US Department of Energy under Contract DE-AC02-05CH11231. This work was also supported by the HIPOP (BIO2011-26694) project of the Spanish Ministry of Economy and Competitiveness (MINECO) (to FJRD), the PEROXICATS (KBBE-2010-4-265397) and INDOX (KBBE-2013-.3.3-04-613549) European projects (to ATM), and the Chilean National Fund for Scientific and Technological Development Grant 1131030 (to LFL).

Gene expression profiles of Escherichia coli K-12 W3110 were compared as a function of steady-state external pH. Cultures were grown to an optical density at 600 nm of 0.3 in potassium-modified Luria-Bertani medium buffered at pH 5.0, 7.0, and 8.7. For each of the three pH conditions, cDNA from RNA of five independent cultures was hybridized to Affymetrix E . coli arrays. Analysis of variance with an α level of 0.001 resulted in 98% power to detect genes showing a twofold difference in expression. Normalized expression indices were calculated for each gene and intergenic region (IG). Differential expression among the three pH classes was observed for 763 genes and 353 IGs. Hierarchical clustering yielded six well-defined clusters of pH profiles, designated Acid High (highest expression at pH 5.0), Acid Low (lowest expression at pH 5.0), Base High (highest at pH 8.7), Base Low (lowest at pH 8.7), Neutral High (highest at pH 7.0, lower in acid or base), and Neutral Low (lowest at pH 7.0, higher at both pH extremes). Flagellar and chemotaxis genes were repressed at pH 8.7 (Base Low cluster), where the cell's transmembrane proton potential is diminished by the maintenance of an inverted pH gradient. High pH also repressed the proton pumps cytochrome o ( cyo ) and NADH dehydrogenases I and II. By contrast, the proton-importing ATP synthase F 1 F o and the microaerophilic cytochrome d ( cyd ), which minimizes proton export, were induced at pH 8.7. These observations are consistent with a model in which high pH represses synthesis of flagella, which expend proton motive force, while stepping up electron transport and ATPase components that keep protons inside the cell. Acid-induced genes, on the other hand, were coinduced by conditions associated with increased metabolic rate, such as oxidative stress. All six pH-dependent clusters included envelope and periplasmic proteins, which directly experience external pH. Overall, this study showed that (i) low pH accelerates acid consumption and proton export, while coinducing oxidative stress and heat shock regulons; (ii) high pH accelerates proton import, while repressing the energy-expensive flagellar and chemotaxis regulons; and (iii) pH differentially regulates a large number of periplasmic and envelope proteins.

The high level of accuracy and sensitivity of next generation sequencing for quantifying genetic material across organismal boundaries gives it tremendous potential for pathogen discovery and diagnosis in human disease. Despite this promise, substantial bacterial contamination is routinely found in existing human-derived RNA-seq datasets that likely arises from environmental sources. This raises the need for stringent sequencing and analysis protocols for studies investigating sequence-based microbial signatures in clinical samples.

We examined gene expression patterns in the lignin-degrading fungus Phanerochaete chrysosporium when it colonizes hybrid poplar ( Populus alba × tremula ) and syringyl (S)-rich transgenic derivatives. A combination of microarrays and liquid chromatography-tandem mass spectrometry (LC-MS/MS) allowed detection of a total of 9,959 transcripts and 793 proteins. Comparisons of P. chrysosporium transcript abundance in medium containing poplar or glucose as a sole carbon source showed 113 regulated genes, 11 of which were significantly higher (>2-fold, P < 0.05) in transgenic line 64 relative to the parental line. Possibly related to the very large amounts of syringyl (S) units in this transgenic tree (94 mol% S), several oxidoreductases were among the upregulated genes. Peptides corresponding to a total of 18 oxidoreductases were identified in medium consisting of biomass from line 64 or 82 (85 mol% S) but not in the parental clone (65 mol% S). These results demonstrate that P. chrysosporium gene expression patterns are substantially influenced by lignin composition.

Corrigendum: Insights into the role of DNA methylation in diatoms by genome-wide profiling in Phaeodactylum tricornutum

DNA cytosine methylation is a widely conserved epigenetic mark in eukaryotes that appears to have critical roles in the regulation of genome structure and transcription. Genome-wide methylation maps have so far only been established from the supergroups Archaeplastida and Unikont. Here we report the first whole-genome methylome from a stramenopile, the marine model diatom Phaeodactylum tricornutum. Around 6% of the genome is intermittently methylated in a mosaic pattern. We find extensive methylation in transposable elements. We also detect methylation in over 320 genes. Extensive gene methylation correlates strongly with transcriptional silencing and differential expression under specific conditions. By contrast, we find that genes with partial methylation tend to be constitutively expressed. These patterns contrast with those found previously in other eukaryotes. By going beyond plants, animals and fungi, this stramenopile methylome adds significantly to our understanding of the evolution of DNA methylation in eukaryotes. Fil: Veluchamy, Alaguraj. Institut de Biologie de l'École Normale Supérieure; Francia Fil: Lin, Xin. Institut de Biologie de l'École Normale Supérieure; Francia. Xiamen University; China Fil: Maumus, Florian. Fil: Rivarola, Maximo Lisandro. Fil: Bhavsar, Jaysheel. Fil: Creasy, Todd. Fil: O'Brien, Kimberly. Fil: Sengamalay, Naomi A.. Fil: Tallon, Luke J.. Fil: Smith, Andrew D.. Fil: Rayko, Edda. Fil: Ahmed, Ikhlak. Fil: Crom, Stéphane Le. Fil: Farrant, Gregory K.. Fil: Sgro, Jean-Yves. Fil: Olson, Sue A.. Fil: Bondurant, Sandra Splinter. Fil: Allen, Andrew. Fil: Rabinowicz, Pablo D.. Fil: Sussman, Michael R.. Fil: Bowler, Chris. Fil: Tirichine, Leïla.

Mechanisms for controlling symbiont populations are critical for maintaining the associations that exist between a host and its microbial partners. We describe here the transcriptional, metabolic, and ultrastructural characteristics of a diel rhythm that occurs in the symbiosis between the squid Euprymna scolopes and the luminous bacterium Vibrio fischeri. The rhythm is driven by the host’s expulsion from its light-emitting organ of most of the symbiont population each day at dawn. The transcriptomes of both the host epithelium that supports the symbionts and the symbiont population itself were characterized and compared at four times over this daily cycle. The greatest fluctuation in gene expression of both partners occurred as the day began. Most notable was an up-regulation in the host of >50 cytoskeleton-related genes just before dawn and their subsequent down-regulation within 6 h. Examination of the epithelium by TEM revealed a corresponding restructuring, characterized by effacement and blebbing of its apical surface. After the dawn expulsion, the epithelium reestablished its polarity, and the residual symbionts began growing, repopulating the light organ. Analysis of the symbiont transcriptome suggested that the bacteria respond to the effacement by up-regulating genes associated with anaerobic respiration of glycerol; supporting this finding, lipid analysis of the symbionts’ membranes indicated a direct incorporation of host-derived fatty acids. After 12 h, the metabolic signature of the symbiont population shifted to one characteristic of chitin fermentation, which continued until the following dawn. Thus, the persistent maintenance of the squid–vibrio symbiosis is tied to a dynamic diel rhythm that involves both partners.

Brown-rot fungi such as Postia placenta are common inhabitants of forest ecosystems and are also largely responsible for the destructive decay of wooden structures. Rapid depolymerization of cellulose is a distinguishing feature of brown-rot, but the biochemical mechanisms and underlying genetics are poorly understood. Systematic examination of the P. placenta genome, transcriptome, and secretome revealed unique extracellular enzyme systems, including an unusual repertoire of extracellular glycoside hydrolases. Genes encoding exocellobiohydrolases and cellulose-binding domains, typical of cellulolytic microbes, are absent in this efficient cellulose-degrading fungus. When P. placenta was grown in medium containing cellulose as sole carbon source, transcripts corresponding to many hemicellulases and to a single putative β -1–4 endoglucanase were expressed at high levels relative to glucose-grown cultures. These transcript profiles were confirmed by direct identification of peptides by liquid chromatography-tandem mass spectrometry (LC-MS/MS). Also upregulated during growth on cellulose medium were putative iron reductases, quinone reductase, and structurally divergent oxidases potentially involved in extracellular generation of Fe(II) and H2O2. These observations are consistent with a biodegradative role for Fenton chemistry in which Fe(II) and H2O2 react to form hydroxyl radicals, highly reactive oxidants capable of depolymerizing cellulose. The P. placenta genome resources provide unparalleled opportunities for investigating such unusual mechanisms of cellulose conversion. More broadly, the genome offers insight into the diversification of lignocellulose degrading mechanisms in fungi. Comparisons with the closely related white-rot fungus Phanerochaete chrysosporium support an evolutionary shift from white-rot to brown-rot during which the capacity for efficient depolymerization of lignin was lost. This work was supported by the U.S. Department of Energy’s Office of Science, Biological and Environmental Research Program, and University of California, Lawrence Berkeley National Laboratory Contract DE-AC02–05CH11231; Lawrence Livermore National Laboratory Contract DE-AC52–07NA27344; Los Alamos National Laboratory Contract DE-AC02–06NA25396; University of Wisconsin Grant DE-FG02–87ER13712; Forest Products Laboratory, U.S. Department of Agriculture, Cooperative State Research, Education, and Extension Services Grant 2007–35504-18257; National Institutes of Health Grant GM060201 (to University of New Mexico); Centro de Investigaciones Biológicas (Madrid) EUproject NMP2–2006-026456; Ministry of Education Czech Republic Grant LC06066.

Background Many E. coli genes show pH-dependent expression during logarithmic growth in acid (pH 5–6) or in base (pH 8–9). The effect of rapid pH change, however, has rarely been tested. Rapid acid treatment could distinguish between genes responding to external pH, and genes responding to cytoplasmic acidification, which occurs transiently following rapid external acidification. It could reveal previously unknown acid-stress genes whose effects are transient, as well as show which acid-stress genes have a delayed response. Results Microarray hybridization was employed to observe the global gene expression of E. coli K-12 W3110 following rapid acidification of the external medium, from pH 7.6 to pH 5.5. Fluorimetric observation of pH-dependent tetR-YFP showed that rapid external acidification led to a half-unit drop in cytoplasmic pH (from pH 7.6 to pH 6.4) which began to recover within 20 s. Following acid treatment, 630 genes were up-regulated and 586 genes were down-regulated. Up-regulated genes included amino-acid decarboxylases (cadA, adiY, gadA), succinate dehydrogenase (sdhABCD), biofilm-associated genes (bdm, gatAB, and ymgABC), and the Gad, Fur and Rcs regulons. Genes with response patterns consistent with cytoplasmic acid stress were revealed by addition of benzoate, a membrane-permeant acid that permanently depresses cytoplasmic pH without affecting external pH. Several genes (yagU, ygiN, yjeI, and yneI) were up-regulated specifically by external acidification, while other genes (fimB, ygaC, yhcN, yhjX, ymgABC, yodA) presented a benzoate response consistent with cytoplasmic pH stress. Other genes (the nuo operon for NADH dehydrogenase I, and the HslUV protease) showed delayed up-regulation by acid, with expression rising by 10 min following the acid shift. Conclusion Transcriptomic profiling of E. coli K-12 distinguished three different classes of change in gene expression following rapid acid treatment: up-regulation with or without recovery, and delayed response to acid. For eight genes showing acid response and recovery (fimB, ygaC, yhcN, yhjX, ymgABC, yodA), responses to the permeant acid benzoate revealed expression patterns consistent with sensing of cytoplasmic pH. The delayed acid response of nuo genes shows that NADH dehydrogenase I is probably induced as a secondary result of acid-associated metabolism, not as a direct response to cytoplasmic acidification.

Background In Escherichia coli, pH regulates genes for amino-acid and sugar catabolism, electron transport, oxidative stress, periplasmic and envelope proteins. Many pH-dependent genes are co-regulated by anaerobiosis, but the overall intersection of pH stress and oxygen limitation has not been investigated. Results The pH dependence of gene expression was analyzed in oxygen-limited cultures of E. coli K-12 strain W3110. E. coli K-12 strain W3110 was cultured in closed tubes containing LBK broth buffered at pH 5.7, pH 7.0, and pH 8.5. Affymetrix array hybridization revealed pH-dependent expression of 1,384 genes and 610 intergenic regions. A core group of 251 genes showed pH responses similar to those in a previous study of cultures grown with aeration. The highly acid-induced gene yagU was shown to be required for extreme-acid resistance (survival at pH 2). Acid also up-regulated fimbriae (fimAC), periplasmic chaperones (hdeAB), cyclopropane fatty acid synthase (cfa), and the "constitutive" Na+/H+ antiporter (nhaB). Base up-regulated core genes for maltodextrin transport (lamB, mal), ATP synthase (atp), and DNA repair (recA, mutL). Other genes showed opposite pH responses with or without aeration, for example ETS components (cyo,nuo, sdh) and hydrogenases (hya, hyb, hyc, hyf, hyp). A hypF strain lacking all hydrogenase activity showed loss of extreme-acid resistance. Under oxygen limitation only, acid down-regulated ribosome synthesis (rpl,rpm, rps). Acid up-regulated the catabolism of sugar derivatives whose fermentation minimized acid production (gnd, gnt, srl), and also a cluster of 13 genes in the gadA region. Acid up-regulated drug transporters (mdtEF, mdtL), but down-regulated penicillin-binding proteins (dacACD, mreBC). Intergenic regions containing regulatory sRNAs were up-regulated by acid (ryeA, csrB, gadY, rybC). Conclusion pH regulates a core set of genes independently of oxygen, including yagU, fimbriae, periplasmic chaperones, and nhaB. Under oxygen limitation, however, pH regulation is reversed for genes encoding electron transport components and hydrogenases. Extreme-acid resistance requires yagU and hydrogenase production. Ribosome synthesis is down-regulated at low pH under oxygen limitation, possibly due to the restricted energy yield of catabolism. Under oxygen limitation, pH regulates metabolism and transport so as to maximize alternative catabolic options while minimizing acidification or alkalinization of the cytoplasm.