Development of a 37 k high-density oligonucleotide microarray: a new tool for functional genome research in rainbow trout

A rainbow trout high-density oligonucleotide microarray was constructed using all tentative consensus (TC) sequences that are publicly available from all international rainbow trout Oncorhynchus mykiss genomic research projects through the Rainbow Trout Gene Index database. The new array contains 60-mer oligonucleotide probes representing 37 394 unique TC sequences and 1417 control spots. The array (4 × 44 format) was manufactured according to the design by Agilent Technologies using the inkjet-based SurePrint technology (design number 016320). The performance of the new microarray platform was evaluated by analysing gene expression associated with rainbow trout, vitellogenesis-induced muscle atrophy. This microarray will open new avenues of research that will aid in the development of novel strategies for genetic improvement for economically important traits benefiting the salmonid aquaculture industries.     

Using high-density exon arrays to profile gene expression in closely related species

Global comparisons of gene expression profiles between species provide significant insight into gene regulation, evolutionary processes and disease mechanisms. In this work, we describe a flexible and intuitive approach for global expression profiling of closely related species, using high-density exon arrays designed for a single reference genome. The high-density probe coverage of exon arrays allows us to select identical sets of perfect-match probes to measure expression levels of orthologous genes. This eliminates a serious confounding factor in probe affinity effects of species-specific microarray probes, and enables direct comparisons of estimated expression indexes across species. Using a newly designed Affymetrix exon array, with eight probes per exon for approximately 315 000 exons in the human genome, we conducted expression profiling in corresponding tissues from humans, chimpanzees and rhesus macaques. Quantitative real-time PCR analysis of differentially expressed candidate genes is highly concordant with microarray data, yielding a validation rate of 21/22 for human versus chimpanzee differences, and 11/11 for human versus rhesus differences. This method has the potential to greatly facilitate biomedical and evolutionary studies of gene expression in nonhuman primates and can be easily extended to expression array design and comparative analysis of other animals and plants.     

Improved statistical analysis of budding yeast TAG microarrays revealed by defined spike-in pools

Saccharomyces cerevisiae knockout collection TAG microarrays are an emergent platform for rapid, genome-wide functional characterization of yeast genes. TAG arrays report abundance of unique oligonucleotide ‘TAG’ sequences incorporated into each deletion mutation of the yeast knockout collection, allowing measurement of relative strain representation across experimental conditions for all knockout mutants simultaneously. One application of TAG arrays is to perform genome-wide synthetic lethality screens, known as synthetic lethality analyzed by microarray (SLAM). We designed a fully defined spike-in pool to resemble typical SLAM experiments and performed TAG microarray hybridizations. We describe a method for analyzing two-color array data to efficiently measure the differential knockout strain representation across two experimental conditions, and use the spike-in pool to show that the sensitivity and specificity of this method exceed typical current approaches.     

Efficient oligonucleotide probe selection for pan-genomic tiling arrays

Background  

Array comparative genomic hybridization is a fast and cost-effective method for detecting, genotyping, and comparing the genomic sequence of unknown bacterial isolates. This method, as with all microarray applications, requires adequate coverage of probes targeting the regions of interest. An unbiased tiling of probes across the entire length of the genome is the most flexible design approach. However, such a whole-genome tiling requires that the genome sequence is known in advance. For the accurate analysis of uncharacterized bacteria, an array must query a fully representative set of sequences from the species' pan-genome. Prior microarrays have included only a single strain per array or the conserved sequences of gene families. These arrays omit potentially important genes and sequence variants from the pan-genome.     

Molecular cloning and characterization of cDNAs coding for apopolysialoglycoproteins in cherry salmon (Oncorhynchus masou) eggs

Recently we have cloned the cDNAs and genomic DNAs for apopolysialoglycoproteins (apoPSGPs) of Salmo gairdneri (rainbow trout) [Sorimachi, H., Emori, Y., Kawasaki, H., Kitajima, K., Inoue, S., Suzuki, K., & Inoue, Y. (1988) J. Biol. Chem. 262, 17678-17684], and the sequence analyses have indicated that the mRNAs for apoPSGPs vary in length and contain different numbers of identical 39-bp repeating units encoding the tridecapeptide (Asp-Asp-Ala-Thr-Ser-Glu-Ala-Ala-Thr-Gly-Pro-Ser-Gly) as well as highly conserved sequences encoding pre-, pro-, and telo-peptide regions. In this study we isolated cDNA clones for yamame (cherry salmon, river resident form; Oncorhynchus masou ishikawai) apoPSGP using a genomic DNA fragment for rainbow trout apoPSGP as a probe. The nucleotide sequence analyses revealed that the structures of mRNAs for yamame apoPSGP including the noncoding regions are essentially identical to those for rainbow trout, showing 90% sequence identity. Within the repeating region, 4 bp out of the 39 were replaced, producing a different tridecapeptide, Asp-Asp-Ala-Thr-Ser-Glu-Ala-Ala-Thr-Gly-Pro-Ser-Ser. This tridecapeptide is unique to yamame and common among all cDNAs obtained from yamame. Genomic Southern blot analysis showed that the yamame apoPSGP genes constituted a multiple gene family with a similar gene organization to that of rainbow trout. Oligodeoxynucleotide probes (18 bases) synthesized based on specific sequences for the yamame repeating unit hybridized only to the yamame DNA and not to the rainbow trout DNA, and vice versa.(ABSTRACT TRUNCATED AT 250 WORDS)     

Optimized design and assessment of whole genome tiling arrays

MOTIVATION: Recent advances in microarray technologies have made it feasible to interrogate whole genomes with tiling arrays and this technique is rapidly becoming one of the most important high-throughput functional genomics assays. For large mammalian genomes, analyzing oligonucleotide tiling array data is complicated by the presence of non-unique sequences on the array, which increases the overall noise in the data and may lead to false positive results due to cross-hybridization. The ability to create custom microarrays using maskless array synthesis has led us to consider ways to optimize array design characteristics for improving data quality and analysis. We have identified a number of design parameters to be optimized including uniqueness of the probe sequences within the whole genome, melting temperature and self-hybridization potential. RESULTS: We introduce the uniqueness score, U, a novel quality measure for oligonucleotide probes and present a method to quickly compute it. We show that U is equivalent to the number of shortest unique substrings in the probe and describe an efficient greedy algorithm to design mammalian whole genome tiling arrays using probes that maximize U. Using the mouse genome, we demonstrate how several optimizations influence the tiling array design characteristics. With a sensible set of parameters, our designs cover 78% of the mouse genome including many regions previously considered 'untilable' due to the presence of repetitive sequence. Finally, we compare our whole genome tiling array designs with commercially available designs. AVAILABILITY: Source code is available under an open source license from http://www.ebi.ac.uk/~graef/arraydesign/.     

Expression profiling by oligonucleotide microarrays spotted on coated polymer slides

We have developed a ready-to-spot polymer microarray slide, which is coated with a uniform layer of reactive electrophilic groups using anthraquinone-mediated photo-coupling chemistry. The slide coating reduces the hydrophobicity of the native polymer significantly, thereby enabling robust and efficient one-step coupling of spotted 5' amino-linked oligonucleotides onto the polymer slide. The utility of the coated polymer slide in gene expression profiling was assessed by fabrication of spotted oligonucleotide microarrays using a collection of 5' amino-linked 70-mer oligonucleotide probes representing 96 yeast genes from Operon. Two-colour hybridizations with labelled cDNA target pools derived from standard grown and heat-shocked wild type yeast cells could reproducibly measure heat shock induced expression of seven different heat shock protein (HSP) genes. Moreover, the observed fold changes were comparable to those reported previously using spotted cDNA arrays and high-density 25-mer oligonucleotide arrays from Affymetrix. The low hybridization signals obtained from the DeltaSSA4 mutant cDNA target, together with the high signal detected in two-colour hybridizations with heat-shocked wild type yeast relative to the DeltaSSA4 mutant strain implies that unspecific binding of cDNA target to the SSA4-specific 70-mer oligonucleotide probes is negligible. Combined, our results indicate that the coated polymer microarray slide represents a robust and cost-effective array platform for pre-spotted oligonucleotide arrays.     

Non-sex specific genes associated with the secondary mitotic period of primordial germ cell proliferation in the gonads of embryonic rainbow trout (Oncorhynchus mykiss)

The purposes of this study were to quantify the secondary proliferation of primordial germ cells (PGCs) in both sexes of rainbow trout, determine if a sex difference in the timing of PGC proliferation and eventual pre‐meiotic number exists, and use microarray data collected during this period to identify genes that are associated with PGC mitosis. The experiments used vasa‐green fluorescent protein (vasa‐GFP) transgenic rainbow trout of known genetic sex that allowed for the identification and collection of PGCs in vivo. An increase was observed in the number of PGCs counted in the gonads of both female and male embryonic vasa‐GFP rainbow trout, from 300 to 700° days (water temperature in °C × days post‐fertilization). For both sexes, a statistically significant (P < 0.05) increase in the PGC number was first noted at either 350 or 400° days of development. By 700° days, a 20–50‐fold increase in germ cell number was apparent. No sex‐specific differences in the timing of PGC proliferation or number were notable in any of the families until 700° days. In conjunction, a custom microarray based on cDNA libraries from embryonic rainbow trout gonads was used to identify genes involved in PGC mitosis. Five genes were discovered: guanine nucleotide binding protein, integral membrane protein 2B, transmembrane protein 47, C‐src tyrosine‐protein kinase, and the decorin precursor protein. All the genes identified have not been previously associated with germ cell mitosis, but are known to be involved with the cell plasma membrane and/or cell signaling pathways. Mol. Reprod. Dev. 78:181–187, 2011. © 2011 Wiley‐Liss, Inc.     

70-Kilodalton heat shock polypeptides from rainbow trout: characterization of cDNA sequences.

RTG-2 cells, a line of fibroblasts from rainbow trout (Salmo gairdnerii), are induced to synthesize a distinct set of heat-shock polypeptides after exposure to elevated temperature or to low concentrations of sodium arsenite. We isolated and characterized two cDNA sequences, THS70.7 and THS70.14, encoding partial information for two distinct species of 70-kilodalton heat shock polypeptide (hsp70) from these cells. These sequences are identical at 73.3% of the nucleotide positions in their regions of overlap, and their degree of sequence conservation at the polypeptide level is 88.1%. The two derived trout hsp70 polypeptide sequences show extensive homology with derived amino acid sequences for hsp70 polypeptides from Drosophila melanogaster and Saccharomyces cerevisiae. Northern blot analysis of RNA from arsenite-induced RTG-2 cells, with the trout hsp70 cDNAs as probes, revealed the presence of three hsp70 mRNA species. Southern blot analysis of trout testis DNA cleaved with various restriction endonucleases revealed a small number of bands hybridizing to the hsp70 cDNAs, suggesting the existence of a small family of hsp70 genes in this species. Finally, trout hsp70 cDNA sequences cross-hybridized with restriction fragments in genomic DNA from HeLa cells, bovine liver, Caenorhabditis elegans, and D. melanogaster.     

Error-pooling-based statistical methods for identifying novel temporal replication profiles of human chromosomes observed by DNA tiling arrays

Statistical analysis on tiling array data is extremely challenging due to the astronomically large number of sequence probes, high noise levels of individual probes and limited number of replicates in these data. To overcome these difficulties, we first developed statistical error estimation and weighted ANOVA modeling approaches to high-density tiling array data, especially the former based on an advanced error-pooling method to accurately obtain heterogeneous technical error of small-sample tiling array data. Based on these approaches, we analyzed the high-density tiling array data of the temporal replication patterns during cell-cycle S phase of synchronized HeLa cells on human chromosomes 21 and 22. We found many novel temporal replication patterns, identifying about 26% of over 1 million tiling array sequence probes with significant differential replication during the four 2-h time periods of S phase. Among these differentially replicated probes, 126941 sequence probes were matched to 417 known genes. The majority of these genes were found to be replicated within one or two consecutive time periods, while the others were replicated at two non-consecutive time periods. Also, coding regions found to be more differentially replicated in particular time periods than noncoding regions in the gene-poor chromosome 21 (25% differentially replicated among genic probes versus 18.6% among intergenic probes), while such a phenomenon was less prominent in gene-rich chromosome 22. A rigorous statistical testing for local proximity of differentially replicated genic and intergenic probes was performed to identify significant stretches of differentially replicated sequence regions. From this analysis, we found that adjacent genes were frequently replicated at different time periods, potentially implying the existence of quite dense replication origins. Evaluating the conditional probability significance of identified gene ontology terms on chromosomes 21 and 22, we detected some over-represented molecular functions and biological processes among these differentially replicated genes, such as the ones relevant to hydrolase, transferase and receptor-binding activities. Some of these results were confirmed showing >70% consistency with cDNA microarray data that were independently generated in parallel with the tiling arrays. Thus, our improved analysis approaches specifically designed for high-density tiling array data enabled us to reliably and sensitively identify many novel temporal replication patterns on human chromosomes.     

A putative beta2-adrenoceptor from the rainbow trout (Oncorhynuchus mykiss). Molecular characterization and pharmacology.

Extensive molecular characterization of mammalian beta-adrenoceptors has revealed complex modes of regulation and interaction. Relatively little attention, however, has focused on adrenoceptors from early branching vertebrates such as fish. Using an RT-PCR approach we have cloned a rainbow trout beta2-adrenoceptor gene that codes for a 409-amino-acid protein with the same seven transmembrane domain structure as its mammalian counterparts. This rainbow trout beta2-adrenoceptor shares a high degree of amino-acid sequence conservation with other vertebrate beta2-adrenoceptors. The conclusion that this sequence is a rainbow trout beta2-adrenoceptor is further supported by phylogenetic analysis of vertebrate beta-adrenoceptor sequences and competitive pharmacological binding data. RNase protection assays demonstrate that the rainbow trout beta2-adrenoceptor gene is highly expressed in the liver and red and white muscle, with lower levels of expression in the gills, heart, kidney and spleen of the rainbow trout. The lack of regulatory phosphorylation sites within the G-protein-binding domain of the rainbow trout beta2-adrenoceptor sequence suggests that the in vivo control of trout beta2-adrenoceptor signaling differs substantially from that of mammals.     

Array-MAPH: a methodology for the detection of locus copy-number changes in complex genomes

High-throughput genome-wide screening methods to detect subtle genomic imbalances are extremely important for diagnostic genetics and genomics. Here, we provide a detailed protocol for a microarray-based technique, applying the principle of multiplex amplifiable probe hybridization (MAPH). Methodology and software have been developed for designing unique PCR-amplifiable sequences (400-600 bp) covering any genomic region of interest. These sequences are amplified, cloned and spotted onto arrays (targets). A mixture of the same sequences (probes) is hybridized to genomic DNA immobilized on a membrane. Bound probes are recovered and quantitatively amplified by PCR, labeled and hybridized to the array. The procedure can be completed in 4-5 working days, excluding microarray preparation. Unlike array-comparative genomic hybridization (array-CGH), test DNA of specifically reduced complexity is hybridized to an array of identical small amplifiable target sequences, resulting in increased hybridization specificity and higher potential for increasing resolution. Array-MAPH can be used for detection of small-scale copy-number changes in complex genomes, leading to genotype-phenotype correlations and the discovery of new genes.