Global amplification of mRNA by template-switching PCR: linearity and application to microarray analysis

Revealing protein-protein interactions at the transcriptome scale by sequencing

We describe PROPER-seq (protein-protein interaction sequencing) to map protein-protein interactions (PPIs) en masse. PROPER-seq first converts transcriptomes of input cells into RNA-barcoded protein libraries, in which all interacting protein pairs are captured through nucleotide barcode ligation, recorded as chimeric DNA sequences, and decoded at once by sequencing and mapping. We applied PROPER-seq to human embryonic kidney cells, T lymphocytes, and endothelial cells and identified 210,518 human PPIs (collected in the PROPER v.1.0 database). Among these, 1,365 and 2,480 PPIs are supported by published co-immunoprecipitation (coIP) and affinity purification-mass spectrometry (AP-MS) data, 17,638 PPIs are predicted by the prePPI algorithm without previous experimental validation, and 100 PPIs overlap human synthetic lethal gene pairs. In addition, four previously uncharacterized interaction partners with poly(ADP-ribose) polymerase 1 (PARP1) (a critical protein in DNA repair) known as XPO1, MATR3, IPO5, and LEO1 are validated in vivo. PROPER-seq presents a time-effective technology to map PPIs at the transcriptome scale, and PROPER v.1.0 provides a rich resource for studying PPIs.

Chemical capping improves template switching and enhances sequencing of small RNAs

Template-switching reverse transcription is widely used in RNA sequencing for low-input and low-quality samples, including RNA from single cells or formalin-fixed paraffin-embedded (FFPE) tissues. Previously, we identified the native eukaryotic mRNA 5′ cap as a key structural element for enhancing template switching efficiency. Here, we introduce CapTS-seq, a new strategy for sequencing small RNAs that combines chemical capping and template switching. We probed a variety of non-native synthetic cap structures and found that an unmethylated guanosine triphosphate cap led to the lowest bias and highest efficiency for template switching. Through cross-examination of different nucleotides at the cap position, our data provided unequivocal evidence that the 5′ cap acts as a template for the first nucleotide in reverse transcriptase-mediated post-templated addition to the emerging cDNA—a key feature to propel template switching. We deployed CapTS-seq for sequencing synthetic miRNAs, human total brain and liver FFPE RNA, and demonstrated that it consistently improves library quality for miRNAs in comparison with a gold standard template switching-based small RNA-seq kit.

Non-templated addition and template switching by Moloney murine leukemia virus (MMLV)-based reverse transcriptases co-occur and compete with each other

Single-cell RNA-Seq (scRNA-Seq) has led to an unprecedented understanding of gene expression and regulation in individual cells. Many scRNA-Seq approaches rely upon the template switching property of Moloney murine leukemia virus (MMLV)-type reverse transcriptases. Template switching is believed to happen in a sequential process involving nontemplated addition of three protruding nucleotides (+CCC) to the 3′-end of the nascent cDNA, which can then anneal to the matching rGrGrG 3′-end of the template-switching oligo (TSO), allowing the reverse transcriptase (RT) to switch templates and continue copying the TSO sequence. In this study, we present a detailed analysis of template switching biases with respect to the RNA template, specifically of the role of the sequence and nature of its 5′-end (capped versus noncapped) in these biases. Our findings confirmed that the presence of a 5′-m⁷G cap enhances template switching efficiency. We also profiled the composition of the nontemplated addition in the absence of TSO and observed that the 5′-end of RNA template influences the terminal transferase activity of the RT. Furthermore, we found that designing new TSOs that pair with the most common nontemplated additions did little to improve template switching efficiency. Our results provide evidence suggesting that, in contrast to the current understanding of the template switching process, nontemplated addition and template switching are concurrent and competing processes.

Strategies for Converting RNA to Amplifiable cDNA for Single-Cell RNA Sequencing Methods

This review describes the features of molecular biology techniques for single-cell RNA sequencing (scRNA-seq), including methods developed in our laboratory. Existing scRNA-seq methods require the conversion of first-strand cDNA to amplifiable cDNA followed by whole-transcript amplification. There are three primary strategies for this conversion: poly-A tagging, template switching, and RNase H-DNA polymerase I-mediated second-strand cDNA synthesis for in vitro transcription. We discuss the merits and limitations of these strategies and describe our Reverse Transcription with Random Displacement Amplification technology that allows for direct first-strand cDNA amplification from RNA without the need for conversion to an amplifiable cDNA. We believe that this review provides all users of single-cell transcriptome technologies with an understanding of the relationship between the quantitative performance of various methods and their molecular features.

Gene copy number variation in natural populations of Plasmodium falciparum in Eastern Africa

BACKGROUND

Gene copy number variants (CNVs), which consist of deletions and amplifications of single or sets of contiguous genes, contribute to the great diversity in the Plasmodium falciparum genome. In vitro studies in the laboratory have revealed their important role in parasite fitness phenotypes such as red cell invasion, transmissibility and cytoadherence. Studies of natural parasite populations indicate that CNVs are also common in the field and thus may facilitate adaptation of the parasite to its local environment.

RESULTS

In a survey of 183 fresh field isolates from three populations in Eastern Africa with different malaria transmission intensities, we identified 94 CNV loci using microarrays. All CNVs had low population frequencies (minor allele frequency < 5%) but each parasite isolate carried an average of 8 CNVs. Nine CNVs showed high levels of population differentiation (F_ST > 0.3) and nine exhibited significant clines in population frequency across a gradient in transmission intensity. The clearest example of this was a large deletion on chromosome 9 previously reported only in laboratory-adapted isolates. This deletion was present in 33% of isolates from a population with low and highly seasonal malaria transmission, and in < 9% of isolates from populations with higher transmission. Subsets of CNVs were strongly correlated in their population frequencies, implying co-selection.

CONCLUSIONS

These results support the hypothesis that CNVs are the target of selection in natural populations of P. falciparum. Their environment-specific patterns observed here imply an important role for them in conferring adaptability to the parasite thus enabling it to persist in its highly diverse ecological environment.

A20 upregulation during treated HIV disease is associated with intestinal epithelial cell recovery and function

TRIAL REGISTRATION

ClinicalTrials.gov Clinical Trial NCT00594880.

TrpA1 activation in peripheral sensory neurons underlies the ionic basis of pain hypersensitivity in response to vinca alkaloids

Chemotherapy induced peripheral neuropathy (CIPN), a side effect of many anti-cancer drugs including the vinca alkaloids, is characterized by a severe pain syndrome that compromises treatment in many patients. Currently there are no effective treatments for this pain syndrome except for the reduction of anti-cancer drug dose. Existing data supports the model that the pain associated with CIPN is the result of anti-cancer drugs augmenting the function of the peripheral sensory nociceptors but the cellular mechanisms underlying the effects of anti-cancer drugs on sensory neuron function are not well described. Studies from animal models have suggested a number of disease etiologies including mitotoxicity, axonal degeneration, immune signaling, and reduced sensory innervations but these outcomes are the result of prolonged treatment paradigms and do not necessarily represent the early formative events associated with CIPN. Here we show that acute exposure to vinca alkaloids results in an immediate pain syndrome in both flies and mice. Furthermore, we demonstrate that exposure of isolated sensory neurons to vinca alkaloids results in the generation of an inward sodium current capable of depolarizing these neurons to threshold resulting in neuronal firing. These neuronal effects of vinca alkaloids require the transient receptor potential ankyrin-1 (TrpA1) channel, and the hypersensitization to painful stimuli in response to the acute exposure to vinca alkaloids is reduced in TrpA1 mutant flies and mice. These findings demonstrate the direct excitation of sensory neurons by CIPN-causing chemotherapy drugs, and identify TrpA1 as an important target during the pathogenesis of CIPN.

Empirical assessment of analysis workflows for differential expression analysis of human samples using RNA-Seq

BACKGROUND

RNA-Seq has supplanted microarrays as the preferred method of transcriptome-wide identification of differentially expressed genes. However, RNA-Seq analysis is still rapidly evolving, with a large number of tools available for each of the three major processing steps: read alignment, expression modeling, and identification of differentially expressed genes. Although some studies have benchmarked these tools against gold standard gene expression sets, few have evaluated their performance in concert with one another. Additionally, there is a general lack of testing of such tools on real-world, physiologically relevant datasets, which often possess qualities not reflected in tightly controlled reference RNA samples or synthetic datasets.

RESULTS

Here, we evaluate 219 combinatorial implementations of the most commonly used analysis tools for their impact on differential gene expression analysis by RNA-Seq. A test dataset was generated using highly purified human classical and nonclassical monocyte subsets from a clinical cohort, allowing us to evaluate the performance of 495 unique workflows, when accounting for differences in expression units and gene- versus transcript-level estimation. We find that the choice of methodologies leads to wide variation in the number of genes called significant, as well as in performance as gauged by precision and recall, calculated by comparing our RNA-Seq results to those from four previously published microarray and BeadChip analyses of the same cell populations. The method of differential gene expression identification exhibited the strongest impact on performance, with smaller impacts from the choice of read aligner and expression modeler. Many workflows were found to exhibit similar overall performance, but with differences in their calibration, with some biased toward higher precision and others toward higher recall.

CONCLUSIONS

There is significant heterogeneity in the performance of RNA-Seq workflows to identify differentially expressed genes. Among the higher performing workflows, different workflows exhibit a precision/recall tradeoff, and the ultimate choice of workflow should take into consideration how the results will be used in subsequent applications. Our analyses highlight the performance characteristics of these workflows, and the data generated in this study could also serve as a useful resource for future development of software for RNA-Seq analysis.

Integrated analysis of the Plasmodium species transcriptome

The genome sequence available for different Plasmodium species is a valuable resource for understanding malaria parasite biology. However, comparative genomics on its own cannot fully explain all the species-specific differences which suggests that other genomic aspects such as regulation of gene expression play an important role in defining species-specific characteristics. Here, we developed a comprehensive approach to measure transcriptional changes of the evolutionary conserved syntenic orthologs during the intraerythrocytic developmental cycle across six Plasmodium species. We show significant transcriptional constraint at the mid-developmental stage of Plasmodium species while the earliest stages of parasite development display the greatest transcriptional variation associated with critical functional processes. Modeling of the evolutionary relationship based on changes in transcriptional profile reveal a phylogeny pattern of the Plasmodium species that strictly follows its mammalian hosts. In addition, the work shows that transcriptional conserved orthologs represent potential future targets for anti-malaria intervention as they would be expected to carry out key essential functions within the parasites. This work provides an integrated analysis of orthologous transcriptome, which aims to provide insights into the Plasmodium evolution thereby establishing a framework to explore complex pathways and drug discovery in Plasmodium species with broad host range.

•

Comparison of variations in mRNA abundance across six different Plasmodium species.

•

Transcriptional conservation and divergence of Plasmodium syntenic orthologs.

•

Pattern of Plasmodium transcriptome evolution are established.

•

Transcriptionally conserved orthologs represent attractive intervention targets.

Malaria remains a major public health concern despite global efforts in the fight against this disease. The intraerythrocytic stage of the malaria parasites is currently in the spotlight for anti-malarial intervention and vaccine targets. The primary goal of this study is to generate a comprehensive and directly comparable transcriptome dataset across multiple Plasmodium species originating from different hosts. We establish that specific pathways and intraerythrocytic stages are more transcriptionally diverged than others, reflecting transcriptional evolutionary diversity. We further propose a panel of transcriptionally conserved genes as potential drug targets.

Trimming of sequence reads alters RNA-Seq gene expression estimates

BACKGROUND

High-throughput RNA-Sequencing (RNA-Seq) has become the preferred technique for studying gene expression differences between biological samples and for discovering novel isoforms, though the techniques to analyze the resulting data are still immature. One pre-processing step that is widely but heterogeneously applied is trimming, in which low quality bases, identified by the probability that they are called incorrectly, are removed. However, the impact of trimming on subsequent alignment to a genome could influence downstream analyses including gene expression estimation; we hypothesized that this might occur in an inconsistent manner across different genes, resulting in differential bias.

RESULTS

To assess the effects of trimming on gene expression, we generated RNA-Seq data sets from four samples of larval Drosophila melanogaster sensory neurons, and used three trimming algorithms--SolexaQA, Trimmomatic, and ConDeTri-to perform quality-based trimming across a wide range of stringencies. After aligning the reads to the D. melanogaster genome with TopHat2, we used Cuffdiff2 to compare the original, untrimmed gene expression estimates to those following trimming. With the most aggressive trimming parameters, over ten percent of genes had significant changes in their estimated expression levels. This trend was seen with two additional RNA-Seq data sets and with alternative differential expression analysis pipelines. We found that the majority of the expression changes could be mitigated by imposing a minimum length filter following trimming, suggesting that the differential gene expression was primarily being driven by spurious mapping of short reads. Slight differences with the untrimmed data set remained after length filtering, which were associated with genes with low exon numbers and high GC content. Finally, an analysis of paired RNA-seq/microarray data sets suggests that no or modest trimming results in the most biologically accurate gene expression estimates.

CONCLUSIONS

We find that aggressive quality-based trimming has a large impact on the apparent makeup of RNA-Seq-based gene expression estimates, and that short reads can have a particularly strong impact. We conclude that implementation of trimming in RNA-Seq analysis workflows warrants caution, and if used, should be used in conjunction with a minimum read length filter to minimize the introduction of unpredictable changes in expression estimates.

Age-related thermal response: the cellular resilience of juveniles

Understanding species' responses to environmental challenges is key to predicting future biodiversity. However, there is currently little data on how developmental stages affect responses and also whether universal gene biomarkers to environmental stress can be identified both within and between species. Using the Antarctic clam, Laternula elliptica, as a model species, we examined both the tissue-specific and age-related (juvenile versus mature adult) gene expression response to acute non-lethal warming (12 h at 3 °C). In general, there was a relatively muted response to this sub-lethal thermal challenge when the expression profiles of treated animals, of either age, were compared with those of 0 °C controls, with none of the "classical" stress response genes up-regulated. The expression profiles were very variable between the tissues of all animals, irrespective of age with no single transcript emerging as a universal biomarker of thermal stress. However, when the expression profiles of treated animals of the different age groups were directly compared, a very different pattern emerged. The profiles of the younger animals showed significant up-regulation of chaperone and antioxidant transcripts when compared with those of the older animals. Thus, the younger animals showed evidence of a more robust cellular response to warming. These data substantiate previous physiological analyses showing a more resilient juvenile population.

Identification of a gene signature of a pre-transformation process by senescence evasion in normal human epidermal keratinocytes

BACKGROUND

Epidemiological data show that the incidence of carcinomas in humans is highly dependent on age. However, the initial steps of the age-related molecular oncogenic processes by which the switch towards the neoplastic state occurs remain poorly understood, mostly due to the absence of powerful models. In a previous study, we showed that normal human epidermal keratinocytes (NHEKs) spontaneously and systematically escape from senescence to give rise to pre-neoplastic emerging cells.

METHODS

Here, this model was used to analyze the gene expression profile associated with the early steps of age-related cell transformation. We compared the gene expression profiles of growing or senescent NHEKs to post-senescent emerging cells. Data analyses were performed by using the linear modeling features of the limma package, resulting in a two-sided t test or F-test based on moderated statistics. The p-values were adjusted for multiple testing by controlling the false discovery rate according to Benjamini Hochberg method.The common gene set resulting of differential gene expression profiles from these two comparisons revealed a post-senescence neoplastic emergence (PSNE) gene signature of 286 genes.

RESULTS

About half of these genes were already reported as involved in cancer or premalignant skin diseases. However, bioinformatics analyses did not highlight inside this signature canonical cancer pathways but metabolic pathways, including in first line the metabolism of xenobiotics by cytochrome P450. In order to validate the relevance of this signature as a signature of pretransformation by senescence evasion, we invalidated two components of the metabolism of xenobiotics by cytochrome P450, AKR1C2 and AKR1C3. When performed at the beginning of the senescence plateau, this invalidation did not alter the senescent state itself but significantly decreased the frequency of PSNE. Conversely, overexpression of AKR1C2 but not AKR1C3 increased the frequency of PSNE.

CONCLUSIONS

To our knowledge, this study is the first to identify reprogrammation of metabolic pathways in normal keratinocytes as a potential determinant of the switch from senescence to pre-transformation.

Differences in the gene expression profiles of haemocytes from schistosome-susceptible and -resistant biomphalaria glabrata exposed to Schistosoma mansoni excretory-secretory products

During its life cycle, the helminth parasite Schistosoma mansoni uses the freshwater snail Biomphalaria glabrata as an intermediate host to reproduce asexually generating cercariae for infection of the human definitive host. Following invasion of the snail, the parasite develops from a miracidium to a mother sporocyst and releases excretory-secretory products (ESPs) that likely influence the outcome of host infection. To better understand molecular interactions between these ESPs and the host snail defence system, we determined gene expression profiles of haemocytes from S. mansoni-resistant or -susceptible strains of B. glabrata exposed in vitro to S. mansoni ESPs (20 μg/ml) for 1 h, using a 5K B. glabrata cDNA microarray. Ninety-eight genes were found differentially expressed between haemocytes from the two snail strains, 57 resistant specific and 41 susceptible specific, 60 of which had no known homologue in GenBank. Known differentially expressed resistant-snail genes included the nuclear factor kappa B subunit Relish, elongation factor 1α, 40S ribosomal protein S9, and matrilin; known susceptible-snail specific genes included cathepsins D and L, and theromacin. Comparative analysis with other gene expression studies revealed 38 of the 98 identified genes to be uniquely differentially expressed in haemocytes in the presence of ESPs, thus identifying for the first time schistosome ESPs as important molecules that influence global snail host-defence cell gene expression profiles. Such immunomodulation may benefit the schistosome, enabling its survival and successful development in the snail host.

Fidelity and representativeness of two isothermal multiple displacement amplification systems to preamplify limiting amounts of total RNA

In this study we investigated the fidelity and representativeness of two novel multiple displacement amplification (MDA) protocols leading to whole transcriptome amplification (WTA). WTA is used to amplify a limiting amount of experimental RNA, allowing its use in downstream applications. Using Phi29 and Bst DNA polymerase-based MDA, henceforth referred to as WTA-Phi and WTA-Bst, respectively, we successfully amplified very low amounts of linearly concatenated cDNA originating from 10 to 100 ng of starting RNA. The average yield obtained from 10 ng was 3.5 and 4.7 μg for WTA-Phi and WTA-Bst, respectively, while 100 ng of starting RNA yielded 7.0 and 12.4 μg for WTA-Phi and WTA-Bst, respectively. Representational distortion of the templates, analyzed via conventional PCR, showed robust amplification of 11 different transcripts when either WTA-Phi or WTA-Bst synthesized templates were used, while some transcripts were not detected from unamplified templates. Loci representation, a measure of amplification consistency, was evaluated using TaqMan RT-qPCR amplification of five different transcripts, yielding values ranging from 96.4 to 189.3 %, comparable to those obtained using genomic target-based MDA systems. The two MDA protocols described in this study efficiently lead to representative WTA, using as little as 10 ng of starting RNA.

Hypoxia impacts large adults first: consequences in a warming world

Future oceans are predicted to contain less oxygen than at present. This is because oxygen is less soluble in warmer water and predicted stratification will reduce mixing. Hypoxia in marine environments is thus likely to become more widespread in marine environments and understanding species-responses is important to predicting future impacts on biodiversity. This study used a tractable model, the Antarctic clam, Laternula elliptica, which can live for 36 years, and has a well-characterized ecology and physiology to understand responses to hypoxia and how the effect varied with age. Younger animals had a higher condition index, higher adenylate energy charge and transcriptional profiling indicated that they were physically active in their response to hypoxia, whereas older animals were more sedentary, with higher levels of oxidative damage and apoptosis in the gills. These effects could be attributed, in part, to age-related tissue scaling; older animals had proportionally less contractile muscle mass and smaller gills and foot compared with younger animals, with consequential effects on the whole-animal physiological response. The data here emphasize the importance of including age effects, as large mature individuals appear to be less able to resist hypoxic conditions and this is the size range that is the major contributor to future generations. Thus, the increased prevalence of hypoxia in future oceans may have marked effects on benthic organisms' abilities to persist and this is especially so for long-lived species when predicting responses to environmental perturbation.

Early differential gene expression in haemocytes from resistant and susceptible Biomphalaria glabrata strains in response to Schistosoma mansoni

The outcome of infection in the host snail Biomphalaria glabrata with the digenean parasite Schistosoma mansoni is determined by the initial molecular interplay occurring between them. The mechanisms by which schistosomes evade snail immune recognition to ensure survival are not fully understood, but one possibility is that the snail internal defence system is manipulated by the schistosome enabling the parasite to establish infection. This study provides novel insights into the nature of schistosome resistance and susceptibility in B. glabrata at the transcriptomic level by simultaneously comparing gene expression in haemocytes from parasite-exposed and control groups of both schistosome-resistant and schistosome-susceptible strains, 2 h post exposure to S. mansoni miracidia, using an novel 5K cDNA microarray. Differences in gene expression, including those for immune/stress response, signal transduction and matrix/adhesion genes were identified between the two snail strains and tests for asymmetric distributions of gene function also identified immune-related gene expression in resistant snails, but not in susceptible. Gene set enrichment analysis revealed that genes involved in mitochondrial electron transport, ubiquinone biosynthesis and electron carrier activity were consistently up-regulated in resistant snails but down-regulated in susceptible. This supports the hypothesis that schistosome-resistant snails recognize schistosomes and mount an appropriate defence response, while in schistosome-susceptible snails the parasite suppresses this defence response, early in infection.

Options available for labelling nucleic acid samples in DNA microarray-based detection methods

DNA microarrays are considered by many researchers to be the platform of choice for the high-throughput analysis of nucleic acids. Since the past two decades, they have been used constantly as powerful tools in differential gene expression, SNP genotyping, DNA sequencing, gene discovery, disease diagnostic and pathways reconstruction. Several methods have been developed to enable samples of limited amounts of RNA to be quantified. Here we evaluate classical and up-to-date assays made available for labelling those samples. This review also sheds light on the recently developed strategies that ensure high sensitivity such as sample and signal amplification, quantum dot, surface plasmom resonance, nanoparticles and cationinc polythiophenes.

Highly parallel genome-wide expression analysis of single mammalian cells

BACKGROUND

We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs.

METHODOLOGY/PRINCIPAL FINDINGS

The method uses tagged priming and template-switching, resulting in the incorporation of universal PCR priming sites at both ends of the synthesized cDNA for global PCR amplification. Coupled with a whole-genome gene expression microarray platform, we routinely obtain expression correlation values of R(2)~0.76-0.80 between individual cells and R(2)~0.69 between 50 pg total RNA replicates. Expression profiles generated from single cells or 50 pg total RNA correlate well with that generated with higher input (1 ng total RNA) (R(2)~0.80). Also, the assay is sufficiently sensitive to detect, in a single cell, approximately 63% of the number of genes detected with 1 ng input, with approximately 97% of the genes detected in the single-cell input also detected in the higher input.

CONCLUSIONS/SIGNIFICANCE

In summary, our method facilitates whole-genome gene expression profiling in contexts where starting material is extremely limiting, particularly in areas such as the study of progenitor cells in early development and tumor stem cell biology.

Gene expression profiling of cancer stem cells in the Du145 prostate cancer cell line

Cancer stem cells (CSCs) that exhibit tumor-initiating properties have been identified in primary prostate cancer and prostate cancer cell lines. CSCs are the root of tumor metastasis and recurrence even in the cases where they constitute a minority of the tumor mass. In this study, putative CSCs were isolated from the Du145 cell line and an Affymetrix microarray was used to investigate their gene expression profile. The results were validated by real-time polymerase chain reaction (PCR). The results of the microarray indicated that 138 genes were upregulated and 93 genes were downregulated in the CSCs. Certain genes that may be significant in the regulation of CSCs from the Du145 cell line were identified. These results may aid the studying of the mechanisms through which CSCs acquire their distinctive properties.

Oligonucleotide microarray expression profiling of contrasting invasive phenotypes in colorectal cancer

This chapter refers to the application of laser-capture microdissection with oligonucleotide microarray analysis. The protocol described has been successfully used to identify differential transcript expression between contrasting colorectal cancer invasive phenotypes. Tissue processing, RNA extraction, quality control, amplification, fluorescent labelling, purification, hybridisation, and elements of data analysis are covered.

Differential gene expression during the moult cycle of Antarctic krill (Euphausia superba)

Background

All crustaceans periodically moult to renew their exoskeleton. In krill this involves partial digestion and resorption of the old exoskeleton and synthesis of new cuticle. Molecular events that underlie the moult cycle are poorly understood in calcifying crustaceans and even less so in non-calcifying organisms such as krill. To address this we constructed an Antarctic krill cDNA microarray in order to generate gene expression profiles across the moult cycle and identify possible activation pathways.

Results

A total of 26 different cuticle genes were identified that showed differential gene expression across the moult cycle. Almost all cuticle genes were up regulated during premoult and down regulated during late intermoult. There were a number of transcripts with significant sequence homology to genes potentially involved in the synthesis, breakdown and resorption of chitin. During early premoult glutamine synthetase, a gene involved in generating an amino acid used in the synthesis of glucosamine, a constituent of chitin, was up regulated more than twofold. Mannosyltransferase 1, a member of the glycosyltransferase family of enzymes that includes chitin synthase was also up regulated during early premoult. Transcripts homologous to a β-N-acetylglucosaminidase (β-NAGase) precursor were expressed at a higher level during late intermoult (prior to apolysis) than during premoult. This observation coincided with the up regulation during late intermoult, of a coatomer subunit epsilon involved in the production of vesicles that maybe used to transport the β-NAGase precursors into the exuvial cleft. Trypsin, known to activate the β-NAGase precursor, was up regulated more than fourfold during premoult. The up regulation of a predicted oligopeptide transporter during premoult may allow the transport of chitin breakdown products across the newly synthesised epi- and exocuticle layers.

Conclusion

We have identified many genes differentially expressed across the moult cycle of krill that correspond with known phenotypic structural changes. This study has provided a better understanding of the processes involved in krill moulting and how they may be controlled at the gene expression level.

Incorporation of non-natural nucleotides into template-switching oligonucleotides reduces background and improves cDNA synthesis from very small RNA samples

Background

The template switching PCR (TS-PCR) method of cDNA synthesis represents one of the most straightforward approaches to generating full length cDNA for sequencing efforts. However, when applied to very small RNA samples, such as those obtained from tens or hundreds of cells, this approach leads to high background and low cDNA yield due to concatamerization of the TS oligo.

Results

In this study, we describe the application of nucleotide isomers that form non-standard base pairs in the template switching oligo to prevent background cDNA synthesis. When such bases are added to the 5' end of the template switching (TS) oligo, they inhibit MMLV-RT from extending the cDNA beyond the TS oligo, thus increasing cDNA yield by reducing formation of concatamers of the TS oligo that are the source of significant background.

Conclusions

Our results demonstrate that this novel approach for cDNA synthesis has valuable utility for application of ultra-high throughput technologies, such as whole transcriptome sequencing using 454 technology, to very small biological samples comprised of tens of cells as might be obtained via approaches like laser microdissection.

Insights into shell deposition in the Antarctic bivalve Laternula elliptica: gene discovery in the mantle transcriptome using 454 pyrosequencing

Background

The Antarctic clam, Laternula elliptica, is an infaunal stenothermal bivalve mollusc with a circumpolar distribution. It plays a significant role in bentho-pelagic coupling and hence has been proposed as a sentinel species for climate change monitoring. Previous studies have shown that this mollusc displays a high level of plasticity with regard to shell deposition and damage repair against a background of genetic homogeneity. The Southern Ocean has amongst the lowest present-day CaCO₃ saturation rate of any ocean region, and is predicted to be among the first to become undersaturated under current ocean acidification scenarios. Hence, this species presents as an ideal candidate for studies into the processes of calcium regulation and shell deposition in our changing ocean environments.

Results

454 sequencing of L. elliptica mantle tissue generated 18,290 contigs with an average size of 535 bp (ranging between 142 bp-5.591 kb). BLAST sequence similarity searching assigned putative function to 17% of the data set, with a significant proportion of these transcripts being involved in binding and potentially of a secretory nature, as defined by GO molecular function and biological process classifications. These results indicated that the mantle is a transcriptionally active tissue which is actively proliferating. All transcripts were screened against an in-house database of genes shown to be involved in extracellular matrix formation and calcium homeostasis in metazoans. Putative identifications were made for a number of classical shell deposition genes, such as tyrosinase, carbonic anhydrase and metalloprotease 1, along with novel members of the family 2 G-Protein Coupled Receptors (GPCRs). A membrane transport protein (SEC61) was also characterised and this demonstrated the utility of the clam sequence data as a resource for examining cold adapted amino acid substitutions. The sequence data contained 46,235 microsatellites and 13,084 Single Nucleotide Polymorphisms(SNPs/INDELS), providing a resource for population and also gene function studies.

Conclusions

This is the first 454 data from an Antarctic marine invertebrate. Sequencing of mantle tissue from this non-model species has considerably increased resources for the investigation of the processes of shell deposition and repair in molluscs in a changing environment. A number of promising candidate genes were identified for functional analyses, which will be the subject of further investigation in this species and also used in model-hopping experiments in more tractable and economically important model aquaculture species, such as Crassostrea gigas and Mytilus edulis.

Evaluation of quantitative rt-PCR using nonamplified and amplified RNA

Quantitative reverse-transcription polymerase chain reaction (qRT-PCR) is a standard assay in molecular medicine for gene expression analysis. Samples from incisional/needle biopsies, laser-microdissected tumor cells and other biologic sources, normally available in clinical cancer studies, generate very small amounts of RNA that are restrictive for expression analysis. As a consequence, an RNA amplification procedure is required to assess the gene expression levels of such sample types. The reproducibility and accuracy of relative gene expression data produced by sensitive methodology as qRT-PCR when cDNA converted from amplified (A) RNA is used as template has not yet been properly addressed. In this study, to properly evaluate this issue, we performed 1 round of linear RNA amplification in 2 breast cell lines (C5.2 and HB4a) and assessed the relative expression of 34 genes using cDNA converted from both nonamplified (NA) and A RNA. Relative gene expression was obtained from beta actin or glyceraldehyde 3-phosphate dehydrogenase normalized data using different dilutions of cDNA, wherein the variability and fold-change differences in the expression of the 2 methods were compared. Our data showed that 1 round of linear RNA amplification, even with suboptimal-quality RNA, is appropriate to generate reproducible and high-fidelity qRT-PCR relative expression data that have similar confidence levels as those from NA samples. The use of cDNA that is converted from both A and NA RNA in a single qRT-PCR experiment clearly creates bias in relative gene expression data.

Loading-related regulation of gene expression in bone in the contexts of estrogen deficiency, lack of estrogen receptor alpha and disuse

Loading-related changes in gene expression in resident cells in the tibia of female mice in the contexts of normality (WT), estrogen deficiency (WT-OVX), absence of estrogen receptor alpha (ERalpha(-/-)) and disuse due to sciatic neurectomy (WT-SN) were established by microarray. Total RNA was extracted from loaded and contra-lateral non-loaded tibiae at selected time points after a single, short period of dynamic loading sufficient to engender an osteogenic response. There were marked changes in the expression of many genes according to context as well as in response to loading within those contexts. In WT mice at 3, 8, 12 and 24 h after loading the expression of 642, 341, 171 and 24 genes, respectively, were differentially regulated compared with contra-lateral bones which were not loaded. Only a few of the genes differentially regulated by loading in the tibiae of WT mice have recognized roles in bone metabolism or have been linked previously to osteogenesis (Opn, Sost, Esr1, Tgfb1, Lrp1, Ostn, Timp, Mmp, Ctgf, Postn and Irs1, BMP and DLX5). The canonical pathways showing the greatest loading-related regulation were those involving pyruvate metabolism, mitochondrial dysfunction, calcium-induced apoptosis, glycolysis/gluconeogenesis, aryl hydrocarbon receptor and oxidative phosphorylation. In the tibiae from WT-OVX, ERalpha(-/-) and WT-SN mice, 440, 439 and 987 genes respectively were differentially regulated by context alone compared to WT. The early response to loading in tibiae of WT-OVX mice involved differential regulation compared to their contra-lateral non-loaded pair of fewer genes than in WT, more down-regulation than up-regulation and a later response. This was shared by WT-SN. In tibiae of ERalpha(-/-) mice, the number of genes differentially regulated by loading was markedly reduced at all time points. These data indicate that in resident bone cells, both basal and loading-related gene expression is substantially modified by context. Many of the genes differentially regulated by the earliest loading-related response were primarily involved in energy metabolism and were not specific to bone.

Organ-specific defence strategies of pepper (Capsicum annuum L.) during early phase of water deficit

Drought is one of the major factors that limits crop production and reduces yield. To understand the early response of plants under nearly natural conditions, pepper plants (Capsicum annuum L.) were grown in a greenhouse and stressed by withholding water for 1 week. Plants adapted to the decreasing water content of the soil by adjustment of their osmotic potential in root tissue. As a consequence of drought, strong accumulation of raffinose, glucose, galactinol and proline was detected in the roots. In contrast, in leaves the levels of fructose, sucrose and also galactinol increased. Due to the water deficit cadaverine, putrescine, spermidine and spermine accumulated in leaves, whereas the concentration of polyamines was reduced in roots. To study the molecular basis of these responses, a combined approach of suppression subtractive hybridisation and microarray technique was performed on the same material. A total of 109 unique ESTs were detected as responsive to drought, while additional 286 ESTs were selected from the bulk of rare transcripts on the array. The metabolic profiles of stressed pepper plants are discussed with respect to the transcriptomic changes detected, while attention is given to the differences between defence strategies of roots and leaves.

Gene expression associated with changes in cold tolerance levels of the Antarctic springtail, Cryptopygus antarcticus

The ability of the Antarctic microarthropod Cryptopygus antarcticus (Collembola, Isotomidae) to survive low temperatures has been well studied at the physiological level, with recent investigations indicating the importance of the moulting process in conferring this ability. This study investigated gene expression in groups of C. antarcticus that have distinct differences in their ability to survive low temperatures. A microarray containing c. 5400 C. antarcticus expressed sequence tags was used to investigate gene expression differences between groups of animals with different supercooling points (SCP), and to low temperatures close to their SCP. By demonstrating the involvement of moult-related genes in the differential survival of two groups of C. antarcticus with distinct SCP profiles, the results of this investigation add support to the suggestion that moulting plays a role in conferring cold tolerance in C. antarcticus.

Time series analysis of the transcriptional responses of Biomphalaria glabrata throughout the course of intramolluscan development of Schistosoma mansoni and Echinostoma paraensei

Successful colonization of a compatible snail host by a digenetic trematode miracidium initiates a complex, proliferative development program requiring weeks to reach culmination in the form of production of cercariae which, once started, may persist for the remainder of the life span of the infected snail. How are such proliferative and invasive parasites able to circumvent host defenses and establish chronic infections? Using a microarray designed to monitor the internal defense and stress-related responses of the freshwater snail Biomphalaria glabrata, we have undertaken a time course study to monitor snail responses following exposure to two different trematode species to which the snail is susceptible: the medically important Schistosoma mansoni, exemplifying sporocyst production in its larval development, or Echinostoma paraensei, representing an emphasis on rediae production in its larval development. We sampled eight time points (0.5, 1, 2, 4, 8, 16 and 32 days p.i.) that cover the period required for cercariae to be produced. Following exposure to S. mansoni, there was a preponderance of up-regulated over down-regulated array features through 2 days p.i. but by 4 days p.i. and thereafter, this pattern was strongly reversed. For E. paraensei, there was a preponderance of down-regulated array features over up-regulated features at even 0.5 days p.i., a pattern that persists throughout the course of infection except for 1 day p.i., when up-regulated array features slightly outnumbered down-regulated features. Examination of particular array features revealed several that were up-regulated by both parasites early in the course of infection and one, fibrinogen related protein 4 (FREP 4), that remained significantly elevated throughout the course of infection with either parasite, effectively serving as a marker of infection. Many defense-related transcripts were persistently down-regulated, including several fibrinogen-containing lectins and homologs of molecules best known from vertebrate phagocytic cells. Our results are consistent with earlier studies suggesting that both parasites are able to interfere with host defense responses, including a tendency for E. paraensei to do so more rapidly and strongly than S. mansoni. They further suggest mechanisms for how trematodes are able to establish the chronic infections necessary for their continued success.

In vitro manipulation of gene expression in larval Schistosoma: a model for postgenomic approaches in Trematoda

With rapid developments in DNA and protein sequencing technologies, combined with powerful bioinformatics tools, a continued acceleration of gene identification in parasitic helminths is predicted, potentially leading to discovery of new drug and vaccine targets, enhanced diagnostics and insights into the complex biology underlying host-parasite interactions. For the schistosome blood flukes, with the recent completion of genome sequencing and comprehensive transcriptomic datasets, there has accumulated massive amounts of gene sequence data, for which, in the vast majority of cases, little is known about actual functions within the intact organism. In this review we attempt to bring together traditional in vitro cultivation approaches and recent emergent technologies of molecular genomics, transcriptomics and genetic manipulation to illustrate the considerable progress made in our understanding of trematode gene expression and function during development of the intramolluscan larval stages. Using several prominent trematode families (Schistosomatidae, Fasciolidae, Echinostomatidae), we have focused on the current status of in vitro larval isolation/cultivation as a source of valuable raw material supporting gene discovery efforts in model digeneans that include whole genome sequencing, transcript and protein expression profiling during larval development, and progress made in the in vitro manipulation of genes and their expression in larval trematodes using transgenic and RNA interference (RNAi) approaches.

Differential transcriptomic responses of Biomphalaria glabrata (Gastropoda, Mollusca) to bacteria and metazoan parasites, Schistosoma mansoni and Echinostoma paraensei (Digenea, Platyhelminthes)

A 70-mer-oligonucleotide-based microarray (1152 features) that emphasizes stress and immune responses factors was constructed to study transcriptomic responses of the snail Biomphalaria glabrata to different immune challenges. In addition to sequences with relevant putative ID and Gene Ontology (GO) annotation, the array features non-immune factors and unknown B. glabrata ESTs for functional gene discovery. The transcription profiles of B. glabrata (3 biological replicates, each a pool of 5 snails) were recorded at 12h post-wounding, exposure to Gram negative or Gram positive bacteria (Escherichia coli and Micrococcus luteus, respectively), or infection with compatible trematode parasites (Schistosoma mansoni or Echinostoma paraensei, 20 miracidia/snail), relative to controls, using universal reference RNA. The data were subjected to Significance Analysis for Microarrays (SAM), with a false positive rate (FPR) <or=10%. Wounding yielded a modest differential expression profile (27 up/21 down) with affected features mostly dissimilar from other treatments. Partially overlapping, yet distinct expression profiles were recorded from snails challenged with E. coli (83 up/20 down) or M. luteus (120 up/42 down), mostly showing up-regulation of defense and stress-related features. Significantly altered expression of selected immune features indicates that B. glabrata detects and responds differently to compatible trematodes. Echinostoma paraensei infection was associated mostly with down-regulation of many (immune-) transcripts (42 up/68 down), whereas S. mansoni exposure yielded a preponderance of up-regulated features (140 up/23 down), with only few known immune genes affected. These observations may reflect the divergent strategies developed by trematodes during their evolution as specialized pathogens of snails to negate host defense responses. Clearly, the immune defenses of B. glabrata distinguish and respond differently to various immune challenges.

Global transcriptional response of pig brain and lung to natural infection by Pseudorabies virus

Background

Pseudorabies virus (PRV) is an alphaherpesviruses whose native host is pig. PRV infection mainly causes signs of central nervous system disorder in young pigs, and respiratory system diseases in the adult.

Results

In this report, we have analyzed native host (piglets) gene expression changes in response to acute pseudorabies virus infection of the brain and lung using a printed human oligonucleotide gene set from Illumina. A total of 210 and 1130 out of 23,000 transcript probes displayed differential expression respectively in the brain and lung in piglets after PRV infection (p-value < 0.01), with most genes displaying up-regulation. Biological process and pathways analysis showed that most of the up-regulated genes are involved in cell differentiation, neurodegenerative disorders, the nervous system and immune responses in the infected brain whereas apoptosis, cell cycle control, and the mTOR signaling pathway genes were prevalent in the infected lung. Additionally, a number of differentially expressed genes were found to map in or close to quantitative trait loci for resistance/susceptibility to pseudorabies virus in piglets.

Conclusion

This is the first comprehensive analysis of the global transcriptional response of the native host to acute alphaherpesvirus infection. The differentially regulated genes reported here are likely to be of interest for the further study and understanding of host viral gene interactions.

Anti-schistosomal intervention targets identified by lifecycle transcriptomic analyses

Background

Novel methods to identify anthelmintic drug and vaccine targets are urgently needed, especially for those parasite species currently being controlled by singular, often limited strategies. A clearer understanding of the transcriptional components underpinning helminth development will enable identification of exploitable molecules essential for successful parasite/host interactions. Towards this end, we present a combinatorial, bioinformatics-led approach, employing both statistical and network analyses of transcriptomic data, for identifying new immunoprophylactic and therapeutic lead targets to combat schistosomiasis.

Methodology/Principal Findings

Utilisation of a Schistosoma mansoni oligonucleotide DNA microarray consisting of 37,632 elements enabled gene expression profiling from 15 distinct parasite lifecycle stages, spanning three unique ecological niches. Statistical approaches of data analysis revealed differential expression of 973 gene products that minimally describe the three major characteristics of schistosome development: asexual processes within intermediate snail hosts, sexual maturation within definitive vertebrate hosts and sexual dimorphism amongst adult male and female worms. Furthermore, we identified a group of 338 constitutively expressed schistosome gene products (including 41 transcripts sharing no sequence similarity outside the Platyhelminthes), which are likely to be essential for schistosome lifecycle progression. While highly informative, statistics-led bioinformatics mining of the transcriptional dataset has limitations, including the inability to identify higher order relationships between differentially expressed transcripts and lifecycle stages. Network analysis, coupled to Gene Ontology enrichment investigations, facilitated a re-examination of the dataset and identified 387 clusters (containing 12,132 gene products) displaying novel examples of developmentally regulated classes (including 294 schistosomula and/or adult transcripts with no known sequence similarity outside the Platyhelminthes), which were undetectable by the statistical comparisons.

Conclusions/Significance

Collectively, statistical and network-based exploratory analyses of transcriptomic datasets have led to a thorough characterisation of schistosome development. Information obtained from these experiments highlighted key transcriptional programs associated with lifecycle progression and identified numerous anti-schistosomal candidate molecules including G-protein coupled receptors, tetraspanins, Dyp-type peroxidases, fucosyltransferases, leishmanolysins and the netrin/netrin receptor complex.

Despite the implementation of focused and well-funded chemotherapeutic control initiatives over the last decade, schistosomiasis remains a significant cause of morbidity and mortality within countries of the developing world. There is, therefore, an urgent need for the rapid translation of genomic information into viable vaccines or new drug classes capable of eradicating the parasitic schistosome worms responsible for this neglected tropical disease. In our effort to identify potential targets for novel chemotherapeutic and immunoprophylactic interventions, we detail a combined bioinformatics approach, comprising exploratory statistical and network analyses, to thoroughly describe the transcriptional progression of Schistosoma mansoni across three environmental niches. Our results indicate that, although schistosomes are masters at host deception and survival, there are numerous exploitable candidate molecules displaying either differential or constitutive expression throughout the parasite's lifecycle. Importantly, some of these transcripts represent gene families not commonly found outside—or expanded within—the phylum Platyhelminthes, and thus represent priority targets. Many of the candidates identified herein will be subjected to ongoing and future hypothesis-led functional investigations. The completion of such specific examinations ultimately will contribute to the successful development of novel control strategies useful in the alleviation of schistosome-induced immunopathologies, morbidities and mortalities.

Comparative transcriptional and genomic analysis of Plasmodium falciparum field isolates

Mechanisms for differential regulation of gene expression may underlie much of the phenotypic variation and adaptability of malaria parasites. Here we describe transcriptional variation among culture-adapted field isolates of Plasmodium falciparum, the species responsible for most malarial disease. It was found that genes coding for parasite protein export into the red cell cytosol and onto its surface, and genes coding for sexual stage proteins involved in parasite transmission are up-regulated in field isolates compared with long-term laboratory isolates. Much of this variability was associated with the loss of small or large chromosomal segments, or other forms of gene copy number variation that are prevalent in the P. falciparum genome (copy number variants, CNVs). Expression levels of genes inside these segments were correlated to that of genes outside and adjacent to the segment boundaries, and this association declined with distance from the CNV boundary. This observation could not be explained by copy number variation in these adjacent genes. This suggests a local-acting regulatory role for CNVs in transcription of neighboring genes and helps explain the chromosomal clustering that we observed here. Transcriptional co-regulation of physical clusters of adaptive genes may provide a way for the parasite to readily adapt to its highly heterogeneous and strongly selective environment.

T7-based linear amplification of low concentration mRNA samples using beads and microfluidics for global gene expression measurements

We have demonstrated in vitro transcription (IVT) of cDNA sequences from purified Jurkat T-cell mRNA immobilized on microfluidic packed beds down to single-cell quantities. The microfluidically amplified antisense-RNA (aRNA) was nearly identical in length and quantity compared with benchtop reactions using the same starting sample quantities. Microarrays were used to characterize the number and population of genes in each sample, allowing comparison of the microfluidic and benchtop processes. For both benchtop and microfluidic assays, we measured the expression of approximately 4000 to 9000 genes for sample amounts ranging from 20 pg to 10 ng (2 to 1000 cell equivalents), corresponding to 41 to 93% of the absolute number of genes detected from a 100 ng total RNA control sample. Concordance of genes detected between methods (benchtop vs. microfluidic) and repeats (microfluidic vs. microfluidic) typically exceeded 90%. Validation of microarray by Real-time PCR of a panel of five genes suggests transcription of genes present is approximately six times more efficient with the microfluidic IVT compared with benchtop processing. Microfluidic IVT introduces no bias to the gene expression profile of the sample and provides more efficient transcription of mRNA sequences present at the single-cell level.

The Maestro (Mro) gene is dispensable for normal sexual development and fertility in mice

The mammalian gonad arises as a bipotential primordium from which a testis or ovary develops depending on the chromosomal sex of the individual. We have previously used DNA microarrays to screen for novel genes controlling the developmental fate of the indifferent embryonic mouse gonad. Maestro (Mro), which encodes a HEAT-repeat protein, was originally identified as a gene exhibiting sexually dimorphic expression during mouse gonad development. Wholemount in situ hybridisation analysis revealed Mro to be expressed in the embryonic male gonad from approximately 11.5 days post coitum, prior to overt sexual differentiation. No significant expression was detected in female gonads at the same developmental stage. In order to address its physiological function, we have generated mice lacking Maestro using gene targeting. Male and female mice homozygous for a Mro null allele are viable and fertile. We examined gonad development in homozygous male embryos in detail and observed no differences when compared to wild-type controls. Immunohistochemical analysis of homozygous mutant testes of adult mice revealed no overt abnormalities. Expression profiling using DNA microarrays also indicated no significant differences between homozygote embryonic male gonads and controls. We conclude that Maestro is dispensable for normal male sexual development and fertility in laboratory mice; however, the Mro locus itself does have utility as a site for insertion of transgenes for future studies in the fields of sexual development and Sertoli cell function.

Biomphalaria glabrata transcriptome: cDNA microarray profiling identifies resistant- and susceptible-specific gene expression in haemocytes from snail strains exposed to Schistosoma mansoni

Background

Biomphalaria glabrata is an intermediate snail host for Schistosoma mansoni, one of the important schistosomes infecting man. B. glabrata/S. mansoni provides a useful model system for investigating the intimate interactions between host and parasite. Examining differential gene expression between S. mansoni-exposed schistosome-resistant and susceptible snail lines will identify genes and pathways that may be involved in snail defences.

Results

We have developed a 2053 element cDNA microarray for B. glabrata containing clones from ORESTES (Open Reading frame ESTs) libraries, suppression subtractive hybridization (SSH) libraries and clones identified in previous expression studies. Snail haemocyte RNA, extracted from parasite-challenged resistant and susceptible snails, 2 to 24 h post-exposure to S. mansoni, was hybridized to the custom made cDNA microarray and 98 differentially expressed genes or gene clusters were identified, 94 resistant-associated and 4 susceptible-associated. Quantitative PCR analysis verified the cDNA microarray results for representative transcripts. Differentially expressed genes were annotated and clustered using gene ontology (GO) terminology and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis. 61% of the identified differentially expressed genes have no known function including the 4 susceptible strain-specific transcripts. Resistant strain-specific expression of genes implicated in innate immunity of invertebrates was identified, including hydrolytic enzymes such as cathepsin L, a cysteine proteinase involved in lysis of phagocytosed particles; metabolic enzymes such as ornithine decarboxylase, the rate-limiting enzyme in the production of polyamines, important in inflammation and infection processes, as well as scavenging damaging free radicals produced during production of reactive oxygen species; stress response genes such as HSP70; proteins involved in signalling, such as importin 7 and copine 1, cytoplasmic intermediate filament (IF) protein and transcription enzymes such as elongation factor 1α and EF-2.

Conclusion

Production of the first cDNA microarray for profiling gene expression in B. glabrata provides a foundation for expanding our understanding of pathways and genes involved in the snail internal defence system (IDS). We demonstrate resistant strain-specific expression of genes potentially associated with the snail IDS, ranging from signalling and inflammation responses through to lysis of proteinacous products (encapsulated sporocysts or phagocytosed parasite components) and processing/degradation of these targeted products by ubiquitination.

Use of genomic DNA as an indirect reference for identifying gender-associated transcripts in morphologically identical, but chromosomally distinct, Schistosoma mansoni cercariae

BACKGROUND

The use of DNA microarray technology to study global Schistosoma gene expression has led to the rapid identification of novel biological processes, pathways or associations. Implementation of standardized DNA microarray protocols across laboratories would assist maximal interpretation of generated datasets and extend productive application of this technology.

METHODOLOGY/PRINCIPAL FINDINGS

Utilizing a new Schistosoma mansoni oligonucleotide DNA microarray composed of 37,632 elements, we show that schistosome genomic DNA (gDNA) hybridizes with less variation compared to complex mixed pools of S. mansoni cDNA material (R = 0.993 for gDNA compared to R = 0.956 for cDNA during 'self versus self' hybridizations). Furthermore, these effects are species-specific, with S. japonicum or Mus musculus gDNA failing to bind significantly to S. mansoni oligonucleotide DNA microarrays (e.g R = 0.350 when S. mansoni gDNA is co-hybridized with S. japonicum gDNA). Increased median fluorescent intensities (209.9) were also observed for DNA microarray elements hybridized with S. mansoni gDNA compared to complex mixed pools of S. mansoni cDNA (112.2). Exploiting these valuable characteristics, S. mansoni gDNA was used in two-channel DNA microarray hybridization experiments as a common reference for indirect identification of gender-associated transcripts in cercariae, a schistosome life-stage in which there is no overt sexual dimorphism. This led to the identification of 2,648 gender-associated transcripts. When compared to the 780 gender-associated transcripts identified by hybridization experiments utilizing a two-channel direct method (co-hybridization of male and female cercariae cDNA), indirect methods using gDNA were far superior in identifying greater quantities of differentially expressed transcripts. Interestingly, both methods identified a concordant subset of 188 male-associated and 156 female-associated cercarial transcripts, respectively. Gene ontology classification of these differentially expressed transcripts revealed a greater diversity of categories in male cercariae. Quantitative real-time PCR analysis confirmed the DNA microarray results and supported the reliability of this platform for identifying gender-associated transcripts.

CONCLUSIONS/SIGNIFICANCE

Schistosome gDNA displays characteristics highly suitable for the comparison of two-channel DNA microarray results obtained from experiments conducted independently across laboratories. The schistosome transcripts identified here demonstrate, for the first time, that gender-associated patterns of expression are already well established in the morphologically identical, but chromosomally distinct, cercariae stage.

Development and application of the active surveillance of pathogens microarray to monitor bacterial gene flux

Background

Human and animal health is constantly under threat by emerging pathogens that have recently acquired genetic determinants that enhance their survival, transmissibility and virulence. We describe the construction and development of an Active Surveillance of Pathogens (ASP) oligonucleotide microarray, designed to 'actively survey' the genome of a given bacterial pathogen for virulence-associated genes.

Results

The microarray consists of 4958 reporters from 151 bacterial species and include genes for the identification of individual bacterial species as well as mobile genetic elements (transposons, plasmid and phage), virulence genes and antibiotic resistance genes. The ASP microarray was validated with nineteen bacterial pathogens species, including Francisella tularensis, Clostridium difficile, Staphylococcus aureus, Enterococcus faecium and Stenotrophomonas maltophilia. The ASP microarray identified these bacteria, and provided information on potential antibiotic resistance (eg sufamethoxazole resistance and sulfonamide resistance) and virulence determinants including genes likely to be acquired by horizontal gene transfer (e.g. an alpha-haemolysin).

Conclusion

The ASP microarray has potential in the clinic as a diagnostic tool, as a research tool for both known and emerging pathogens, and as an early warning system for pathogenic bacteria that have been recently modified either naturally or deliberately.

Gene expression profiling in porcine maternal infanticide: a model for puerperal psychosis

The etiology of mental disorders remains largely unclear. Complex interactions between genetic and environmental factors are key to the development of such disorders. Puerperal psychosis is the most extreme form of postnatal mood disorder in women. Similarly, parturition in the pig can trigger extreme behavioral disturbances, including maternal infanticide. In this study, we have used a targeted cDNA microarray approach using the pig as a model to understand the genes and genetic pathways that are involved in these processes. Two subtracted cDNA libraries from porcine hypothalamus were constructed, which were enriched for genes that were over-expressed and under-expressed in the aberrant behavioral phenotype, compared to the matched control. In addition to this, a normalized library was constructed from hypothalamus and pituitary samples taken from pigs in a variety of reproductive states. The libraries were partially sequenced and combined represented approximately 5,159 different genes. Microarray analysis determined differences in gene expression between hypothalamus samples from nine matched pairs of infanticidal versus control animals, using a common reference design. Microarray analysis of variance (MAANOVA) identified 52 clones as being differentially expressed (P <or= 0.002) in the infanticide phenotype, a second analysis with friendly statistics package for microarray analysis (FSPMA) identified 9 genes in common to MAANOVA, and a further 16 genes. A rapid cross-species screen onto a human oligonucleotide array confirmed 3 genes and highlighted 61 more potential candidates. Some of these genes and the pathways in which they are involved were also implicated in a parallel QTL study on maternal infanticide.

Cold hardening processes in the Antarctic springtail, Cryptopygus antarcticus: clues from a microarray

The physiology of the Antarctic microarthropod, Cryptopygus antarcticus, has been well studied, particularly with regard to its ability to withstand low winter temperatures. However, the molecular mechanisms underlying this phenomenon are still poorly understood. 1180 sequences (Expressed Sequence Tags or ESTs) were generated and analysed, from populations of C. antarcticus. This represents the first publicly available sequence data for this species. A sub-set (672 clones) were used to generate a small microarray to examine the differences in gene expression between summer acclimated cold tolerant and non-cold tolerant springtails. Although 60% of the clones showed no sequence similarity to annotated genes in the datasets, of those where putative function could be inferred via database homology, there was a clear pattern of up-regulation of structural proteins being associated with the cold tolerant group. These structural proteins mainly comprised cuticle proteins and provide support for the recent theory that summer SCP variation within Collembola species could be a consequence of moulting, with moulting population having lowered SCPs.

Microarray based analysis of temperature and oxidative stress induced messenger RNA in Schistosoma mansoni

The body's defense against schistosome infection can take many forms. For example, upon developing acute schistosomiasis, patients often have fever coinciding with larval maturation, migration and early oviposition. As the infection becomes established, the parasite comes under oxidative stress generated by the host immune system. The most common treatment for schistosomiasis is the anti-helminthic drug praziquantel. Its effectiveness, however, is limited due to its inability to kill schistosomes 2-4 weeks post-infection. Clearly there is a need for new anti-schistosomal drugs. We hypothesize that gene products expressed as part of a protective response against heat and/or oxidative stress are potential therapeutic targets for future drug development. Using a 12,166 element oligonucleotide microarray to characterize Schistosoma mansoni genes induced by heat and oxidative stress we found that 1878 S. mansoni elements were significantly induced by heat stress. These included previously reported heat-shock genes expressing homologs of HSP40, HSP70 and HSP86. One thousand and one elements were induced by oxidative stress including those expressing homologs of superoxide dismutase, glutathione peroxidase and aldehyde dehydrogenase. Seventy-two elements were common to both stressors and could potentially be exploited in the development of novel anti-schistosomal therapeutics.

The trypanosome transcriptome is remodelled during differentiation but displays limited responsiveness within life stages

BACKGROUND

Trypanosomatids utilise polycistronic transcription for production of the vast majority of protein-coding mRNAs, which operates in the absence of gene-specific promoters. Resolution of nascent transcripts by polyadenylation and trans-splicing, together with specific rates of mRNA turnover, serve to generate steady state transcript levels that can differ in abundance across several orders of magnitude and can be developmentally regulated. We used a targeted oligonucleotide microarray, representing the strongly developmentally-regulated T. brucei membrane trafficking system and approximately 10% of the Trypanosoma brucei genome, to investigate both between-stage, or differentiation-dependent, transcriptome changes and within-stage flexibility in response to various challenges.

RESULTS

6% of the gene cohort are developmentally regulated, including several small GTPases, SNAREs, vesicle coat factors and protein kinases both consistent with and extending previous data. Therefore substantial differentiation-dependent remodeling of the trypanosome transcriptome is associated with membrane transport. Both the microarray and qRT-PCR were then used to analyse transcriptome changes resulting from specific gene over-expression, knockdown, altered culture conditions and chemical stress. Firstly, manipulation of Rab5 expression results in co-ordinate changes to clathrin protein expression levels and endocytotic activity, but no detectable changes to steady-state mRNA levels, which indicates that the effect is mediated post-transcriptionally. Secondly, knockdown of clathrin or the variant surface glycoprotein failed to perturb transcription. Thirdly, exposure to dithiothreitol or tunicamycin revealed no evidence for a classical unfolded protein response, mediated in higher eukaryotes by transcriptional changes. Finally, altered serum levels invoked little transcriptome alteration beyond changes to expression of ESAG6/7, the transferrin receptor.

CONCLUSION

While trypanosomes regulate mRNA abundance to effect the major changes accompanying differentiation, a given differentiated state appears transcriptionally inflexible. The implications of the absence of a transcriptome response in trypanosomes for both virulence and models of life cycle progression are discussed.

The T7-primer is a source of experimental bias and introduces variability between microarray platforms

Eberwine(-like) amplification of mRNA adds distinct 6–10 bp nucleotide stretches to the 5′ end of amplified RNA transcripts. Analysis of over six thousand microarrays reveals that probes containing motifs complementary to these stretches are associated with aberrantly high signals up to a hundred fold the signal observed in unaffected probes. This is not observed when total RNA is used as target source. Different T7 primer sequences are used in different laboratories and platforms and consequently different T7 primer bias is observed in different datasets. This will hamper efforts to compare data sets across platforms.

High-sensitivity detection methods for low-abundance RNA species: applications for functional genomics research

Gene expression analysis has facilitated a more complete understanding of the molecular biology of cellular processes and how variations of RNA expression are useful for the classification of various diseases. Furthermore, recent analysis of a variety of noncoding RNAs, such as microRNAs, has demonstrated that these RNAs play an important role in many cellular events, including cell differentiation and death, and may also serve as biological markers for disease. Besides helping in the understanding of diseases, RNA analysis is used in drug discovery, patient prognosis and treatment evaluation. One obstacle left to overcome is the amount of material required for the analysis, particularly when trying to extract information from precious, limited, clinical samples. Here we review the many approaches scientists take to either amplify the amount of RNA or amplify the signal generated from small amounts of RNA.

Integrating whole transcriptome assays on a lab-on-a-chip for single cell gene profiling

To correlate gene expression profiles to fundamental biological processes such as cell growth, differentiation and migration, it is essential to work at the single cell level. Gene expression analysis always starts with the relatively low efficient reverse transcription (RT) of RNA into complementary DNA (cDNA), an essential step as unprocessed RNAs will not be analysed further. In this paper, we present a novel method for RT that uses microfluidics to manipulate nanolitre volumes. We compare our method to conventional protocols performed in microlitre volumes. More specifically, reverse transcription was performed either in a polydimethylsiloxane (PDMS) rotary mixer or in a tube, using a single cell amount of mouse brain RNA (10 pg), and was followed by a template-switching PCR (TS-PCR) amplification step. We demonstrate that, using the microfluidic protocol, 74% of the genes expressed in mouse brain were detected, while only 4% were found with the conventional approach. We next profiled single neuronal progenitors. Using our microfluidic approach, i.e. performing cell capture, lysis and reverse transcription on-chip followed by TS-PCR amplification in tube, a mean of 5000 genes were detected in each neuron, which corresponds to the expected number of genes expressed in a single cell. This demonstrates the outstanding sensitivity of the microfluidic method.

Hypogonadal Mouse, a Model to Study the Effects of the Endogenous Lack of Gonadotropins on Apoptosis1

Testicular apoptosis is involved in the regulation of germ cell numbers, allowing optimal sperm production. Apoptosis has been described to occur in response to the absence of hormonal stimulation of the testis. Here we investigate the effect of the physiological lack of gonadotropins from birth using the hypogonadal (homozygous for the mutant allele Gnrh1(hpg)) mouse as a model. We pursued a concerted strategy using microarray analysis and RT-PCR to assess transcript levels, TUNEL to quantify the incidence of apoptosis, and Western blotting to assess the respective contribution of the extrinsic and intrinsic apoptotic pathways. Our results indicate a large increase in apoptosis of both somatic and germ cell compartments in the hpg testis, affecting Sertoli cells as well as germ cells of all ages. We confirmed our observations of Sertoli cell apoptosis using anti-Mullerian inhibiting substance staining and staining for cleaved fodrin alpha. In the somatic compartment, apoptosis is primarily regulated via the membrane receptor (extrinsic) apoptotic pathway, while in the germ cell compartment, regulation occurs via both the mitochondrial (intrinsic) and membrane receptor (extrinsic) apoptotic pathways, the latter potentially in a stage-specific manner. This study is the first report of spermatogonial apoptosis in response to gonadotropin deficiency as well as the first report of Sertoli cell apoptosis in response to gonadotropin deficiency in the mouse.