In vitro identification and in silico utilization of interspecies sequence similarities using GeneChip® technology
BMC Genomics volume 6, Article number: 62 (2005)
Genomic approaches in large animal models (canine, ovine etc) are challenging due to insufficient genomic information for these species and the lack of availability of corresponding microarray platforms. To address this problem, we speculated that conserved interspecies genetic sequences can be experimentally detected by cross-species hybridization. The Affymetrix platform probe redundancy offers flexibility in selecting individual probes with high sequence similarities between related species for gene expression analysis.
Gene expression profiles of 40 canine samples were generated using the human HG-U133A GeneChip (U133A). Due to interspecies genetic differences, only 14 ± 2% of canine transcripts were detected by U133A probe sets whereas profiling of 40 human samples detected 49 ± 6% of human transcripts. However, when these probe sets were deconstructed into individual probes and examined performance of each probe, we found that 47% of human probes were able to find their targets in canine tissues and generate a detectable hybridization signal. Therefore, we restricted gene expression analysis to these probes and observed the 60% increase in the number of identified canine transcripts. These results were validated by comparison of transcripts identified by our restricted analysis of cross-species hybridization with transcripts identified by hybridization of total lung canine mRNA to new Affymetrix Canine GeneChip®.
The experimental identification and restriction of gene expression analysis to probes with detectable hybridization signal drastically increases transcript detection of canine-human hybridization suggesting the possibility of broad utilization of cross-hybridizations of related species using GeneChip technology.
Genome-wide analyses of multiple organisms in a variety of experimental biological systems are powerful tools employed in biomedical research. However, microarray platforms for numerous species, particularly large mammals, have yet to be fully developed. Large mammals, often preferred species for modeling of various pathophysiological conditions and environmental responses, demonstrate better concordance to humans with regard to toxicological effects compared to small animal models (rodents) [1, 2]. Further, large animal models are indispensable in disease modeling that requires longer live span , in organ transplant research , and in studies where multiple samples are required to be collected . For organ- and system-specific applications, such as the respiratory system, canine and ovine lungs are physiologically and hemodynamically more closely related to human lungs than rodent lungs and allow regional evaluation of air, blood, and exudate distribution. Thus, canine [6, 7] and ovine [8, 9] models are the most commonly utilized in pathophysiological modeling of the respiratory system.
Despite these obvious advantages of the large mammal model, the lack of large mammal genomic sequence data presents significant limitation for genome-wide analysis of gene expression profiles by microarray techniques. To overcome this limitation, cross-species RNA hybridization has been utilized as one potential solution [10–13]. This approach is based on the notion that evolutionarily-conserved genetic sequences between mammals will be efficiently detected during cross-species hybridization. Despite a paucity of literature focused on the use of Affymetrix GeneChip oligonucleotide array for cross-species hybridizations [14, 15], we speculated that this array platform could potentially serve as a powerful tool for cross-species genome analysis. GeneChips carry 11–16 oligonucleotide probes per target gene, increasing the opportunity to potentially identify genes with stretches of high homology to specific probes. To test this hypothesis, we cross-hybridized total RNA obtained from canine lung tissues to the human HG-U133A GeneChip. Analysis of generated expression profiles using Affymetrix MicroArray Suite (MAS 5.0) showed that approximately 14% of canine mRNA hybridizes to the human probe sets and generates hybridization signal similar to that of human mRNA. Analysis of available canine orthologues corresponding to these well performing human probes revealed high sequence similarities between human and canine counterparts. Based on these findings we speculated that even canine genes with only partial sequence similarities to their human counterparts can be identified by MAS 5.0 if the analysis is restricted to probes with high sequence similarity to their canine counterparts. Antipova et al. recently reported that utilization of a single well performing Affymetrix probe pair can be sufficient for detection of its corresponding target , further supporting our hypothesis.
A straight forward approach to test this hypothesis would be to align the target sequences from the human HG_U133A GeneChip to the known canine sequences using Basic Local Alignment Search Tool (BLAST), followed by identification of human probes with high sequence similarity to their canine counterparts. Unfortunately, due to limited canine genomic information, only 3% of human target sequences from the HG_U133A GeneChip have been matched to known canine orthologues, which renders this approach inadequate. However, we speculated that since canine orthologues with high sequence similarities to human probes generate detectable hybridization signals, the reverse should be also true: the detectable hybridization signal should reflect high sequence similarity. Therefore, we examined individual probes instead of complete probe sets. Probes with detectable hybridization signals were considered adequate representatives of their corresponding canine orthologues and were thus employed for final signal intensity calculations. The MAS 5.0 analysis when restricted to probes with detectable hybridization signals, demonstrated a 4-fold increase in canine transcript detection as compared to unrestricted analysis. Our results strongly suggest that identification of probes with detectable cross-species hybridization signal and subsequent restriction of gene expression analysis to these probes can be applied to multiple related species for which gene arrays are not commercially available.
The hybridization pattern of the protein translocation complex beta gene involved in cell maintenance (SEC61B, Gene Ontology #8151) is assumed to be constantly expressed in the majority of tissues and was selected for demonstration of the probe level analysis. This analysis was conducted using signal intensity values generated by hybridization of human and canine SEC61B mRNAs to HG-U133A GeneChip. The probe level analysis deconstructs the entire Affymetrix probe set into its individual components and allows one to evaluate the performance of each individual probe. Detailed output of human SEC61B hybridization to 203133_at probe set (Figures 1A and 2A) shows that signal intensity value of PM (perfect match or total signal) obtained for 11th probe pair is higher than that of MM (mismatch or background signal). Therefore, correction for the background will produce positive signal intensity value, and this probe pair will contribute to the "Present" (detectable transcript) classification by MAS 5.0. However, cross-species hybridization of canine mRNA to the same probe pair does not produce sufficient signal intensity value (Figures 1C and 2B), possibly due to differences in target sequence between canine and human. This will, therefore, shift the overall probe set classification towards "Absent" (non-detectable transcript) classification. At the same time, there are probes that failed to detect both human and canine SEC61B. This could be due to either secondary structure formation (Figure 1B and 1D)  or the presence of sequences in the hybridizing mRNA mixture that possess homology to the 11th MM probe, thus raising the background signal over the specific signal. The pattern of SEC61B hybridization shows that, among 11 probe pairs of the 203133_at probe set, 3 probes (#2, #3, and #5) failed to detect its target in human tissues but the detectable hybridization of the remaining 8 probes prevailed (Figure 2A), resulting in classification of the SEC61B transcript by MAS 5.0 as "Present". In contrast, a reduction in the number of probes that were able to detect canine SEC61B down to only 6 probes (#1, ##6–10, Figure 2B) caused MAS 5.0 to classify this gene as "Absent". As shown in Figure 2C, the failure of the 5 non-hybridizing canine SEC61B probes was mainly attributed to differences in gene sequences between these two species.
This observation held true for most of the evaluated targets, leading us to speculation that among the 11 human probes of a given probe set there will be probes with high sequence similarity to its canine target and that identification and utilization of these probes will increase sensitivity of cross-species hybridization. Moreover, when we evaluated performance of all 9,871,960 HG-U133A probes from 40 cross-species hybridizations, we found that almost half of these probes (4,688,566) generated a detectable hybridization signal (PM-MM>20). Therefore, we speculated that excellent and poorly performing probes should be equally distributed in the probe sets as well, and that masking of the poorly performing half will increase transcript detecting ability of these modified probe sets. The test run of MAS 5.0 analysis using masking of the 5 poorly performing probes of 203133_at probe set converted the initial "Absent" call for canine SEC61B gene to a "Present" call and will formulate the basis for further investigation. Since the canine genome is not yet completed, we were unable to select probes with high canine/human sequence similarity solely based on gene sequence information. Instead, we evaluated the overall performance of each individual probe based on their hybridization signal. We speculated that sequences of human probes with a reasonable cross-species hybridization signal (PM-MM>20) most likely possess a high similarity to their canine targets.
The probe level analysis was performed for 40 expression profiles generated by hybridization of 8 different human tissues and 40 expression profiles generated by hybridization of canine lung tissues to the human HG-U133A GeneChip. The resulting hybridization signal pattern of the 203133_at probe set for both species is presented in Figure 2D and 2E. As expected, differences in the canine SEC61B mRNA sequence region corresponding to the 11th probe (Figure 2C) did not allow canine target detection in 39 out of 40 hybridizations (Figure 2E). These data suggest that once a probe poorly hybridizes to its target in several initial experiments, it will most likely fail to detect its target in the substantial number of consequent hybridizations. To investigate whether this pattern persists throughout whole gene array, we analyzed hybridization patterns of all 246799 probes of the HG-U133A GeneChip. The resulting histogram for the same species hybridization (Figure 3A) demonstrates that a large fraction (~44% cases) of probes performed well during 38 or more hybridizations (95%). This trend was also true for poorly performing probes, where a considerable fraction (~11%) probes failed to detect their target during 38 or more hybridizations. In canine/human hybridization the number of poorly performing probes was higher (24%) while fraction of well performing probes lowered to 38% (Figure 3B). The analysis of the overall net effect of canine/human sequence differences on the hybridization pattern was conducted by correction for poorly performing probes during human samples hybridization. This revealed a 7% decrease in well performing and 13% increase in poorly performing probe fractions (Figure 3C). For detailed evaluation of direct effects of interspecies sequential differences on the cross-species hybridization, we selected probe sets where all probes generated detectable hybridization signal during human sample hybridization. Interestingly, among these probes, the fraction of probes that detected their corresponding canine targets during 38 or more cross-species hybridization was surprisingly high (~72%), with a very low (~3%) fraction of probes that failed to detect their targets in canine tissues (Figure 4). To evaluate correlation between performance of an individual probe and its alignment to the corresponding canine target, the number of mismatches were identified and counted for 20 randomly selected canine genes. This analysis demonstrated that, on average, probes in both excellent and poorly performing fractions had 1 and 2.5 mismatched bases, respectively, with the mismatch rate significantly lower in excellent performing probes (Figure 4, insert). The presence of a mismatch in probes with detectable hybridization signal was not unexpected and correlates with previously reported observation that relatively long (≥ 16 bp) stretch of the perfect match could produce stable hybridization and signal generation .
The criterion for masking of a particular probe was based on the performance of a given probe throughout 40 GeneChips. We found that masking of probes which poorly performed (PM-MM≤20) in 15 or more hybridizations (~40% of all hybridizations) produced the highest number of "Present" calls (Figure 5A) therefore this percentage was selected as a cut off value for masking. We speculated that probes that failed to detect its target in more than 40% of hybridizations lack homology to their canine counterparts at least in part due to sequence divergence. To test this hypothesis we conducted a sequence analysis of randomly selected samples and find that out of the 42 tested probes 27 probes aligned against canine sequence with one or more mismatches and 22 of these probes (approximately 82%) were masked by our procedure.
We also analyzed the effects of our procedure on the probe set size (Figure 5B). The main fractions of probe sets adjusted to human (77% of all probe sets) or canine (81% of all probe sets) hybridizations are represented by 5–10 and 4–8 probe pairs, respectively. Our masking procedure left 2172 (10%) probe-sets unaltered and completely eliminated 7 (0.03%) probe sets based on the human hybridization pattern. Analysis of canine hybridization identified 184 (0.8%) highly homologous canine targets, which hybridize to the entire corresponding human probe sets, and eliminated 11 probe-sets (0.05%) where none of the probes recognizes its target in canine tissues.
We next examined relative location of poorly performing human probes on the canine genes. It has been reported previously that human/canine homology is prone to be higher within the coding region . Our findings are consistent with this report and show that our approach masked more probes towards the 3' untranslated region of target sequences (Figure 5C).
The effect of our masking procedure was also evaluated for transcript abundance detection. The combined (40 GeneChips) output of regular MAS 5.0 analysis was compared with outputs produced using masking files generated for human and canine hybridizations (Figure 6). In human hybridization, our algorithm reduced "Absent" calls by 3 fold and increased "Present" calls by more than 1.5 fold. The application of the masking file generated for canine hybridization had a more modest effect than human-specific masking. The transcript detection during the hybridization of canine mRNA to the human HG-U133A chip was much less sensitive (14%) than that of hybridization of human mRNA (48%). Our canine-specific masking recovered approximately 60% (Figure 6) of transcripts which, otherwise, were called "Absent" by MAS 5.0. The anticipated increase in transcript detection by application of human-specific masking is attributed to the masking of imperfectly designed probes. As shown in Figure 1, the 11th probe of the 203133_at probe set failed to detect its target in both canine and human tissues. Therefore, masking of such a probe will improve transcript detection in both canine and human sample analyses. Similarly, the canine-specific masking of human hybridization also increased the number of detectable transcripts as compared to non-masking results (Figure 6A).
Finally, we compared transcripts detected by our masking procedure with those identified by hybridization of canine lung tissues to the newly developed Affymetrix canine GeneChip. We have identified 996 orthologous transcripts common to both HG-U133A and canine GeneChips and compared their detection in both arrays using "Present"/"Absent" classifications by MAS 5.0 or GCOS 1.1, respectively. This analysis showed that our masking approach increased the fraction of correctly identified transcripts by 46% and decreased the false negative fraction by 53% as compared to unmasked transcript detection (Figure 7).
The widespread use of large mammal models in respiratory research, coupled to the limited availability of corresponding gene microarrays, motivated us to investigate the feasibility of utilizing gene arrays designed for related organisms. In this work, we hybridized canine lung tissues to human Affymetrix HG-U133A GeneChip, and demonstrated that oligonucleotide probes designed for identification of human targets were able to generate detectable hybridization signal for canine tissues (Figure 6B). However, the number of identified transcripts was 3–4 times lower than that of the transcripts which were identified during human sample hybridization (Figure 6A). Our analysis indicated that the low sensitivity of the majority of the HG-U133A probe sets to canine transcripts was due to partial genomic differences between human and canine sequences (Figure 2B and 2C). The abundance of probes in Affymetrix probe sets (11–16 probes per set) presented an opportunity for reduction of probe sets to probes that generate detectable cross-species hybridization signal and therefore, improving transcript detecting ability of such modified probe sets. Our probe level analysis allows evaluation and identification of probes that failed to detect their targets in canine tissues. We demonstrated that this poor performance was not random, but rather linked to either imperfect probe design (Figure 1B) or interspecies differences in genomic sequences (Figure 4). Another indirect linkage of probe/target sequence similarity and probe performance was demonstrated by inclination of poorly performing probes toward the 3' end of a target (Figure 5C), which is more prone to diversity . The masking of this poorly performing probes resulted in an increase of "Present" calls in both species (Figure 5A), where the increase in human transcript detection was attributed mainly to masking of imperfectly designed probes. Interestingly, the increase of "Present" calls was more prominent in canine hybridizations when we masked the most poorly performing probes that were not able to detect their targets in more that 30 hybridizations. However, thereafter the difference in number of "Present" calls between canine and human remains constant (Figure 5A). This comparable increase in the amount of "Present" calls is most likely species-independent and attributed to probe set optimizations, which is reflected by changes in probe set sizes. We demonstrated that masking probes that were not detecting their targets in 40% cases generated maximum of "Present" calls in both species (Figure 5A) and speculated that probe pairs that poorly performed in more than 40% of hybridizations lack homology to their canine counterparts. This was confirmed by alignment analysis of these probes to their canine counterparts (Figure 2), which showed that 82% of masked probes produced one or more mismatches with their canine counterparts. The remaining 18% can be attributed to a sequencing error, a single nucleotide polymorphism, or longer than 16 bp stretches of perfect matches that can produce detectable hybridization signal . We also identified probes that were perfectly aligned to their canine counterparts yet perform poorly. These probes represented approximately 7% of tested probe pairs, a value similar to that previously reported  and can be attributed to the secondary structure or to the presence of yet unknown sequence that recognizes MM probe as a perfect match and the resulting increase in background signal obscure the hybridization signal.
We demonstrated that modification of probe sets by masking poorly performing probes drastically increased detection of canine transcripts by the human HG-U133A GeneChip (Figure 6B). To validate these results, we compared transcripts identified by our masking approach with those identified by hybridization of canine lung tissues RNA to the Affymetrix canine-specific GeneChip. This comparison demonstrated that modification of probe sets prior to MAS 5.0 analysis drastically increases correlations between transcripts identified by same- and cross-species hybridizations (Figure 7). The slight increase (from 3% to 10%) in false positive transcript fraction introduced by our procedure is a common feature of high density oligonucleotide array analyses  and is well compensated by the approximately 10-fold decrease in the false negative transcript fraction from 59% to 6% (Figure 7).
The studies presented here show that the Affymetrix microarray platform can be successfully used for cross-species hybridization of related organisms. A simplified version of this method (probes that failed to generate detectable signal in all hybridizations were masked) has already been successfully applied to small (7 GeneChips) data sets  and produced biologically sensible results. We expect that method presented here will be applicable not only for canine/human hybridization but also for other combinations of cross-species hybridizations. Whenever there is homology between species, there is an opportunity for identification probes with detectable hybridization signal using probe level analysis described here. This method will accelerate genome-wide analysis in experimental models for which species-specific gene arrays are not currently available.
Animal preparation and sample collection
All experimental procedures were approved by the Johns Hopkins University Animal Care and Use Committee. Four mongrel dogs (weight 21.3 ± 1.5 kg) were anesthetized with 25 mk/kg pentobarbital i.v. and instrumented with femoral arterial and venous catheters. Anesthesia was maintained with additional pentobarbital (5 mg/kg iv every hour and when indicated) and muscle relaxation provided by pancuronium (3 mg bolus and 0.5 mg hourly iv). A 39 or 41 French double-lumen endobronchial tube (Mallinkrodt, St. Louis, Missouri) was placed via a tracheostomy and left lung was then mildly injured by repeated lavage with warmed saline. One animal was sacrificed by exsanguination after supplemental pentobarbital (10 mg/kg i.v.) right after lavage, the chest opened, and 10 control tissue samples were taken from 5 corresponding regions in both lungs (apex- dependent and non-dependent, base- non dependent, mid, and dependent). Other 3 animals were ventilated for 5 hours, then sacrificed as described above and 30 samples (10 from each animal) were collected following the control animal pattern. The lung tissue samples were immersed in RNAlater (Ambion, Austin, Texas), snap-frozen and stored at -80°C until processed for RNA isolation. The protocol for cross-species hybridization sample preparation has been described by our group previously [20, 21]. Briefly, total RNA was isolated from individual lung samples described above using the Trizol reagent and Qiagen RNeasy columns. Individual cDNAs were prepared from each RNA isolate using reverse transcriptase (GIBCO-BRL SuperScript). Each cDNA was subsequently used as a template to make biotin-labeled cRNA using an in vitro transcription reaction (Enzo), resulting in a single cRNA for each lung sample. Each cRNA was hybridized with an individual Affymetrix HG-U133A oligonucleotide array, which was subsequently processed and scanned according to the manufacturer's instructions. All arrays were hybridized in the same batch to avoid variability in hybridization conditions. The sample collection and preparation for the same specie hybridization was adapted from the described above protocol and applied for Affymetrix Canine GeneChip hybridization by McVerry et al .
Human tissues and cell lines sampling and RNA isolation
All human samples and analyses posted on the HopGene server were generated with IRB approval (#B0112210102). Pathologic samples of human myocardium 1 from unused donor heart, 2 from patients with ischemic cardiomyopathy, and 2 from patients with nonischemic cardiomyopathy were collected at time of transplantation. Peripheral blood mononuclear cells were collected form 1 healthy person, 2 smokers, and 2 patients with allergic asthma and separated by Ficoll-hypaque density gradient centrifugation. Airway brushing cells were collected in normal saline form 1 healthy person, 3 smokers, and 1 patient with allergic asthma and separated by centrifugation. Bronchoalveolar lavage (BAL) was performed during bronchoscopy of 2 healthy individuals, 2 smokers, and 2 patients with allergic asthma by instilling 100 ml of normal saline, BAL fluid was aspirated into tubes and cells collected by centrifugation. The muscle tissue was collected from 5 individuals with different stages of obesity. All human tissues were snap-frozen in liquid nitrogen upon collection. These frozen samples (~50 mg from each sample) were directly solubilized in chaotropic solubilization buffer using a Brinkman Polytron tissue disruptor. Larger tissue fragments (>100 mg) were pulverized into frozen powder with a mortar and pestle, pre-chilled to liquid nitrogen temperature, and then the frozen powder was solubilized with the Polytron. RNA was purified using Trizol LS (Life Technologies) and an additional RNA purification step was conducted using the RNAeasy purification kit (Qiagen Inc., Valuencia, CA).
Human smooth muscle cell cultures were treated with ApoC-1 for 10 min (n = 2), 24 h (n = 2) or left intact (n = 1); human pulmonary artery endothelial cells were exposed to hyperoxia for 6 h (n = 1), 24 h (n = 1), 48 h (n = 1), or normoxia (n = 2); IB3-1 bronchial epithelial cells were treated with 1 mM PBA for 24 h (n = 1) or left intact (n = 4). Cells were rinsed twice with ice cold Hank's buffer then scraped into 3ml of ice cold TRIzol (Gibco BRL). Total RNAs was purified using Trizol LS (Life Technologies) and an additional RNA clean-up step was conducted using the RNAeasy purification kit (Qiagen Inc., Valuencia, CA). Approximately 10 μg of purified, total RNA was used for analyses.
RNA preparation and hybridization
Purified total RNA was reverse transcribed to first-strand cDNA using a hybrid primer consisting of oligo-dT and T7 RNA polymerase promoter sequences. The single-stranded cDNA was then converted to double-stranded cDNA. Complementary DNA corresponding to 5–10 μg of total RNA was used in a cRNA amplification step using T7 RNA polymerase and two biotinylated nucleotide precursors. The resulting biotinylated cRNA was fragmented to a size of approximately 50 bp. Approximately 20–30 μg of cRNA from each tissue was hybridized to corresponding (HG_U133A or Canine) GeneChips (Affymetrix, Santa Clara, CA). The bound cRNA was visualized by binding of streptavidin/phycoerythrin conjugates to the hybridized GeneChip, followed by laser scanning of bound phycoerythrin. These scan results are posted on HopGene ftp server and are freely accessible at ftp://188.8.131.52/pga/rawdata/public/BMC_Genomics_DG_1
Expression Data Analysis
Data analysis was conducted using Affymetrix MicroArray Suit 5.0 (MAS 5.0) and GeneChip Operating Software (GCOS 1.1) (Affymetrix Inc., Santa Clara, CA). Signal intensity values for hybridization of human samples or heterohybridization of canine samples to HG-U133A (40 GeneChips for each specie), and hybridization of canine lung tissues to Affymetrix canine GeneChips (12 GeneChips) were generated using user-definable parameters set to Affymetrix default values and total chip fluorescence intensities scaled to 150. Each GeneChip was analyzed with and without application of custom generated probe masking files. The probe masking function is imbedded into MAS 5.0 and GCOS 1.1 and can be easy evoked and linked to a list of probes selected for masking. The overall schema of analysis of 40 human HG-U133A GeneChips, which were hybridized with canine mRNA is outlined in Figure 8.
Creating "Perfect match"/"Mismatch" probe level matrix
Each HG-U133A GeneChip contains 22,215 (excluding internal controls) probe sets, each representing one target sequence. The majority of the probe sets (98 %), contain 11 probe pairs with each probe pair consisting of two 25-mer oligonucleotides. One of these oligonucleotides is a "perfect match" (PM) to a target sequence and another is a "mismatch" (MM) with an intentionally mutated middle nucleotide (Figure 1A). Therefore, each probe set contains 11 PM probes and 11 corresponding MM probes for calculation of total and nonspecific hybridization signals, respectively (reader is referred to Affymetrix website for detailed explanations, http://www.affymetrix.com). The MAS 5.0 program analyzes the differences in the amount of hybridized mRNA to each PM (total signal) and MM probes (background signal) across entire probe set, and computes the relative abundance of a transcript and its significance (p value). Transcripts with p < 0.04 are called "Present," with 0.04 ≤ p ≤ 0.06 called "Marginal," and with p > 0.06 called "Absent". The evaluation of the effects of our masking approach on the change in the transcript detection call was based solely on the Affymetrix transcript abundance classification without modification of the original algorithm, which makes employment of our masking files universal among Affymetrix users.
Expression profiles analyzed in this report were generated using total mRNA extracted from different parts of canine lungs exposed to 15 different experimental conditions. Signal intensities files (CEL files) generated by MAS 5.0 and GCOS 1.1 for HG-U133A and Canine GeneChip, respectively were imported into the R-based Bioconductor analytical package http://www.hopkins-genomics.org/resources/resourcesSoftware.html and processed using affy library . Intensity signals for each experimental group were corrected for background and normalized using default mas5 package settings, which closely simulate the MAS 5.0 algorithm. The raw intensities of PM and MM values for each of 246799 probes per chip (excluding internal Affymetrix controls) were extracted from mas5 generated probe level matrix  and converted into the flat text file with 80 by 246799 entries (Figure 8). The R-language scripts for extraction of PM and MM values are available upon request. Further manipulations and analyses of extracted data were conducted using HopGene designed algorithms written in Python 2.2. Statistical parts of the script employed the Python Statistical Module developed by NSG group http://www.nmr.mgh.harvard.edu/Neural_Systems_Group/gary/python.html
Correlation of the human-probe/canine-target similarity and hybridization signal
To evaluate relation between probe/target sequence similarity and hybridization signal all probe pairs were classified as efficient (PM-MM>20, the minimum practical signal value ) or inefficient (PM-MM≤20). This classification was applied to individual probe pairs throughout 40 canine/human hybridizations and 40 human/human hybridizations. Probe level analysis of human/human hybridization was aimed at identification of non-functional due to poor design or inappropriate hybridization conditions probe pairs (Figure 1B). These probes were used for identification of the net effect of canine/human genetic differences, that is probes that were inefficient in both human/human and canine/human hybridizations (Figure 4) were eliminated from further analysis. Target sequences (FASTA format) representing canine net effect probes (Figure 4) were retrieved from NetAffx http://www.affymetrix.com/analysis/index.affx and queried against NCBI non-redundant (nr) database for their corresponding canine orthologues using MegaBLAST http://www.ncbi.nlm.nih.gov/BLAST. Since the MegaBLAST batch query tool limits the number of submitting sequences, we divided this human dataset into 3 batches (~7000 sequences each) and used Canis familiaris term as the search restriction. BLASTN 2.2.8 [Jan-05-2004] analysis employed sequences form GenBank, EMBL, DDBJ, PDB databases, excluding EST, STS, GSS, or phase 0, 1 or 2 HTGS sequences and showed that 662 target sequences on the HG_U133A GeneChip had 1122 canine counterparts with >75% sequence homology. For manual evaluation 20 target sequences were randomly selected from the pool of sequences that contained probes efficient in less that 5% hybridizations (hybridized only to 0–2 chips, Figure 4). Since each target sequence contains 11 probes the total of 220 probes was evaluated and 22 probes that efficiently hybridized only to 0–2 GeneChips (inefficient group) and 31 probes that efficiently hybridized to 38–40 (95%) chips (efficient group) were aligned to their canine counterparts. Alignments of these 53 probes to their canine counterparts were generated using CLUSTALW http://clustalw.genome.ad.jp and nucleotide mismatches were counted and accordingly assigned to the "efficient" or "inefficient" group. Since mismatch distribution was heavily skewed towards small numbers of mismatches (1–2) the correlation between probe performance and target sequence homology (number of identified mismatches) was evaluated using nonparametric Mann-Whitney test with p < 0.05 used as a significance cut-off (Figure 4, insert).
Creating probe pair masking files
Hybridization of individual probe pairs (246799total) to each GeneChip (40 total) was evaluated and inefficient probe pairs that produced hybridization signal below minimum practical signal value  (PM-MM≤20) were identified. The flat text file containing two columns was generated, with the first column containing probe set ID entry and the second column containing order numbers of probe pairs to be masked. The empty template masking file was generated using MAS 5.0 masking tool and subsequently populated with data from the flat text file. The first masking file was built from probe pairs that were inefficient in all 40 hybridizations. The second masking file contained probe pairs that were inefficient in ≥35 hybridizations. The subsequent masking files were generated with 5 chosen as an increment of inefficient hybridizations with the last masking file comprised of probe pairs that were inefficient in the only one hybridization. Resulting 9 masking files (40, 35, 30, 25, 20, 15, 10, 5, 1) were applied during MAS 5.0 analysis of canine and human mRNA hybridizations to the HG-U133A GeneChips. The highest percentages of the "Present" call (Figure 5) were generated around the point at which probes that were inefficient in ≥15 hybridizations (approximately 40%) were masked and was considered an optimum condition for a given dataset. In further discussions the percentages of the described above masking points were conveniently approximated.
Evaluating effects of masking approach
To simplify the evaluation of our masking procedure, only probe-sets containing 11 probe pairs, which comprises ~98% of all HG_U133A probe-sets, were assessed. The remaining 2% of probe-sets served as internal controls (20 probe pairs per set) and probe sets built from 14–16 probe pairs were inherited from the HG_U95AVs2 GeneChip. Our masking probe pair list for 11-member probe-sets contained 156,174 probe pairs (~65% of total). All counting and poorly performing probe selecting protocols were scripted in Python 2.2 http://www.python.org and are available upon request.
Concordance of the masking HG_U133A probes with their poor homology to the corresponding canine counterparts was evaluated using analysis of 5 human/canine alignments including one presented in Figure 2C and four randomly selected alignments from three representative pools of probe sets: probe sets where all probes were detectable, probe sets where none of the probes were detectable and probe sets represented by both detectable and non-detectable probes ftp://184.108.40.206/pga/rawdata/public/BMC_Genomics_DG_1/Examples/.
Affymetrix Canine and HG_U133A GeneChip comparison
Identification of homologous to Affymetrix Canine GeneChip target sequences on human HG_U133A GeneChip was conducted using MegaBLAST. The Affymetrix Canine GeneChip's target sequences in FASTA format were retrieved from NetAffx, split into 12 batches (~2000 sequences each) to overcome MegaBLAST limitations and queried against nr database using Homo sapiens term for filtering. GeneBank accession numbers of human orthologues identified by MegaBLAST for canine chip targets were matched against those on HG-U133A GeneChip using Microsoft Access 2000. Detection calls for common orthologues between 40 cross-hybridized HG-U133A GeneChips and 12 same-species hybridized canine GeneChips were compared. The transcripts that were efficiently detected in ≥ 40% (described above optimum condition) of tested GeneChips were considered as expressed. Then detection of these expressed in canine lung tissues transcripts was compared with their corresponding detection in original and masked cross-hybridization of canine lung tissues to human HG-U133A GeneChip. Transcripts detected by Canine GeneChip but missed by HG-U133A cross-hybridization were "false negative" and those not detected by Canine GeneChip but shown by HG-U133A heterohybridization were "false positive".
Olson H, Betton G, Robinson D, Thomas K, Monro A, Kolaja G, Lilly P, Sanders J, Sipes G, Bracken W, Dorato M, Van Deun K, Smith P, Berger B, Heller A: Concordance of the toxicity of pharmaceuticals in humans and in animals. Regul Toxicol Pharmacol. 2000, 32: 56-67. 10.1006/rtph.2000.1399.
Olson H, Betton G, Stritar J, Robinson D: The predictivity of the toxicity of pharmaceuticals in humans from animal data--an interim assessment. Toxicol Lett. 1998, 102-103: 535-538. 10.1016/S0378-4274(98)00261-6.
Hein WR, Griebel PJ: A road less travelled: large animal models in immunological research. Nat Rev Immunol. 2003, 3: 79-84. 10.1038/nri977.
Pirenne J: Contribution of large animal models to the development of clinical intestinal transplantation. Acta Gastroenterol Belg. 1999, 62: 221-225.
Blache D, Adam CL, Martin GB: The mature male sheep: a model to study the effects of nutrition on the reproductive axis. Reprod Suppl. 2002, 59: 219-233.
Gu MN, Xiao JF, Huang YR, Zhou W, Fu WJ, Chen YM, Piao YJ, Chen PY, Zeng FY: Study of direct lung injury by seawater in canine models. Di Yi Jun Yi Da Xue Xue Bao. 2003, 23: 201-205.
Hubloue I, Biarent D, Abdel Kafi S, Bejjani G, Melot C, Naeije R, Leeman M: Endothelin receptor blockade in canine oleic acid-induced lung injury. Intensive Care Med. 2003, 29: 1003-1006.
Enkhbaatar P, Murakami K, Shimoda K, Mizutani A, McGuire R, Schmalstieg F, Cox R, Hawkins H, Jodoin J, Lee S, Traber L, Herndon D, Traber D: Inhibition of neuronal nitric oxide synthase by 7-nitroindazole attenuates acute lung injury in an ovine model. Am J Physiol Regul Integr Comp Physiol. 2003, 285: R366-72.
Murakami K, Bjertnaes LJ, Schmalstieg FC, McGuire R, Cox RA, Hawkins HK, Herndon DN, Traber LD, Traber DL: A novel animal model of sepsis after acute lung injury in sheep. Crit Care Med. 2002, 30: 2083-2090. 10.1097/00003246-200209000-00022.
Chismar JD, Mondala T, Fox HS, Roberts E, Langford D, Masliah E, Salomon DR, Head SR: Analysis of result variability from high-density oligonucleotide arrays comparing same-species and cross-species hybridizations. Biotechniques. 2002, 33: 516-8, 520, 522 passim.
Moody DE, Zou Z, McIntyre L: Cross-species hybridisation of pig RNA to human nylon microarrays. BMC Genomics. 2002, 3: 27-10.1186/1471-2164-3-27.
Wang Z, Dooley TP, Curto EV, Davis RL, VandeBerg JL: Cross-species application of cDNA microarrays to profile gene expression using UV-induced melanoma in Monodelphis domestica as the model system. Genomics. 2004, 83: 588-599. 10.1016/j.ygeno.2003.10.007.
Medhora M, Bousamra M, Zhu D, Somberg L, Jacobs ER: Upregulation of collagens detected by gene array in a model of flow-induced pulmonary vascular remodeling. Am J Physiol Heart Circ Physiol. 2002, 282: H414-22.
Higgins MA, Berridge BR, Mills BJ, Schultze AE, Gao H, Searfoss GH, Baker TK, Ryan TP: Gene expression analysis of the acute phase response using a canine microarray. Toxicol Sci. 2003, 74: 470-484. 10.1093/toxsci/kfg142.
Ji W, Zhou W, Gregg K, Yu N, Davis S: A method for cross-species gene expression analysis with high-density oligonucleotide arrays. Nucleic Acids Res. 2004, 32: e93-10.1093/nar/gnh084.
Antipova AA, Tamayo P, Golub TR: A strategy for oligonucleotide microarray probe reduction. Genome Biol. 2002, 3: RESEARCH0073-10.1186/gb-2002-3-12-research0073.
Mei R, Hubbell E, Bekiranov S, Mittmann M, Christians FC, Shen MM, Lu G, Fang J, Liu WM, Ryder T, Kaplan P, Kulp D, Webster TA: Probe selection for high-density oligonucleotide arrays. Proc Natl Acad Sci U S A. 2003, 100: 11237-11242. 10.1073/pnas.1534744100.
Rajagopalan D: A comparison of statistical methods for analysis of high density oligonucleotide array data. Bioinformatics. 2003, 19: 1469-1476. 10.1093/bioinformatics/btg202.
Grigoryev DN, Ma SF, Irizarry RA, Ye SQ, Quackenbush J, Garcia JG: Orthologous gene-expression profiling in multi-species models: search for candidate genes. Genome Biol. 2004, 5: R34-10.1186/gb-2004-5-5-r34.
Simon BA, Easley BR, Grigoryev DN, Ma SF, Ye SQ, Lavoie T, Garcia JGN: Microarray analysis of regional cellular responses to local mechanical stress in experimental acute lung injury. (Submitted).
Ye SQ, Simon BA, Maloney JP, Zambelli-Weiner A, Gao L, Grant A, Easley RB, McVerry BJ, Tuder RM, Standiford T, Brower RG, Barnes KC, Garcia JG: Pre-B-Cell Colony-enhancing Factor as a Potential Novel Biomarker in Acute Lung Injury. Am J Respir Crit Care Med. 2005, 171: 361-370. 10.1164/rccm.200404-563OC.
McVerry BJ, Peng X, Hassoun PM, Sammani S, Simon BA, Garcia JG: Sphingosine 1-phosphate reduces vascular leak in murine and canine models of acute lung injury. Am J Respir Crit Care Med. 2004, 170: 987-993. 10.1164/rccm.200405-684OC.
Irizarry RA, Gautier L, Cope L: An R package for analyses of Affymetrix oligonucleotide arrays. The Analysis of Gene Expression Data: Methods and Software. Edited by: Parmigiani G, Garrett ES, Irizarry RA and Zeger SL. 2003, New York, Springer
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.
Doniger SW, Salomonis N, Dahlquist KD, Vranizan K, Lawlor SC, Conklin BR: MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biol. 2003, 4: R7-10.1186/gb-2003-4-1-r7.
We thank James MacDonald for helpful suggestions on Bioconductor application, and Nicholas Shank for organizing and publishing, pertinent to this paper, materials on the HopGene website. This work was supported by the NHLBI – sponsored HopGene Program in Genomics Application (HL-69340), SCCOR (HL-073994), and individual (DNG) NRSA grant (F32 HL74590-01A1).
DNG originated the probe level analysis of cross-species hybridization, carried out the bioinformatics studies and drafted the manuscript. SFM led the human component of the project and conducted human sample selection and hybridization analysis. BAS led the canine component of the project, originated and designed canine ventilation model, provided canine samples for hybridization, and revised the manuscript. RAI participated in the design of the study, established Bioconductor database for new canine GeneChip and consulted on the statistical analysis. SQY optimized protocols for canine samples hybridization to human and canine GeneChips, and assured the quality of generated microarray data. JGNC conceived of the study, supervised and coordinated the project, revised and advised on the manuscript preparation. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Grigoryev, D.N., Ma, SF., Simon, B.A. et al. In vitro identification and in silico utilization of interspecies sequence similarities using GeneChip® technology. BMC Genomics 6, 62 (2005). https://doi.org/10.1186/1471-2164-6-62
- Probe Pair
- Transcript Detection
- Canine Tissue
- Detectable Hybridization Signal
- 11th Probe