50-mer Probes Research Articles

Improved LINE-1 Detection through Pattern Matching by Increasing Probe Length.

Long Interspersed Element-1 (LINE-1 or L1) is an autonomous transposable element that accounts for 17% of the human genome. Strong correlations between abnormal L1 expression and diseases, particularly cancer, have been documented by numerous studies. L1PD (LINE-1 Pattern Detection) had been previously created to detect L1s by using a fixed pre-determined set of 50-mer probes and a pattern-matching algorithm. L1PD uses a novel seed-and-pattern-match strategy as opposed to the well-known seed-and-extend strategy employed by other tools. This study discusses an improved version of L1PD that shows how increasing the size of the k-mer probes from 50 to 75 or to 100 yields better results, as evidenced by experiments showing higher precision and recall when compared to the 50-mers. The probe-generation process was updated and the corresponding software is now shared so that users may generate probes for other reference genomes (with certain limitations). Additionally, L1PD was applied to other non-human genomes, such as dogs, horses, and cows, to further validate the pattern-matching strategy. The improved version of L1PD proves to be an efficient and promising approach for L1 detection.

Resolving the polymorphism-in-probe problem is critical for correct interpretation of expression QTL studies

Polymorphisms in the target mRNA sequence can greatly affect the binding affinity of microarray probe sequences, leading to false-positive and false-negative expression quantitative trait locus (QTL) signals with any other polymorphisms in linkage disequilibrium. We provide the most complete solution to this problem, by using the latest genome and exome sequence reference data to identify almost all common polymorphisms (frequency >1% in Europeans) in probe sequences for two commonly used microarray panels (the gene-based Illumina Human HT12 array, which uses 50-mer probes, and exon-based Affymetrix Human Exon 1.0 ST array, which uses 25-mer probes). We demonstrate the impact of this problem using cerebellum and frontal cortex tissues from 438 neuropathologically normal individuals. We find that although only a small proportion of the probes contain polymorphisms, they account for a large proportion of apparent expression QTL signals, and therefore result in many false signals being declared as real. We find that the polymorphism-in-probe problem is insufficiently controlled by previous protocols, and illustrate this using some notable false-positive and false-negative examples in MAPT and PRICKLE1 that can be found in many eQTL databases. We recommend that both new and existing eQTL data sets should be carefully checked in order to adequately address this issue.

Legionella pneumophila 6S RNA optimizes intracellular multiplication

Legionella pneumophila is a Gram-negative opportunistic human pathogen that infects and multiplies in a broad range of phagocytic protozoan and mammalian phagocytes. Based on the observation that small regulatory RNAs (sRNAs) play an important role in controlling virulence-related genes in several pathogenic bacteria, we attempted to identify sRNAs expressed by L. pneumophila. We used computational prediction followed by experimental verification to identify and characterize sRNAs encoded in the L. pneumophila genome. A 50-mer probe microarray was constructed to test the expression of predicted sRNAs in bacteria grown under a variety of conditions. This strategy successfully identified 22 expressed RNAs, out of which 6 were confirmed by northern blot and RACE. One of the identified sRNAs is highly expressed in postexponential phase, and computational prediction of its secondary structure reveals a striking similarity to the structure of 6S RNA, a widely distributed prokaryotic sRNA, known to regulate the activity of sigma(70)-containing RNA polymerase. A 70-mer probe microarray was used to identify genes affected by L. pneumophila 6S RNA in stationary phase. The 6S RNA positively regulates expression of genes encoding type IVB secretion system effectors, stress response genes such as groES and recA, as well as many genes involved in acquisition of nutrients and genes with unknown or hypothetical functions. Deletion of 6S RNA significantly reduced L. pneumophila intracellular multiplication in both protist and mammalian host cells, but had no detectable effect on growth in rich media.

Microarray analysis of Paramecium bursaria chlorella virus 1 transcription.

Paramecium bursaria chlorella virus 1 (PBCV-1), a member of the family Phycodnaviridae, is a large double-stranded DNA, plaque-forming virus that infects the unicellular green alga Chlorella sp. strain NC64A. The 330-kb PBCV-1 genome is predicted to encode 365 proteins and 11 tRNAs. To monitor global transcription during PBCV-1 replication, a microarray containing 50-mer probes to the PBCV-1 365 protein-encoding genes (CDSs) was constructed. Competitive hybridization experiments were conducted by using cDNAs from poly(A)-containing RNAs obtained from cells at seven time points after virus infection. The results led to the following conclusions: (i) the PBCV-1 replication cycle is temporally programmed and regulated; (ii) 360 (99%) of the arrayed PBCV-1 CDSs were expressed at some time in the virus life cycle in the laboratory; (iii) 227 (62%) of the CDSs were expressed before virus DNA synthesis begins; (iv) these 227 CDSs were grouped into two classes: 127 transcripts disappeared prior to initiation of virus DNA synthesis (considered early), and 100 transcripts were still detected after virus DNA synthesis begins (considered early/late); (v) 133 (36%) of the CDSs were expressed after virus DNA synthesis begins (considered late); and (vi) expression of most late CDSs is inhibited by adding the DNA replication inhibitor, aphidicolin, prior to virus infection. This study provides the first comprehensive evaluation of virus gene expression during the PBCV-1 life cycle.

De Novo Identification of Single Nucleotide Mutations in Caenorhabditis elegans Using Array Comparative Genomic Hybridization

Array comparative genomic hybridization (aCGH) has been used primarily to detect copy-number variants between two genomes. Here we report using aCGH to detect single nucleotide mutations on oligonucleotide microarrays with overlapping 50-mer probes. This technique represents a powerful method for rapidly detecting novel homozygous single nucleotide mutations in any organism with a sequenced reference genome.

Rapid High Resolution Single Nucleotide Polymorphism–Comparative Genome Hybridization Mapping in Caenorhabditis elegans

We have developed a significantly improved and simplified method for high-resolution mapping of phenotypic traits in Caenorhabditis elegans using a combination of single nucleotide polymorphisms (SNPs) and oligo array comparative genome hybridization (array CGH). We designed a custom oligonucleotide array using a subset of confirmed SNPs between the canonical wild-type Bristol strain N2 and the Hawaiian isolate CB4856, populated with densely overlapping 50-mer probes corresponding to both N2 and CB4856 SNP sequences. Using this method a mutation can be mapped to a resolution of approximately 200 kb in a single genetic cross. Six mutations representing each of the C. elegans chromosomes were detected unambiguously and at high resolution using genomic DNA from populations derived from as few as 100 homozygous mutant segregants of mutant N2/CB4856 heterozygotes. Our method completely dispenses with the PCR, restriction digest, and gel analysis of standard SNP mapping and should be easy to extend to any organism with interbreeding strains. This method will be particularly powerful when applied to difficult or hard-to-map low-penetrance phenotypes. It should also be possible to map polygenic traits using this method.

Targeted next-generation sequencing by specific capture of multiple genomic loci using low-volume microfluidic DNA arrays

We report a flexible method for selective capture of sequence fragments from complex, eukaryotic genome libraries for next-generation sequencing based on hybridization to DNA microarrays. Using microfluidic array architecture and integrated hardware, the process is amenable to complete automation and does not introduce amplification steps into the standard library preparation workflow, thereby avoiding bias of sequence distribution and fragment lengths. We captured a discontiguous human genomic target region of 185 kb using a tiling design with 50mer probes. Analysis by high-throughput sequencing using an Illumina/Solexa 1G Genome Analyzer revealed 2150-fold enrichment with mean per base coverage between 4.6 and 107.5-fold for the individual target regions. This method represents a flexible and cost-effective approach for large-scale resequencing of complex genomes.

On-Chip Hybridization Kinetics for Optimization of Gene Expression Experiments

DNA microarray technology is a powerful tool for getting an overview of gene expression in biological samples. Although the successful use of microarray-based expression analysis was demonstrated in a number of applications, the main problem with this approach is the fact that expression levels deduced from hybridization experiments do not necessarily correlate with RNA concentrations. Moreover oligonucleotide probes corresponding to the same gene can give different hybridization signals. Apart from cross-hybridizations and differential splicing, this could be due to secondary structures of probes or targets. In addition, for low-copy genes, hybridization equilibrium may be reached after hybridization times much longer than the one commonly used (overnight, i.e., 15 h). Thus, hybridization signals could depend on kinetic properties of the probe, which may vary between different oligonucleotide probes immobilized on the same microarray. To validate this hypothesis, on-chip hybridization kinetics and duplex thermostability analysis were performed using oligonucleotide microarrays containing 50-mer probes corresponding to 10 mouse genes. We demonstrate that differences in hybridization kinetics between the probes exist and can influence the interpretation of expression data. In addition, we show that using on-chip hybridization kinetics, quantification of targets is feasible using calibration curves.

Efficient high-resolution deletion discovery inCaenorhabditis elegansby array comparative genomic hybridization

We have developed array Comparative Genomic Hybridization for Caenorhabditis elegans as a means of screening for novel induced deletions in this organism. We designed three microarrays consisting of overlapping 50-mer probes to annotated exons and micro-RNAs, the first with probes to chromosomes X and II, the second with probes to chromosome II alone, and a third to the entire genome. These arrays were used to reliably detect both a large (50 kb) multigene deletion and a small (1 kb) single-gene deletion in homozygous and heterozygous samples. In one case, a deletion breakpoint was resolved to fewer than 50 bp. In an experiment designed to identify new mutations we used the X:II and II arrays to detect deletions associated with lethal mutants on chromosome II. One is an 8-kb deletion targeting the ast-1 gene on chromosome II and another is a 141-bp deletion in the gene C06A8.1. Others span large sections of the chromosome, up to >750 kb. As a further application of array Comparative Genomic Hybridization in C. elegans we used the whole-genome array to detect the extensive natural gene content variation (almost 2%) between the N2 Bristol strain and the strain CB4856, a strain isolated in Hawaii and JU258, a strain isolated in Madeira.

Utility of RT-PCR-based Dot-blot Hybridization for Detecting and Genotyping Echoviruses

We attempted to detect and identify virus types quickly by improving an RT-PCR-based dot-blot hybridization test for echoviruses, important human pathogens mainly causing aseptic meningitis. This test was applied to reference viruses of seven echovirus serotypes prevalent in Korea (E6, 7, 9, 11, 13, 25, and 30) and seventy isolates of echovirus isolated in Korea between 2002 and 2004. The primers for target DNA and hybridization probes (25mer, 50mer, and 70mer) were designed within the VP1 region of the echovirus. In RT-PCR, a nonradioactive digoxigenin-DNA labeling mix was added instead of dNTP to initiate PCR. The PCR product was then hybridized against 25mer, 50mer, and 70mer probe DNA spotted on nylon membranes and the reaction was observed. To investigate the optimal conditions for hybridization, various concentrations of target DNA (0.1, 1, 10, and 100 ng/㎕), size of probe DNA (25mer, 50mer, and 70mer), concentrations of probe DNA (10~50 pM), and reaction time were included. In the test zone, the optimal condition in terms of time and cost was a reaction time of 1 h with 10 ng/㎕ target DNA concentration and 10 pM of a 50mer probe. We found 100% diagnosis of the serotypes for seven reference echoviruses and 90% (63/70) sensitivity for clinical isolates. Also, tests with this probe for reactivity with seven reference echoviruses by using DNA chips showed that diagnostic identification was possible without other serotype cross-reactivity. Therefore, efficiency analysis of probe and target DNA on clinical specimens by using dot-blot analysis indicated that this system can be applied to the prestages of the DNA chip and that the dot blot analysis itself can be used in applications to develop a tool for diagnosing specific viral serotypes.

Diversity of Biogeochemical Cycling Genes from Puget Sound Sediments using DNA Microarrays

Oligonucleotide-based microarray permits the simultaneous analysis of thousands of genes on a single chip, so that a better picture of the interactions among thousands of genes can be investigated at the same time. Our oligo microchips contained 763 50-mer probes that scan the region of different functional genes encoding amoA, pmoA, nirS, nirK, nifH, and dsrAB. These genes code for key enzymes in the ecosystem processes of nitrification, methane oxidation, denitrification, nitrogen fixation and sulfur reduction, respectively. We used these oligochips to characterize the distribution of the above genes from Puget Sound sediments at different depths. The composition and distribution of genes from shallower sediments (depths 00.5 cm, 2.0-2.5 cm, 5.0-5.5 cm, and 25.0-25.5 cm) were highly similar but were different from those collected at deeper depths (depths 50-50.5 cm and 84.0-84.5 cm). The deeper sediments present a different community structure with a markedly lower diversity than the shallower depths. Analysis of positive hybridization signals also revealed presence of genes common to all samples. The majority of these genes were similar to those retrieved from various environments (i.e. soils, groundwater, river water, strotomites, marine sediments, and estuarine sediments). Parallel coordinate display showed that the most dominant functional guilds are those that are involved in nitrogen cycling. Our results also indicated that this technology has potential as a tool in revealing a comprehensive "snapshot" of the functional gene composition in marine sediments, although more work is needed to understand the biological meaning of each detectable hybridization signal.

Identification of PML-RARa Target Genes Using Microarray and ChIP-on-Chip Analysis.

PML-RARa is very toxic to early myeloid cell lines that express neutrophil elastase. As a consequence, myeloid cell lines that stably express PML-RARa must undergo significant adaptations to avoid toxicity (eg. loss of expression of neutrophil elastase (Lane and Ley, MCB, 2005). To identify the direct target genes of PML-RARa, we decided to bypass the issue of adaptation by profiling naive U937 cells that transiently expressed PML-RARa. To define “immediate-early” expression changes caused by PML-RARa, we harvested U937 cells electroporated with expression vectors containing eGFP-PML-RARa (known to possess all the measurable functional properties of PML-RARa) vs. eGFP alone, at the earliest time point when adequate numbers of eGFP positive cells could be collected via high speed cell sorting (6 hours). RNA obtained from these cells was subjected to linear amplification and hybridized with Affymetrix U133 Plus 2.0 arrays. Probesets that demonstrated the largest differences in expression between cells transfected with eGFP vs. eGFP-PML-RARa were PML and RARa themselves (15-fold increase for PML, 9-fold for RARa), providing an important validation for the experimental system. Differentially expressed genes were selected by three filters for further study: (i) ANOVA statistic (P < 0.05); (ii) “present” call and >150 average expression units in at least one sample; and (iii) Reproducible (n=2) >2-fold differences between eGFP-PML-RARa vs. eGFP vector transfected cells. Only 280 probesets (out of 54,613 total) met all of these criteria. Among them, 37 genes were involved in signal transduction, 15 in cell proliferation and differentiation, 7 in apoptosis, 14 in cell cycle, 15 in transport, 27 in protein metabolism and modification, 13 in oncogenesis, and 13 in immunity and defense (Panther classification system). 15/15 genes selected from a variety of categories were further validated with Q-PCR in duplicate using replicate samples, and all differences were statistically significant (p<0.05). Q-PCR validated genes that were up-regulated at 6 hours by eGFP-PML-RARa included CD36 (196% increase over eGFP vector control), LATS2 (204%), Sox4 (395%), COPEB (325%), HoxA5 (294%), IFITM3 (139%), AIF1 (197%), and PLD1 (315%). Q-PCR validated genes that were down-regulated by eGFP-PML-RARa at 6 hours included BTG2 (60% the level of eGFP vector control), KLF7(52%), CDKN1A (p21)(39%), Pim1 (58%), SLC23A2 (36%), HoxA13 (51%), and PDCD1LG1 (42%). Among these genes, Sox4, Pim1, BTG2, and PDCD1LG1 have previously been implicated in cancer pathogenesis by retroviral insertional mutagenesis (http://RTCGD.ncifcrf.gov ). To identify promoters that were directly in contact with the eGFP-PML-RARa fusion protein, we performed “ChIP-on-chip” analysis using chromatin immunoprecipitation coupled with a NimbleGen promoter tiling array containing 15 50-mer probes (~100 bp spacing) between positions −1500 and +500 of the transcription initiation sites of a curated subset of 24,275 known human genes (HG17, Build 35). By subtracting non-specific hybridization signals from the eGFP control chip, we identified and validated several promoter regions as direct binding targets of eGFP-PML-RARa, including HoxA13 and CDKN1A (p21). This experimental system provides a new approach for the identification of the direct transcriptional targets of PML-RARa, which may lead to new insights into APL pathogenesis and treatment.

Empirical Establishment of Oligonucleotide Probe Design Criteria

Criteria for the design of gene-specific and group-specific oligonucleotide probes were established experimentally via an oligonucleotide array that contained perfect match (PM) and mismatch probes (50-mers and 70-mers) based upon four genes. The effects of probe-target identity, continuous stretch, mismatch position, and hybridization free energy on specificity were tested. Little hybridization was observed at a probe-target identity of < or =85% for both 50-mer and 70-mer probes. PM signal intensities (33 to 48%) were detected at a probe-target identity of 94% for 50-mer oligonucleotides and 43 to 55% for 70-mer probes at a probe-target identity of 96%. When the effects of sequence identity and continuous stretch were considered independently, a stretch probe (>15 bases) contributed an additional 9% of the PM signal intensity compared to a nonstretch probe (< or =15 bases) at the same identity level. Cross-hybridization increased as the length of continuous stretch increased. A 35-base stretch for 50-mer probes or a 50-base stretch for 70-mer probes had approximately 55% of the PM signal. Little cross-hybridization was observed for probes with a minimal binding free energy greater than -30 kcal/mol for 50-mer probes or -40 kcal/mol for 70-mer probes. Based on the experimental results, a set of criteria are suggested for the design of gene-specific and group-specific oligonucleotide probes, and the experimentally established criteria should provide valuable information for new software and algorithms for microarray-based studies.

Fabrication of DNA microarrays on nanoengineered polymeric ultrathin film prepared by self-assembly of polyelectrolyte multilayers.

Microarray-based technology is in need of flexible and cost-effective chemistry for fabrication of oligonucleotide microarrays. We have developed a novel method for the fabrication of oligonucleotide microarrays with unmodified oligonucleotide probes on nanoengineered three-dimensional thin films that are deposited on glass slides by consecutive layer-to-layer adsorption of polyelectrolytes. Unmodified oligonucleotide probes were spotted and immobilized on these multilayered polyelectrolyte thin films (PET) by electrostatic adsorption and entrapment on the porous structure of the PET film. The PET provides higher probe binding capacity and thus higher hybridization signal than that of the traditional two-dimensional aminosilane and poly-L-lysine coated slides. Immobilized probe densities of 3.4 x 10(12)/cm2 were observed for microarray spots on PET with unmodified 50-mer oligonucleotide probes, which is comparable to the immobilized probe densities of alkyamine-modified 50-mer probes end-tethered on an aldehyde-functionalized slide. The study of hybridization efficiency showed that 90% of immobilized probes on PET film are accessible to target DNA to form duplex format in hybridization. The DNA microarray fabricated on PET film has wider dynamic range (about 3 orders of magnitude) and lower detection limit (0.5 nM) than the conventional amino- and aldehyde-functionalized slides. Oligonucleotide microarrays fabricated on these PET-coated slides also had consistent spot morphology. In addition, discrimination of single nucleotide polymorphism of 16S rRNA genes was achieved with the PET-based oligonucleotide microarrays. The PET microarrays constructed by our self-assembly process is cost-effective, versatile, and well suited for immobilizing many types of biological active molecules so that a wide variety of microarray formats can be developed.

Monte Carlo simulations of free chains in end-linked polymer networks

The structural properties of end-linked polymer networks prepared in the presence of inert linear chain solvent were investigated with Monte Carlo simulations using the three-dimensional bond fluctuation model on a simple cubic lattice. Networks of 50-mer precursor chains were prepared in a solvent of 50-mer inert linear chains with a series of concentrations and two ratios, r, of cross-link sites to chain ends. The networks were formed under both stoichiometric (r=1) and optimal (r=1.2) conditions for minimizing the network imperfections and soluble material and maximizing the elastic material. A maximum is observed in the fraction of elastic material at small degrees of dilution and is explained in terms of entanglement effects. The conformational behavior of a small concentration of linear 50-mer probe chains trapped in end-linked networks of mesh sizes ranging from 10- to 50-mer was also studied. The radius of gyration of the linear chains was found to decrease with decreasing mesh size of the host network, in agreement with a theoretical scaling relationship; but the magnitude of the effect is small.