Probe selection for DNA microarrays using OligoWiz

Pleiotropic Roles for the Plasmodium berghei RNA Binding Protein UIS12 in Transmission and Oocyst Maturation

Colonization of the mosquito host by Plasmodium parasites is achieved by sexually differentiated gametocytes. Gametocytogenesis, gamete formation and fertilization are tightly regulated processes, and translational repression is a major regulatory mechanism for stage conversion. Here, we present a characterization of a Plasmodium berghei RNA binding protein, UIS12, that contains two conserved eukaryotic RNA recognition motifs (RRM). Targeted gene deletion resulted in viable parasites that replicate normally during blood infection, but form fewer gametocytes. Upon transmission to Anopheles stephensi mosquitoes, both numbers and size of midgut-associated oocysts were reduced and their development stopped at an early time point. As a consequence, no salivary gland sporozoites were formed indicative of a complete life cycle arrest in the mosquito vector. Comparative transcript profiling in mutant and wild-type infected red blood cells revealed a decrease in transcript abundance of mRNAs coding for signature gamete-, ookinete-, and oocyst-specific proteins in uis12(-) parasites. Together, our findings indicate multiple roles for UIS12 in regulation of gene expression after blood infection in good agreement with the pleiotropic defects that terminate successful sporogony and onward transmission to a new vertebrate host.

Multiplexed primer extension sequencing: A targeted RNA-seq method that enables high-precision quantitation of mRNA splicing isoforms and rare pre-mRNA splicing intermediates

The study of pre-mRNA splicing has been greatly aided by the advent of RNA sequencing (RNA-seq), which enables the genome-wide detection of discrete splice isoforms. Quantification of these splice isoforms requires analysis of splicing informative sequencing reads, those that unambiguously map to a single splice isoform, including exon-intron spanning alignments corresponding to retained introns, as well as exon-exon junction spanning alignments corresponding to either canonically- or alternatively-spliced isoforms. Because most RNA-seq experiments are designed to produce sequencing alignments that uniformly cover the entirety of transcripts, only a comparatively small number of splicing informative alignments are generated for any given splice site, leading to a decreased ability to detect and/or robustly quantify many splice isoforms. To address this problem, we have recently described a method termed Multiplexed Primer Extension sequencing, or MPE-seq, which uses pools of reverse transcription primers to target sequencing to user selected loci. By targeting reverse transcription to pre-mRNA splice junctions, this approach enables a dramatic enrichment in the fraction of splicing informative alignments generated per splicing event, yielding an increase in both the precision with which splicing efficiency can be measured, and in the detection of splice isoforms including rare splicing intermediates. Here we provide a brief review of the shortcomings associated with RNA-seq that drove our development of MPE-seq, as well as a detailed protocol for implementation of MPE-seq.

Target-enrichment sequencing for detailed characterization of small RNAs

Identification of important, functional small RNA (sRNA) species is currently hampered by the lack of reliable and sensitive methods to isolate and characterize them. We have developed a method, termed target-enrichment of sRNAs (TEsR), that enables targeted sequencing of rare sRNAs and diverse precursor and mature forms of sRNAs not detectable by current standard sRNA sequencing methods. It is based on the amplification of full-length sRNA molecules, production of biotinylated RNA probes, hybridization to one or multiple targeted RNAs, removal of nontargeted sRNAs and sequencing. By this approach, target sRNAs can be enriched by a factor of 500-30,000 while maintaining strand specificity. TEsR enriches for sRNAs irrespective of length or different molecular features, such as the presence or absence of a 5' cap or of secondary structures or abundance levels. Moreover, TEsR allows the detection of the complete sequence (including sequence variants, and 5' and 3' ends) of precursors, as well as intermediate and mature forms, in a quantitative manner. A well-trained molecular biologist can complete the TEsR procedure, from RNA extraction to sequencing library preparation, within 4-6 d.

Phospho-site mutants of the RNA Polymerase II C-terminal domain alter subtelomeric gene expression and chromatin modification state in fission yeast

Eukaryotic gene expression requires that RNA Polymerase II (RNAP II) gain access to DNA in the context of chromatin. The C-terminal domain (CTD) of RNAP II recruits chromatin modifying enzymes to promoters, allowing for transcription initiation or repression. Specific CTD phosphorylation marks facilitate recruitment of chromatin modifiers, transcriptional regulators, and RNA processing factors during the transcription cycle. However, the readable code for recruiting such factors is still not fully defined and how CTD modifications affect related families of genes or regional gene expression is not well understood. Here, we examine the effects of manipulating the Y₁S₂P₃T₄S₅P₆S₇ heptapeptide repeat of the CTD of RNAP II in Schizosaccharomyces pombe by substituting non-phosphorylatable alanines for Ser2 and/or Ser7 and the phosphomimetic glutamic acid for Ser7. Global gene expression analyses were conducted using splicing-sensitive microarrays and validated via RT-qPCR. The CTD mutations did not affect pre-mRNA splicing or snRNA levels. Rather, the data revealed upregulation of subtelomeric genes and alteration of the repressive histone H3 lysine 9 methylation (H3K9me) landscape. The data further indicate that H3K9me and expression status are not fully correlated, suggestive of CTD-dependent subtelomeric repression mechansims that act independently of H3K9me levels.

Multiplexed Intact-Tissue Transcriptional Analysis at Cellular Resolution

In recently developed approaches for high-resolution imaging within intact tissue, molecular characterization over large volumes has been largely restricted to labeling of proteins. But volumetric nucleic acid labeling may represent a far greater scientific and clinical opportunity, enabling detection of not only diverse coding RNA variants but also non-coding RNAs. Moreover, scaling immunohistochemical detection to large tissue volumes has limitations due to high cost, limited renewability/availability, and restricted multiplexing capability of antibody labels. With the goal of versatile, high-content, and scalable molecular phenotyping of intact tissues, we developed a method using carbodiimide-based chemistry to stably retain RNAs in clarified tissue, coupled with amplification tools for multiplexed detection. The resulting technology enables robust measurement of activity-dependent transcriptional signatures, cell-identity markers, and diverse non-coding RNAs in rodent and human tissue volumes. The growing set of validated probes is deposited in an online resource for nucleating related developments from across the scientific community.

Adaptation of Lactococcus lactis to high growth temperature leads to a dramatic increase in acidification rate

Lactococcus lactis is essential for most cheese making, and this mesophilic bacterium has its growth optimum around 30 °C. We have, through adaptive evolution, isolated a mutant TM29 that grows well up to 39 °C, and continuous growth at 40 °C is possible if pre-incubated at a slightly lower temperature. At the maximal permissive temperature for the wild-type, 38 °C, TM29 grows 33% faster and has a 12% higher specific lactate production rate than its parent MG1363, which results in fast lactate accumulation. Genome sequencing was used to reveal the mutations accumulated, most of which were shown to affect thermal tolerance. Of the mutations with more pronounced effects, two affected expression of single proteins (chaperone; riboflavin transporter), two had pleiotropic effects (RNA polymerase) which changed the gene expression profile, and one resulted in a change in the coding sequence of CDP-diglyceride synthase. A large deletion containing 10 genes was also found to affect thermal tolerance significantly. With this study we demonstrate a simple approach to obtain non-GMO derivatives of the important L. lactis that possess properties desirable by the industry, e.g. thermal robustness and increased rate of acidification. The mutations we have identified provide a genetic basis for further investigation of thermal tolerance.

Development of a DNA-based microarray for the detection of zoonotic pathogens in rodent species

The demand for diagnostic tools that allow simultaneous screening of samples for multiple pathogens is increasing because they overcome the limitations of other methods, which can only screen for a single or a few pathogens at a time. Microarrays offer the advantages of being capable to test a large number of samples simultaneously, screening for multiple pathogen types per sample and having comparable sensitivity to existing methods such as PCR. Array design is often considered the most important process in any microarray experiment and can be the deciding factor in the success of a study. There are currently no microarrays for simultaneous detection of rodent-borne pathogens. The aim of this report is to explicate the design, development and evaluation of a microarray platform for use as a screening tool that combines ease of use and rapid identification of a number of rodent-borne pathogens of zoonotic importance. Nucleic acid was amplified by multiplex biotinylation PCR prior to hybridisation onto microarrays. The array sensitivity was comparable to standard PCR, though less sensitive than real-time PCR. The array presented here is a prototype microarray identification system for zoonotic pathogens that can infect rodent species.

•

We have developed a microarray to detect zoonotic pathogens in rodent species.

•

The design stage of a microarray experiment is crucial for a successful experiment.

•

We examined the difference between amplification methods prior to hybridisation.

Genome-Wide Transcriptional Profiling of Clostridium perfringens SM101 during Sporulation Extends the Core of Putative Sporulation Genes and Genes Determining Spore Properties and Germination Characteristics

The formation of bacterial spores is a highly regulated process and the ultimate properties of the spores are determined during sporulation and subsequent maturation. A wide variety of genes that are expressed during sporulation determine spore properties such as resistance to heat and other adverse environmental conditions, dormancy and germination responses. In this study we characterized the sporulation phases of C. perfringens enterotoxic strain SM101 based on morphological characteristics, biomass accumulation (OD600), the total viable counts of cells plus spores, the viable count of heat resistant spores alone, the pH of the supernatant, enterotoxin production and dipicolinic acid accumulation. Subsequently, whole-genome expression profiling during key phases of the sporulation process was performed using DNA microarrays, and genes were clustered based on their time-course expression profiles during sporulation. The majority of previously characterized C. perfringens germination genes showed upregulated expression profiles in time during sporulation and belonged to two main clusters of genes. These clusters with up-regulated genes contained a large number of C. perfringens genes which are homologs of Bacillus genes with roles in sporulation and germination; this study therefore suggests that those homologs are functional in C. perfringens. A comprehensive homology search revealed that approximately half of the upregulated genes in the two clusters are conserved within a broad range of sporeforming Firmicutes. Another 30% of upregulated genes in the two clusters were found only in Clostridium species, while the remaining 20% appeared to be specific for C. perfringens. These newly identified genes may add to the repertoire of genes with roles in sporulation and determining spore properties including germination behavior. Their exact roles remain to be elucidated in future studies.

Investigating nanoparticle properties in plasmonic nanoarchitectures with DNA by surface plasmon resonance imaging

The optimal optical coupling conditions between SPWs and nanoparticle LSPs can be achieved by overlapping the source wavelength with the wavelength of excitation of LSPs.

A comparison of key aspects of gene regulation in Streptomyces coelicolor and Escherichia coli using nucleotide-resolution transcription maps produced in parallel by global and differential RNA sequencing

Streptomyces coelicolor is a model for studying bacteria renowned as the foremost source of natural products used clinically. Post-genomic studies have revealed complex patterns of gene expression and links to growth, morphological development and individual genes. However, the underlying regulation remains largely obscure, but undoubtedly involves steps after transcription initiation. Here we identify sites involved in RNA processing and degradation as well as transcription within a nucleotide-resolution map of the transcriptional landscape. This was achieved by combining RNA-sequencing approaches suited to the analysis of GC-rich organisms. Escherichia coli was analysed in parallel to validate the methodology and allow comparison. Previously, sites of RNA processing and degradation had not been mapped on a transcriptome-wide scale for E. coli. Through examples, we show the value of our approach and data sets. This includes the identification of new layers of transcriptional complexity associated with several key regulators of secondary metabolism and morphological development in S. coelicolor and the identification of host-encoded leaderless mRNA and rRNA processing associated with the generation of specialized ribosomes in E. coli. New regulatory small RNAs were identified for both organisms. Overall the results illustrate the diversity in mechanisms used by different bacterial groups to facilitate and regulate gene expression.

Bacillus pumilus reveals a remarkably high resistance to hydrogen peroxide provoked oxidative stress

Bacillus pumilus is characterized by a higher oxidative stress resistance than other comparable industrially relevant Bacilli such as B. subtilis or B. licheniformis. In this study the response of B. pumilus to oxidative stress was investigated during a treatment with high concentrations of hydrogen peroxide at the proteome, transcriptome and metabolome level. Genes/proteins belonging to regulons, which are known to have important functions in the oxidative stress response of other organisms, were found to be upregulated, such as the Fur, Spx, SOS or CtsR regulon. Strikingly, parts of the fundamental PerR regulon responding to peroxide stress in B. subtilis are not encoded in the B. pumilus genome. Thus, B. pumilus misses the catalase KatA, the DNA-protection protein MrgA or the alkyl hydroperoxide reductase AhpCF. Data of this study suggests that the catalase KatX2 takes over the function of the missing KatA in the oxidative stress response of B. pumilus. The genome-wide expression analysis revealed an induction of bacillithiol (Cys-GlcN-malate, BSH) relevant genes. An analysis of the intracellular metabolites detected high intracellular levels of this protective metabolite, which indicates the importance of bacillithiol in the peroxide stress resistance of B. pumilus.

Discrimination of single base mismatched oligonucleotides related to the rpoB gene of Mycobacterium tuberculosis using a surface plasmon resonance biosensor

Single base mismatched oligonucleotides related to the rpoB gene of Mycobacterium tuberculosis, the mutations of which cause drug resistance of the infectious agent, were detected and discriminated using a surface plasmon resonance biosensor system. Thiol-modified oligonucleotides of the selected sequence (the probe) and 1-mercapto-6-hexanol were immobilized on a gold sensor surface. Hybridization between immobilized probe P2 and perfectly matched target T2 as well as a single base mismatched target TN was investigated in buffer solutions of various stringencies. Discrimination of perfectly matched and single base mismatched targets is achieved due to a difference in the level of their hybridization with the immobilized probe depending on stringency of the buffer solution. In 0.5×SSC buffer solution (7.5 mM sodium citrate, pH 7, containing 75 mM NaCl), sensor response at T2 injection into the measuring sensor cell was 16 times that at TN injection. The experimental results on surface hybridization between the studied oligonucleotides demonstrated a good correlation with theoretical calculations of thermodynamic parameters of these interactions in the solution. The described approach could be proposed as a basis for creating a biosensor for real-time label-free diagnostics of drug-resistant tuberculosis.

Osmotic regulation of hepatic betaine metabolism

Betaine critically contributes to the control of hepatocellular hydration and provides protection of the liver from different kinds of stress. To investigate how the hepatocellular hydration state affects gene expression of enzymes involved in the metabolism of betaine and related organic osmolytes, we used quantitative RT-PCR gene expression studies in rat hepatoma cells as well as metabolic and gene expression profiling in primary hepatocytes of both wild-type and 5,10-methylenetetrahydrofolate reductase (MTHFR)-deficient mice. Anisotonic incubation caused coordinated adaptive changes in the expression of various genes involved in betaine metabolism, in particular of betaine homocysteine methyltransferase, dimethylglycine dehydrogenase, and sarcosine dehydrogenase. The expression of betaine-degrading enzymes was downregulated by cell shrinking and strongly induced by an increase in cell volume under hypotonic conditions. Metabolite concentrations in the culture system changed accordingly. Expression changes were mediated through tyrosine kinases, cyclic nucleotide-dependent protein kinases, and JNK-dependent signaling. Assessment of hepatic gene expression using a customized microarray chip showed that hepatic betaine depletion in MTHFR(-/-) mice was associated with alterations that were comparable to those induced by cell swelling in hepatocytes. In conclusion, the adaptation of hepatocytes to changes in cell volume involves the coordinated regulation of betaine synthesis and degradation and concomitant changes in intracellular osmolyte concentrations. The existence of such a well-orchestrated response underlines the importance of cell volume homeostasis for liver function and of methylamine osmolytes such as betaine as hepatic osmolytes.

Combining bulk segregation analysis and microarrays for mapping of the pH trait in melon

The availability of sequence information for many plants has opened the way to advanced genetic analysis in many non-model plants. Nevertheless, exploration of genetic variation on a large scale and its use as a tool for the identification of traits of interest are still rare. In this study, we combined a bulk segregation approach with our own-designed microarrays to map the pH locus that influences fruit pH in melon. Using these technologies, we identified a set of markers that are genetically linked to the pH trait. Further analysis using a set of melon cultivars demonstrated that some of these markers are tightly linked to the pH trait throughout our germplasm collection. These results validate the utility of combining microarray technology with a bulk segregation approach in mapping traits of interest in non-model plants.

Surface plasmon resonance imaging for nucleic acid detection

Surface plasmon resonance imaging (SPRI) is a powerful tool for simple, fast and cheap nucleic acid detection. Great efforts have been made during the last decade with the aim of developing even more sensitive and specific SPRI-based methods to be used for the direct detection of DNA and RNA. Here, after a description of the fundamentals of SPRI, the state of the art of recent platform and assay developments is presented, with special attention given to advances in SPRI signal enhancement procedures.

Single nucleotide polymorphism detection by optical DNA-based sensing coupled with whole genomic amplification

The work presented here deals with the optimization of a strategy for detection of single nucleotide polymorphisms based on surface plasmon resonance imaging. First, a sandwich-like assay was designed, and oligonucleotide sequences were computationally selected in order to study optimized conditions for the detection of the rs1045642 single nucleotide polymorphism in the gene ABCB1. Then the strategy was optimized on a surface plasmon resonance imaging biosensor using synthetic DNA sequences in order to evaluate the best conditions for the detection of a single mismatching base. Finally, the assay was tested on DNA extracted from human blood which was subsequently amplified using a whole genome amplification kit. The direct detection of the polymorphism was successfully achieved. The biochip was highly regenerable and reusable for up to 20 measurements. Furthermore, coupling these promising results with the multiarray assay, we can foresee applying this biosensor in clinical research extended to concurrent analysis of different polymorphisms.

Divergent gene expression in the conserved dauer stage of the nematodes Pristionchus pacificus and Caenorhabditis elegans

Background

An organism can respond to changing environmental conditions by adjusting gene regulation and by forming alternative phenotypes. In nematodes, these mechanisms are coupled because many species will form dauer larvae, a stress-resistant and non-aging developmental stage, when exposed to unfavorable environmental conditions, and execute gene expression programs that have been selected for the survival of the animal in the wild. These dauer larvae represent an environmentally induced, homologous developmental stage across many nematode species, sharing conserved morphological and physiological properties. Hence it can be expected that some core components of the associated transcriptional program would be conserved across species, while others might diverge over the course of evolution. However, transcriptional and metabolic analysis of dauer development has been largely restricted to Caenorhabditis elegans. Here, we use a transcriptomic approach to compare the dauer stage in the evolutionary model system Pristionchus pacificus with the dauer stage in C. elegans.

Results

We have employed Agilent microarrays, which represent 20,446 P. pacificus and 20,143 C. elegans genes to show an unexpected divergence in the expression profiles of these two nematodes in dauer and dauer exit samples. P. pacificus and C. elegans differ in the dynamics and function of genes that are differentially expressed. We find that only a small number of orthologous gene pairs show similar expression pattern in the dauers of the two species, while the non-orthologous fraction of genes is a major contributor to the active transcriptome in dauers. Interestingly, many of the genes acquired by horizontal gene transfer and orphan genes in P. pacificus, are differentially expressed suggesting that these genes are of evolutionary and functional importance.

Conclusion

Our data set provides a catalog for future functional investigations and indicates novel insight into evolutionary mechanisms. We discuss the limited conservation of core developmental and transcriptional programs as a common aspect of animal evolution.

A genomic bias for genotype-environment interactions in C. elegans

The phenotype of an organism is determined by its genotype and environment. An interaction between these two arises from the differential effect of the environment on gene expression in distinct genotypes; however, the genomic properties identifying these are not well understood. Here we analyze the transcriptomes of five C. elegans strains (genotype) cultivated in five growth conditions (environment), and find that highly regulated genes, as distinguished by intergenic lengths, motif concentration, and expression levels, are particularly biased toward genotype-environment interactions. Sequencing these strains, we find that genes with expression variation across genotypes are enriched for promoter single-nucleotide polymorphisms (SNPs), as expected. However, genes with genotype-environment interactions do not significantly differ from background in terms of their promoter SNPs. Collectively, these results indicate that the highly regulated nature of particular genes predispose them for exhibiting genotype-environment interaction as a consequence of changes to upstream regulators. This observation may provide a deeper understanding into the origin of the extraordinary gene expression diversity present in even closely related species.

Developmental milestones punctuate gene expression in the Caenorhabditis embryo

A fundamental question in developmental biology relates to the connection between morphological stages and their underlying molecular activity. Here we demonstrate that, at the molecular level, embryonic development in five Caenorhabditis species proceeds through two distinct milestones in which the transcriptome is resistant to differences in species-specific developmental timings. By comparing the complete protein-coding transcriptomes of individually timed embryos across ten morphological markers, we found that developmental milestones can be characterized by their expression dynamics and activation of key developmental regulators. This approach led us to discover the nematode phylotypic stage and to show that in chordates and arthropods it is represented as two distinct stages, suggesting that animal body plans might evolve by uncoupling and elaboration on formerly synchronous processes.

Self-custom-made SFP arrays for nonmodel organisms

Successful genetic mapping is dependent upon a high-density set of markers. Therefore, tools for high-throughput discovery of genetic variation are essential. The most abundant genetic marker is the single-nucleotide polymorphism (SNP). However, except for model organisms, genomic information is still limited. Although high-throughput genomic sequencing technologies are becoming relatively inexpensive, only low-throughput genetic markers are accessible (e.g., simple sequence repeats). The use of sequencing for the discovery and screening of high-density genetic variation in whole populations is still expensive. Alternatively, hybridization of genomic DNA (gDNA) on a reference (either genome or transcriptome) is an efficient approach for genetic screening without knowing the alleles in advance (Borevitz et al. Proc Natl Acad Sci USA 104:12057-12062). We describe a protocol for the design of probes for a high-throughput genetic-marker discovery microarray, termed single feature polymorphism (SFP) array. Starting with consensus cDNA sequences (UniGenes), we use OligoWiz to design T (m)-optimized 50-bp long oligonucleotide probes (Ophir et al. BMC Genomics 11:269, 2010). This design is similar to expression arrays and we point out the differences.

Transcriptome dynamics of Pseudomonas putida KT2440 under water stress

Water deprivation can be a major stressor to microbial life in surface and subsurface soil. In unsaturated soils, the matric potential (Ψ(m)) is often the main component of the water potential, which measures the thermodynamic availability of water. A low matric potential usually translates into water forming thin liquid films in the soil pores. Little is known of how bacteria respond to such conditions, where, in addition to facing water deprivation that might impair their metabolism, they have to adapt their dispersal strategy as swimming motility may be compromised. Using the pressurized porous surface model (PPSM), which allows creation of thin liquid films by controlling Ψ(m), we examined the transcriptome dynamics of Pseudomonas putida KT2440. We identified the differentially expressed genes in cells exposed to a mild matric stress (-0.4 MPa) for 4, 24, or 72 h. The major response was detected at 4 h before gradually disappearing. Upregulation of alginate genes was notable in this early response. Flagellar genes were not downregulated, and the microarray data even suggested increasing expression as the stress prolonged. Moreover, we tested the effect of polyethylene glycol 8000 (PEG 8000), a nonpermeating solute often used to simulate Ψ(m), on the gene expression profile and detected a different profile than that observed by directly imposing Ψ(m). This study is the first transcriptome profiling of KT2440 under directly controlled Ψ(m) and also the first to show the difference in gene expression profiles between a PEG 8000-simulated and a directly controlled Ψ(m).

The effects of mitochondrial genotype on hypoxic survival and gene expression in a hybrid population of the killifish, Fundulus heteroclitus

The physiological link between oxygen availability and mitochondrial function is well established. However, whether or not fitness variation is associated with mitochondrial genotypes in the field remains a contested topic in evolutionary biology. In this study, we draw on a population of the teleost fish, Fundulus heteroclitus, where functionally distinct subspecies hybridize, likely as a result of past glacial events. We had two specific aims: (i) to determine the effect of mtDNA genotype on survivorship of male and female fish under hypoxic stress and (ii) to determine the effect of hypoxic stress, sex and mtDNA genotype on gene expression. We found an unexpected and highly significant effect of sex on survivorship under hypoxic conditions, but no significant effect of mtDNA genotype. Gene expression analyses revealed hundreds of transcripts differentially regulated by sex and hypoxia. Mitochondrial transcripts and other predicted pathways were among those influenced by hypoxic stress, and a transcript corresponding to the mtDNA control region was the most highly suppressed transcript under the conditions of hypoxia. An RT-PCR experiment on the control region was consistent with microarray results. Effects of mtDNA sequence variation on genome expression were limited; however, a potentially important epistasis between mtDNA sequence and expression of a nuclear-encoded mitochondrial translation protein was discovered. Overall, these results confirm that mitochondrial regulation is a major component of hypoxia tolerance and further suggest that purifying selection has been the predominant selective force on mitochondrial genomes in these two subspecies.

Comparative genomics and transcriptomics of Propionibacterium acnes

The anaerobic gram-positive bacterium Propionibacterium acnes is a human skin commensal that is occasionally associated with inflammatory diseases. Recent work has indicated that evolutionary distinct lineages of P. acnes play etiologic roles in disease while others are associated with maintenance of skin homeostasis. To shed light on the molecular basis for differential strain properties, we carried out genomic and transcriptomic analysis of distinct P. acnes strains. We sequenced the genome of the P. acnes strain 266, a type I-1a strain. Comparative genome analysis of strain 266 and four other P. acnes strains revealed that overall genome plasticity is relatively low; however, a number of island-like genomic regions, encoding a variety of putative virulence-associated and fitness traits differ between phylotypes, as judged from PCR analysis of a collection of P. acnes strains. Comparative transcriptome analysis of strains KPA171202 (type I-2) and 266 during exponential growth revealed inter-strain differences in gene expression of transport systems and metabolic pathways. In addition, transcript levels of genes encoding possible virulence factors such as dermatan-sulphate adhesin, polyunsaturated fatty acid isomerase, iron acquisition protein HtaA and lipase GehA were upregulated in strain 266. We investigated differential gene expression during exponential and stationary growth phases. Genes encoding components of the energy-conserving respiratory chain as well as secreted and virulence-associated factors were transcribed during the exponential phase, while the stationary growth phase was characterized by upregulation of genes involved in stress responses and amino acid metabolism. Our data highlight the genomic basis for strain diversity and identify, for the first time, the actively transcribed part of the genome, underlining the important role growth status plays in the inflammation-inducing activity of P. acnes. We argue that the disease-causing potential of different P. acnes strains is not only determined by the phylotype-specific genome content but also by variable gene expression.

Universal oligonucleotide microarray for sub-typing of Influenza A virus

A universal microchip was developed for genotyping Influenza A viruses. It contains two sets of oligonucleotide probes allowing viruses to be classified by the subtypes of hemagglutinin (H1-H13, H15, H16) and neuraminidase (N1-N9). Additional sets of probes are used to detect H1N1 swine influenza viruses. Selection of probes was done in two steps. Initially, amino acid sequences specific to each subtype were identified, and then the most specific and representative oligonucleotide probes were selected. Overall, between 19 and 24 probes were used to identify each subtype of hemagglutinin (HA) and neuraminidase (NA). Genotyping included preparation of fluorescently labeled PCR amplicons of influenza virus cDNA and their hybridization to microarrays of specific oligonucleotide probes. Out of 40 samples tested, 36 unambiguously identified HA and NA subtypes of Influenza A virus.

[Oligonucleotide microarray for subtyping of influenza virus A neuraminidase]

Microarray for influenza A neuraminidase subtyping was presented. Selection of oligoprobes proceeded in two steps. First step included selection of peptides specific for each subtype of neuraminidase. At the second step oligoprobes were calculated using found peptides structures with the subsequent additional selection of the most specific and representative probes. From 19 to 24 probes were used for determination of each subtype of neuraminidase. Microchip testing for 19 samples with the most widespread types (N1 and N2) specifies in unequivocal definition 18 of them and only one isolate has not been identified.

Cost-effective multiplexing before capture allows screening of 25 000 clinically relevant SNPs in childhood acute lymphoblastic leukemia

Genetic variants, including single-nucleotide polymorphisms (SNPs), are key determiners of interindividual differences in treatment efficacy and toxicity in childhood acute lymphoblastic leukemia (ALL). Although up to 13 chemotherapeutic agents are used in the treatment of this cancer, it remains a model disease for exploring the impact of genetic variation due to well-characterized cytogenetics, drug response pathways and precise monitoring of minimal residual disease. Here, we have selected clinically relevant genes and SNPs through literature screening, and on the basis of associations with key pathways, protein-protein interactions or downstream partners that have a role in drug disposition and treatment efficacy in childhood ALL. This allows exploration of pathways, where one of several genetic variants may lead to similar clinical phenotypes through related molecular mechanisms. We have designed a cost-effective, high-throughput capture assay of ∼25,000 clinically relevant SNPs, and demonstrated that multiple samples can be tagged and pooled before genome capture in targeted enrichment with a sufficient sequencing depth for genotyping. This multiplexed, targeted sequencing method allows exploration of the impact of pharmacogenetics on efficacy and toxicity in childhood ALL treatment, which will be of importance for personalized chemotherapy.

Comparative profiling of the transcriptional response to iron restriction in six serotypes of Actinobacillus pleuropneumoniae with different virulence potential

Background

Comparative analysis of gene expression among serotypes within a species can provide valuable information on important differences between related genomes. For the pig lung pathogen Actinobacillus pleuropneumoniae, 15 serotypes with a considerable variation in virulence potential and immunogenicity have been identified. This serotypic diversity can only partly be explained by amount of capsule and differences in the RTX toxin genes in their genomes. Iron acquisition in vivo is an important bacterial function and in pathogenic bacteria, iron-limitation is often a signal for the induction of virulence genes. We used a pan-genomic microarray to study the transcriptional response to iron restriction in vitro in six serotypes of A. pleuropneumoniae (1, 2, 3, 5b, 6, and 7), representing at least two levels of virulence.

Results

In total, 45 genes were significantly (p < 0.0001) up-regulated and 67 genes significantly down-regulated in response to iron limitation. Not previously observed in A. pleuropneumoniae was the up-regulation of a putative cirA-like siderophore in all six serotypes. Three genes, recently described in A. pleuropneumoniae as possibly coding for haemoglobin-haptoglobin binding proteins, displayed significant serotype related up-regulation to iron limitation. For all three genes, the expression appeared at its lowest in serotype 3, which is generally considered one of the least virulent serotypes of A. pleuropneumoniae. The three genes share homology with the hmbR haemoglobin receptor of Neisseria meningitidis, a possible virulence factor which contributes to bacterial survival in rats.

Conclusions

By comparative analysis of gene expression among 6 different serotypes of A. pleuropneumoniae we identified a common set of presumably essential core genes, involved in iron regulation. The results support and expand previous observations concerning the identification of new potential iron acquisition systems in A. pleuropneumoniae, showing that this bacterium has evolved several strategies for scavenging the limited iron resources of the host. The combined effect of iron-depletion and serotype proved to be modest, indicating that serotypes of both moderate and high virulence at least in vitro are reacting almost identical to iron restriction. One notable exception, however, is the haemoglobin-haptoglobin binding protein cluster which merits further investigation.

Comparative transcriptomics implicates mechanisms of evolved pollution tolerance in a killifish population

Wild populations of the killifish Fundulus heteroclitus resident in heavily contaminated North American Atlantic coast estuaries have recently and independently evolved dramatic, heritable, and adaptive pollution tolerance. We compared physiological and transcriptome responses to embryonic polychlorinated biphenyl (PCB) exposures between one tolerant population and a nearby sensitive population to gain insight into genomic, physiological and biochemical mechanisms of evolved tolerance in killifish, which are currently unknown. The PCB exposure concentrations at which developmental toxicity emerged, the range of developmental abnormalities exhibited, and global as well as specific gene expression patterns were profoundly different between populations. In the sensitive population, PCB exposures produced dramatic, dose-dependent toxic effects, concurrent with the alterations in the expression of many genes. For example, PCB-mediated cardiovascular system failure was associated with the altered expression of cardiomyocyte genes, consistent with sarcomere mis-assembly. In contrast, genome-wide expression was comparatively refractory to PCB induction in the tolerant population. Tolerance was associated with the global blockade of the aryl hydrocarbon receptor (AHR) signalling pathway, the key mediator of PCB toxicity, in contrast to the strong dose-dependent up-regulation of AHR pathway elements observed in the sensitive population. Altered regulation of signalling pathways that cross-talk with AHR was implicated as one candidate mechanism for the adaptive AHR signalling repression and the pollution tolerance that it affords. In addition to revealing mechanisms of PCB toxicity and tolerance, this study demonstrates the value of comparative transcriptomics to explore molecular mechanisms of stress response and evolved adaptive differences among wild populations.

Expectations, validity, and reality in gene expression profiling

OBJECTIVE

To provide a critical overview of gene expression profiling methodology and discuss areas of future development.

RESULTS

Gene expression profiling has been used extensively in biological research and has resulted in significant advances in the understanding of the molecular mechanisms of complex disorders, including cancer, heart disease, and metabolic disorders. However, translating this technology into genomic medicine for use in diagnosis and prognosis faces many challenges. In addition, gene expression profile analysis is frequently controversial, because its conclusions often lack reproducibility and claims of effective dissemination into translational medicine have, in some cases, been remarkably unjustified. In the last decade, a large number of methodological and technical solutions have been offered to overcome the challenges.

STUDY DESIGN AND SETTING

We consider the strengths, limitations, and appropriate applications of gene expression profiling techniques, with particular reference to the clinical relevance.

CONCLUSION

Some studies have demonstrated the ability and clinical utility of gene expression profiling for use as diagnostic, prognostic, and predictive molecular markers. The challenges of gene expression profiling lie with the standardization of analytic approaches and the evaluation of the clinical merit in broader heterogeneous populations by prospective clinical trials.

Functional genomics of physiological plasticity and local adaptation in killifish

Evolutionary solutions to the physiological challenges of life in highly variable habitats can span the continuum from evolution of a cosmopolitan plastic phenotype to the evolution of locally adapted phenotypes. Killifish (Fundulus sp.) have evolved both highly plastic and locally adapted phenotypes within different selective contexts, providing a comparative system in which to explore the genomic underpinnings of physiological plasticity and adaptive variation. Importantly, extensive variation exists among populations and species for tolerance to a variety of stressors, and we exploit this variation in comparative studies to yield insights into the genomic basis of evolved phenotypic variation. Notably, species of Fundulus occupy the continuum of osmotic habitats from freshwater to marine and populations within Fundulus heteroclitus span far greater variation in pollution tolerance than across all species of fish. Here, we explore how transcriptome regulation underpins extreme physiological plasticity on osmotic shock and how genomic and transcriptomic variation is associated with locally evolved pollution tolerance. We show that F. heteroclitus quickly acclimate to extreme osmotic shock by mounting a dramatic rapid transcriptomic response including an early crisis control phase followed by a tissue remodeling phase involving many regulatory pathways. We also show that convergent evolution of locally adapted pollution tolerance involves complex patterns of gene expression and genome sequence variation, which is confounded with body-weight dependence for some genes. Similarly, exploiting the natural phenotypic variation associated with other established and emerging model organisms is likely to greatly accelerate the pace of discovery of the genomic basis of phenotypic variation.

chipD: a web tool to design oligonucleotide probes for high-density tiling arrays

chipD is a web server that facilitates design of DNA oligonucleotide probes for high-density tiling arrays, which can be used in a number of genomic applications such as ChIP-chip or gene-expression profiling. The server implements a probe selection algorithm that takes as an input, in addition to the target sequences, a set of parameters that allow probe design to be tailored to specific applications, protocols or the array manufacturer's requirements. The algorithm optimizes probes to meet three objectives: (i) probes should be specific; (ii) probes should have similar thermodynamic properties; and (iii) the target sequence coverage should be homogeneous and avoid significant gaps. The output provides in a text format, the list of probe sequences with their genomic locations, targeted strands and hybridization characteristics. chipD has been used successfully to design tiling arrays for bacteria and yeast. chipD is available at http://chipd.uwbacter.org/.

High-throughput marker discovery in melon using a self-designed oligo microarray

BACKGROUND

Genetic maps constitute the basis of breeding programs for many agricultural organisms. The creation of these maps is dependent on marker discovery. Melon, among other crops, is still lagging in genomic resources, limiting the ability to discover new markers in a high-throughput fashion. One of the methods used to search for molecular markers is DNA hybridization to microarrays. Microarray hybridization of DNA from different accessions can reveal differences between them--single-feature polymorphisms (SFPs). These SFPs can be used as markers for breeding purposes, or they can be converted to conventional markers by sequencing. This method has been utilized in a few different plants to discover genetic variation, using Affymetrix arrays that exist for only a few organisms. We applied this approach with some modifications for marker discovery in melon.

RESULTS

Using a custom-designed oligonucleotide microarray based on a partial EST collection of melon, we discovered 6184 putative SFPs between the parents of our mapping population. Validation by sequencing of 245 SFPs from the two parents showed a sensitivity of around 79%. Most SFPs (81%) contained single-nucleotide polymorphisms. Testing the SFPs on another mapping population of melon confirmed that many of them are conserved.

CONCLUSION

Thousands of new SFPs that can be used for genetic mapping and molecular-assisted breeding in melon were discovered using a custom-designed oligo microarray. A portion of these SFPs are conserved and can be used in different breeding populations. Although improvement of the discovery rate is still needed, this approach is applicable to many agricultural systems with limited genomic resources.

A virulent parent with probiotic progeny: comparative genomics of Escherichia coli strains CFT073, Nissle 1917 and ABU 83972

Escherichia coli is a highly versatile species encompassing a diverse spectrum of strains, i.e. from highly virulent isolates causing serious infectious diseases to commensals and probiotic strains. Although much is known about bacterial pathogenicity in E. coli, the understanding of which genetic determinants differentiates a virulent from an avirulent strain still remains limited. In this study we designed a new comparative genomic hybridization microarray based on 31 sequenced E. coli strains and used it to compare two E. coli strains used as prophylactic agents (i.e. Nissle 1917 and 83972) with the highly virulent uropathogen CFT073. Only relatively minor genetic variations were found between the isolates, suggesting that the three strains may have originated from the same virulent ancestral parent. Interestingly, Nissle 1917 (a gut commensal strain) was more similar to CFT073 with respect to genotype and phenotype than 83972 (an asymptomatic bacteriuria strain). The study indicates that genetic variations (e.g. mutations) and expression differences, rather than genomic content per se, contribute to the divergence in disease-causing ability between these strains. This has implications for the use of virulence factors in epidemiological research, and emphasizes the need for more comparative genomic studies of closely related strains to compare their virulence potential.

Evaluating oligonucleotide properties for DNA microarray probe design

Most current microarray oligonucleotide probe design strategies are based on probe design factors (PDFs), which include probe hybridization free energy (PHFE), probe minimum folding energy (PMFE), dimer score, hairpin score, homology score and complexity score. The impact of these PDFs on probe performance was evaluated using four sets of microarray comparative genome hybridization (aCGH) data, which included two array manufacturing methods and the genomes of two species. Since most of the hybridizing DNA is equimolar in CGH data, such data are ideal for testing the general hybridization properties of almost all candidate oligonucleotides. In all our data sets, PDFs related to probe secondary structure (PMFE, hairpin score and dimer score) are the most significant factors linearly correlated with probe hybridization intensities. PHFE, homology and complexity score are correlating significantly with probe specificities, but in a non-linear fashion. We developed a new PDF, pseudo probe binding energy (PPBE), by iteratively fitting dinucleotide positional weights and dinucleotide stacking energies until the average residue sum of squares for the model was minimized. PPBE showed a better correlation with probe sensitivity and a better specificity than all other PDFs, although training data are required to construct a PPBE model prior to designing new oligonucleotide probes. The physical properties that are measured by PPBE are as yet unknown but include a platform-dependent component. A practical way to use these PDFs for probe design is to set cutoff thresholds to filter out bad quality probes. Programs and correlation parameters from this study are freely available to facilitate the design of DNA microarray oligonucleotide probes.

Genome-wide approaches to monitor pre-mRNA splicing

Pre-mRNA processing is an essential control-point in the gene expression pathway of eukaryotic organisms. The budding yeast Saccharomyces cerevisiae offers a powerful opportunity to examine the regulation of this pathway. In this chapter, we will describe methods that have been developed in our lab and others to examine pre-mRNA splicing from a genome-wide perspective in yeast. Our goal is to provide all of the necessary information--from microarray design to experimental setup to data analysis--to facilitate the widespread use of this technology.

Dynamic probe selection for studying microbial transcriptome with high-density genomic tiling microarrays

Background

Current commercial high-density oligonucleotide microarrays can hold millions of probe spots on a single microscopic glass slide and are ideal for studying the transcriptome of microbial genomes using a tiling probe design. This paper describes a comprehensive computational pipeline implemented specifically for designing tiling probe sets to study microbial transcriptome profiles.

Results

The pipeline identifies every possible probe sequence from both forward and reverse-complement strands of all DNA sequences in the target genome including circular or linear chromosomes and plasmids. Final probe sequence lengths are adjusted based on the maximal oligonucleotide synthesis cycles and best isothermality allowed. Optimal probes are then selected in two stages - sequential and gap-filling. In the sequential stage, probes are selected from sequence windows tiled alongside the genome. In the gap-filling stage, additional probes are selected from the largest gaps between adjacent probes that have already been selected, until a predefined number of probes is reached. Selection of the highest quality probe within each window and gap is based on five criteria: sequence uniqueness, probe self-annealing, melting temperature, oligonucleotide length, and probe position.

Conclusions

The probe selection pipeline evaluates global and local probe sequence properties and selects a set of probes dynamically and evenly distributed along the target genome. Unique to other similar methods, an exact number of non-redundant probes can be designed to utilize all the available probe spots on any chosen microarray platform. The pipeline can be applied to microbial genomes when designing high-density tiling arrays for comparative genomics, ChIP chip, gene expression and comprehensive transcriptome studies.

An oligonucleotide microarray to characterize multidrug resistant plasmids

Many of the Enterobacteriaceae carry multiple drug resistance (MDR) genes on large plasmids of replicon type Inc A/C and Inc H1. It is important to understand the transmission of these MDR plasmids because the genes they carry can affect the outcome of antimicrobial therapy. The aim of this study was to design a microarray with oligonucleotide probes for every gene in the six Inc A/C and one Inc H1 plasmids of interest while representing all redundant sequences only once. The microarray is printed in triplicate with 493 unique oligonucleotide probes 70 nucleotides in length. Salmonella enterica and Escherichia coli control strains and test plasmids (in the parent strain and transformed into a known E. coli background strain) were hybridized to the plasmid microarray. This hybridization arrays presents a rapid and cost effective method for high-density screening of isolates to evaluate the gene content of Inc A/C and H1 plasmids and will show how plasmids can change content with transmission.

Optimization of the BLASTN substitution matrix for prediction of non-specific DNA microarray hybridization

DNA microarray measurements are susceptible to error caused by non-specific hybridization between a probe and a target (cross-hybridization), or between two targets (bulk-hybridization). Search algorithms such as BLASTN can quickly identify potentially hybridizing sequences. We set out to improve BLASTN accuracy by modifying the substitution matrix and gap penalties. We generated gene expression microarray data for samples in which 1 or 10% of the target mass was an exogenous spike of known sequence. We found that the 10% spike induced 2-fold intensity changes in 3% of the probes, two-third of which were decreases in intensity likely caused by bulk-hybridization. These changes were correlated with similarity between the spike and probe sequences. Interestingly, even very weak similarities tended to induce a change in probe intensity with the 10% spike. Using this data, we optimized the BLASTN substitution matrix to more accurately identify probes susceptible to non-specific hybridization with the spike. Relative to the default substitution matrix, the optimized matrix features a decreased score for A–T base pairs relative to G–C base pairs, resulting in a 5–15% increase in area under the ROC curve for identifying affected probes. This optimized matrix may be useful in the design of microarray probes, and in other BLASTN-based searches for hybridization partners.

Web-based design of peptide microarrays using microPepArray Pro

Peptide microarrays (peptide arrays) have increasingly become an important research tool for studying protein detection, profiling, and protein-protein interactions, and they have the potential to foster high-throughput protein analysis as DNA arrays did for genomics research a decade ago. Recently, technologies have emerged that allow flexible synthesis of high-density peptide arrays based on specific application needs (e.g., phosphopeptide microarrays). To fully unleash the power of this promising research tool, significant efforts are required to develop computational and informatics resources that facilitate the experimental design and data analysis for a wide range of peptide array-based applications. The design of peptide arrays is inherently more complex than that of DNA arrays. We herein introduce microPepArray Pro, a Web-based general-purpose peptide array design program. microPepArray Pro features strong content design capabilities and maximized user control. The program suits the needs of a diversity of design tasks, works with a variety of peptide array configurations, and is highly expandable: new functionalities can be developed and added to microPepArray Pro with relative ease.

Custom design and analysis of high-density oligonucleotide bacterial tiling microarrays

Background

High-density tiling microarrays are a powerful tool for the characterization of complete genomes. The two major computational challenges associated with custom-made arrays are design and analysis. Firstly, several genome dependent variables, such as the genome's complexity and sequence composition, need to be considered in the design to ensure a high quality microarray. Secondly, since tiling projects today very often exceed the limits of conventional array-experiments, researchers cannot use established computer tools designed for commercial arrays, and instead have to redesign previous methods or create novel tools.

Principal Findings

Here we describe the multiple aspects involved in the design of tiling arrays for transcriptome analysis and detail the normalisation and analysis procedures for such microarrays. We introduce a novel design method to make two 280,000 feature microarrays covering the entire genome of the bacterial species Escherichia coli and Neisseria meningitidis, respectively, as well as the use of multiple copies of control probe-sets on tiling microarrays. Furthermore, a novel normalisation and background estimation procedure for tiling arrays is presented along with a method for array analysis focused on detection of short transcripts. The design, normalisation and analysis methods have been applied in various experiments and several of the detected novel short transcripts have been biologically confirmed by Northern blot tests.

Conclusions

Tiling-arrays are becoming increasingly applicable in genomic research, but researchers still lack both the tools for custom design of arrays, as well as the systems and procedures for analysis of the vast amount of data resulting from such experiments. We believe that the methods described herein will be a useful contribution and resource for researchers designing and analysing custom tiling arrays for both bacteria and higher organisms.

Model-based probe set optimization for high-performance microarrays

A major challenge in microarray design is the selection of highly specific oligonucleotide probes for all targeted genes of interest, while maintaining thermodynamic uniformity at the hybridization temperature. We introduce a novel microarray design framework (Thermodynamic Model-based Oligo Design Optimizer, TherMODO) that for the first time incorporates a number of advanced modelling features: (i) A model of position-dependent labelling effects that is quantitatively derived from experiment. (ii) Multi-state thermodynamic hybridization models of probe binding behaviour, including potential cross-hybridization reactions. (iii) A fast calibrated sequence-similarity-based heuristic for cross-hybridization prediction supporting large-scale designs. (iv) A novel compound score formulation for the integrated assessment of multiple probe design objectives. In contrast to a greedy search for probes meeting parameter thresholds, this approach permits an optimization at the probe set level and facilitates the selection of highly specific probe candidates while maintaining probe set uniformity. (v) Lastly, a flexible target grouping structure allows easy adaptation of the pipeline to a variety of microarray application scenarios. The algorithm and features are discussed and demonstrated on actual design runs. Source code is available on request.