Using comparative genomic hybridization to survey genomic sequence divergence across species: a proof-of-concept from Drosophila
- Suzy CP Renn†1Email author,
- Heather E Machado†1,
- Albyn Jones2,
- Kosha Soneji3,
- Rob J Kulathinal4 and
- Hans A Hofmann5
© Renn et al; licensee BioMed Central Ltd. 2010
Received: 8 December 2009
Accepted: 29 April 2010
Published: 29 April 2010
Genome-wide analysis of sequence divergence among species offers profound insights into the evolutionary processes that shape lineages. When full-genome sequencing is not feasible for a broad comparative study, we propose the use of array-based comparative genomic hybridization (aCGH) in order to identify orthologous genes with high sequence divergence. Here we discuss experimental design, statistical power, success rate, sources of variation and potential confounding factors. We used a spotted PCR product microarray platform from Drosophila melanogaster to assess sequence divergence on a gene-by-gene basis in three fully sequenced heterologous species (D. sechellia, D. simulans, and D. yakuba). Because complete genome assemblies are available for these species this study presents a powerful test for the use of aCGH as a tool to measure sequence divergence.
We found a consistent and linear relationship between hybridization ratio and sequence divergence of the sample to the platform species. At higher levels of sequence divergence (< 92% sequence identity to D. melanogaster) ~84% of features had significantly less hybridization to the array in the heterologous species than the platform species, and thus could be identified as "diverged". At lower levels of divergence (≥ 97% identity), only 13% of genes were identified as diverged. While ~40% of the variation in hybridization ratio can be accounted for by variation in sequence identity of the heterologous sample relative to D. melanogaster, other individual characteristics of the DNA sequences, such as GC content, also contribute to variation in hybridization ratio, as does technical variation.
Here we demonstrate that aCGH can accurately be used as a proxy to estimate genome-wide divergence, thus providing an efficient way to evaluate how evolutionary processes and genomic architecture can shape species diversity in non-model systems. Given the increased number of species for which microarray platforms are available, comparative studies can be conducted for many interesting lineages in order to identify highly diverged genes that may be the target of natural selection.
Comparison of genomic DNA sequence among closely related strains or species is a powerful approach with which to identify heterogeneity in evolutionary processes such as selection, mutation rates, and rates of introgression, as well as to unmask phylogenetic relationships. However, even with the recent advances in DNA sequencing technology and rapidly dropping costs, complete genome sequence data are not readily available for many closely related eukaryotes that serve as model systems for organismal evolution [but see [1, 2]]. As an alternative, comparative genomic hybridization (CGH) offers a means to estimate sequence divergence.
Although the use of genomic DNA (gDNA) hybridization for phylogenetic analyses and genome-wide estimation of sequence similarity date to long before vast amounts of sequence data became available [e.g. [3, 4]], this approach has experienced a renaissance with the development of genomic tools, specifically microarrays. On a relatively coarse level, array-based CGH (aCGH) has been widely applied to identify chromosomal aberrations underlying cancer [for review see ]. When gDNA isolated from a tumor is competitively hybridized against gDNA isolated from normal tissue, genomic regions that have been deleted in the genome of the tumor cells will fail to hybridize to the array features while genomic regions that have been duplicated (amplified) in the genome of the tumor cells will hybridize at a ratio of 2:1 (or greater). At a finer level of resolution, modifications of this technique have allowed microarray-based genotyping of single nucleotide polymorphisms within and between populations [e.g. Arabidopsis: , stickleback fish: ]. Array-based techniques can also be applied to genome-scale comparisons between closely related species (or strains) in order to conduct a (nearly) complete analysis of sequence divergence on a gene-by-gene basis.
Unlike microarrays designed for genotyping known polymorphisms [reviewed by ] or re-sequencing [human: , Arabidopsis: ], microarrays designed for gene expression studies can also be used to compare the genomic content (in coding sequence) of closely related species. In a typical experiment, gDNA from the platform species (from which the microarray was constructed) is compared on the array to gDNA from another (heterologous) species of interest. This technique has been used to reveal genomic regions likely involved in an organism's ability to inhabit a specific environment [Chlamydia trachomatis tissue specificity: , Sinorhizobium meliloti root symbiont: , Clostridium difficile host specificity: ], pathogenicity [Yersinia pesits: [14, 15], Mycobacterium tuberculosis: , Vibrio cholerae: ], genomic duplications and deletions associated with population divergence and speciation [Anopheles gambiae: [18, 19]], and genomic regions that differentiate humans from other primate species [20, 21]. While most studies rely only on presence or absence metrics, a few studies have suggested that the relationship between hybridization signal ratio using aCGH and nucleotide identity is roughly log-linear [11, 22]. Using this relatively inexpensive approach, it is possible to identify rapidly evolving genes [Paxillus involutus: ] and in some cases lend insight to phylogenetic relationships [Shewanella: , Salmonella: , Saccharomyces: ]. While the majority of these examples derive from studies in microbes, the technique is amenable to genomes of any size. It must be noted, of course, that array-based comparisons do not reveal the actual genomic sequence for the novel genes of interest. Instead, an estimate of sequence identity is obtained at a price and effort far below that of whole genome sequencing.
In the present study, we examine the relationship between hybridization ratio and sequence divergence using a cDNA microarray constructed for D. melanogaster. The availability of complete genome assemblies  for Drosophila melanogaster as well as three other Drosophilid species, D. simulans, D. sechellia (both 2-3 MY diverged relative to D. melanogaster) and D. yakuba (10-15 MY diverged relative to D. melanogaster) provides us with a unique opportunity to demonstrate the degree to which hybridization ratio reflects underlying sequence divergence. It is not our goal to devise an explicit model to explain variation in hybridization due to other sequence characteristics; rather, we demonstrate success rate, discuss the effects of a few characteristics contributing to variation in hybridization kinetics, and provide an example for the use of aCGH that can be applied to non-model organisms.
We show that sequence divergence between orthologous genes can be successfully detected for closely and not so closely related species. Approximately 40% of the variation in gDNA hybridization ratios can be explained by sequence divergence, as measured by nucleotide dissimilarity between sequences. Other sequence-specific characteristics also explain part of the variation in hybridization ratio, and become more prominent with increased sequence divergence. Similarly, technical variation increases with increasing sequence divergence; however, this last source of variation can be overcome with increased replication. We demonstrate the potential for functional analysis and the generation of testable hypotheses based upon hybridization statistics for Gene Ontology annotation.
Results and Discussion
Detection of reduced hybridization
Features identified as diverged from D. melanogaster.
P < 0.1 FDR
P < 0.05 FDR
P < 0.01 FDR
Detection of sequence divergence
Relationship between sequence divergence and hybridization ratio
The majority of aCGH studies, even in microbes, aim to identify only presence or absence of orthologs. Such studies generally employ one of two assignment strategies. The first strategy employs a cut-off threshold derived from a separate characterized strain or other published results [32, 33]. The second strategy analyzes each hybridization dataset according to its intrinsic experimental variability [e.g. GACK: , GENCOM: ] in order to determine presence or absence. However, beyond a binary assignment, a more descriptive relationship between sequence divergence and hybridization ratio is possible. For example, Kim et al.  defined a transition zone for those genes thought to be present but highly diverged. For two genes, studied in seven microbial species, a linear relationship was found between the log hybridization ratio and percent divergence . Similar studies conducted with species of the genus Bartonella or strains of Chlamydia also found a linear relationship between log hybridization ratio and percent sequence divergence and noted deviation from the linearity of the relationship for orthologs with < 75 %ID, of which there were few in our study.
Sources of variation and statistical power
While technical variation in hybridization ratios exists even for within-species experiments (where sample and probe sequences are almost identical), previous work has suggested that this variation increases as sequence identity decreases [e.g. [22, 37]]; however, this relationship has not been fully investigated. Here, with the incorporation of additional replicates and the consideration of variation due to hybridization kinetics, we expand our understanding of the sources of variation that cause the hybridization ratio for individual features to deviate from the value predicted by the linear regression.
Increased technical replication leads to an increased power to detect array features that show a difference in the hybridization strength for D. yakuba vs. D. melanogaster
The magnitude of the variation due to technical error is measured by the SD of the fitted value for the hybridization ratio of each feature. The magnitude of the variation due to DNA hybridization kinetics is measured by the standard deviation (SD) of the residuals from the regression line (estimated as the median absolute residual divided by 0.6745). Based on the 8-array D. yakuba vs. D. melanogaster dataset, we found that hybridization ratios for features representing conserved orthologs (greater than ~95 %ID) were more affected by technical error, whereas those of more diverged orthologs (less than ~95 %ID) were more strongly influenced by DNA hybridization kinetics (Figure 4B). Both technical error and variation due to DNA hybridization kinetics increase with greater sequence divergence.
The regression model for %ID and sequence characteristics vs. hybridization ratio [hybridization ratio ~%ID * (%GC + length + % align))] includes the interaction of GC content, probe length, and percent alignment of D. melanogaster probe sequence to D. yakuba sequence.
It is also intuitive that percent alignment length would be positively correlated with the D. yakuba vs. D. melanogaster hybridization ratio, as it is another measure of divergence that is not taken into consideration by %ID. A complete alignment indicates low divergence, whereas an incomplete alignment indicates substantial divergence. The absence of a strong correlation between relative alignment length and hybridization ratio may be due to insertions in the heterologous sequence that result in reduced hybridization strength and inflated relative alignment length (even above 100%). Our goal in discussing these variables is to demonstrate the complexity of hybridization kinetics without devising a precise metric. For that reason we did not include other possible explanatory variables, such as insertions, deletions, GC content of the heterologous sample, or even the presence of paralogs, in our discussion. These factors undoubtedly will contribute to studies in non-model species; however, without full genome sequence, little benefit would be gained from an explicit model. We see aCGH as an efficient and inexpensive method for identifying highly diverged genes among species for which little or no genomic sequence information is available. Direct sequence analysis (in multiple individuals) would be necessary to further investigate any genes of interest identified by this method (e.g., synonymous vs. non-synonymous substitutions).
Example analysis 1: constraint within Gene Ontology categories
We tested for under-representation of GO categories among the genes determined by aCGH to be diverged in D. yakuba (D. melanogaster-biased). Of the 7009 D. melanogaster-biased features with FlyBase gene IDs (representing 6094 genes), approximately half had GO annotations. Several GO categories were found to be under-represented, indicating increased sequence conservation and possible constraint of these processes, functions, and components in the D. yakuba lineage (see Additional file 1; representatives with asterisk on Figure 6). Molecular functions that were found to be under-represented in the set of diverged genes include GTPase and ATPase activities, NADH dehydrogenase activity, hydrogen ion transport activity, and acid-amino acid ligase activity as well as general ligase activity. Underrepresented biological processes include translation, ATP synthesis-coupled electron transport, oxidative phosphorylation, protein transport, actin-filament-based processes, cell-cell signaling, including synaptic transmission, and several processes associated with development (instar larval or pupal development, organ development and organ morphogenesis, eye development, and sensory organ development more generally). Finally, underrepresented cellular components include those associated with the proton-transporting ATPase, cytosol, ribosome, lipid particle, clathrin-coated vesicle, cytoplasmic vesicle, and the synapse. Overall, 111 out of 4874 GO categories that were tested (2.3%) were under-represented (i.e. more conserved) than would have been expected by chance (P < 0.01). As demonstrated in this example analysis, categories involved in very fundamental processes (ie. protein synthesis and transport, basic metabolic activity, development, and synaptic function) are likely to be more constrained.
In this example, we used GO annotations specifically for Drosophila genes. However, this strategy can also be applied to any non-traditional model organism (or group of species) for which there is an array platform available, based on DNA sequence similarity between array features and annotations (which are primarily derived from standard model organisms). Despite the ascertainment bias inherent to GO terms, which is largely due to the nature of hypothesis-driven research conducted in model organisms, this approach allows rigorous statistical analysis for over- and under-representation of particular molecular functions, biological processes, and cellular components [for example, see ]. Also, applying GO terms reduces the complexity of the data, avoids experimenter bias in cross-referencing between experiments and species, and facilitates comparisons of experimental results obtained in different organisms and/or with different platforms.
Example analysis 2: highly diverged genes
Genes that have high divergence among species can be identified as those with substantially reduced hybridization ratios. We focus on array features with greater than four-fold reduction in hybridization ratio for D. yakuba. Among these 1686 highly D. melanogaster biased features, representing 145 genes with GO annotations, six GO categories were found to be overrepresented, including peptidase activity, serine hydrolase activity, sensory perception of taste, and three parent GO terms (P < 0.01). The enrichment for D. melanogaster-biased features suggests that genes in these categories are either highly diverged or are deleted in D. yakuba.
Representative GO annotations for genes highly diverged in all test species.
transcription factor activity
choline O-acetyltransferase activity
protein serine/threonine kinase activity
Casein kinase II alpha subunit
protein serine/threonine kinase activity
Protein kinase 61C
serine-type endopeptidase inhibitor activity
Serine protease inhibitor 3
serine-type endopeptidase inhibitor activity
serine-type endopeptidase inhibitor activity
ionotropic glutamate receptor activity
Ionotropic receptor 41a
guanyl-nucleotide exchange factor activity
structural constituent of cytoskeleton
voltage-gated calcium channel activity
ubiquitin-dependent protein catabolic process
G-protein coupled receptor protein signaling
Salivary gland secretion 3
cysteine-type peptidase activity
cysteine-type peptidase activity
zinc ion binding
zinc ion binding
electron carrier activity
cell wall catabolic process
Accessory gland-specific peptide 26Aa
The results presented here demonstrate that aCGH can robustly detect genes that are highly diverged in a given species compared with one for which a microarray platform is available. Our use of a D. melanogaster microarray to estimate sequence divergence on a gene-by-gene basis for three fully sequenced heterologous species allowed for a proof-of-principle for this approach and allowed us to explore the success and biases inherent in this technique. We found a consistent and linear relationship between array hybridization ratio and sequence divergence between the sample and the platform species. The limit of detectable sequence divergence depends on the power of the experiment. While the power can be increased with additional technical replicates, there will still be a subset of diverged genes that escape detection due to other specific hybridization kinetics of the sample DNA to the array feature. This technique is generally applicable, even though thresholds, correlation strengths, and appropriate divergence distances may differ for different array platforms (cDNAs, long or short oligonucleotides). As the number of microarray platforms available for non-traditional model species is continually increasing [axolotl: , dolphin: , butterfly: [41, 42], coral: , many fish species, for review: , crustaceans: [45–47]], researchers focusing on these and related model systems will continue to benefit from the relatively low cost of array hybridizations. Rapid advances in next-generation sequencing technology notwithstanding [2, 48], aCGH provides an effective alternative to de novo sequencing of a large number of complex eukaryotic genomes.
We used a Drosophila melanogaster microarray with ~22,000 features containing PCR products (~500 base pairs long) generated from custom primers designed to predict open reading frames [; GEO platform number GPL6056]. The microarray was printed on poly-L-lysine slides (Thermo Scientific) in a 48 pin format using an OmniGrid-100 arrayer (GeneMachines). Following hydration, snap drying and UV cross-linking, the slides were blocked with succinic anhydride and sodium borate in 1-Methyl-2-Pyrrolidinone, rinsed, dried according to standard procedure  and stored dry until used.
Sample Preparation and aCGH
Isogenic Drosophila melanogaster, D. simulans, D. sechellia and D. yakuba strains (Dmel\y;cn;bw;sp, Dsim\w , Dsec\Robertson3C, Dyak Tai18E2) were obtained from the Tuscon Drosophila stock center (now known as the San Diego Drosophila Species Stock Center). Genomic DNA was isolated from ~100 Drosophila males of each stock according to a standard ProteinaseK/Phenol:Chloroform protocol. DNA quantity and purity was assayed (via Nanodrop 1000) prior to and after DNA size reduction using a Hydroshear (Genome Solutions/Digilab) with a standard orifice set to maximal possible shearing speed (13) for 20 cycles (maximal shearing speed varies with individual orifice). This treatment resulted in fragments of 500 bp - 2 kb as determined by gel electrophoresis, visualized with ethidium bromide. Two micrograms of sheared genomic DNA was fluorescently labeled through incorporation of Cy3 or Cy5 labeled dCTP (Amersham) in a Klenow fragment (Invitrogen; Bioprime) reaction of 35 microliters according to manufacturer's protocol. Labeled sample DNA was purified by size exclusion on YM-30 filters (Eppendorf) and appropriate samples were combined. Hybridizations proceeded for ~16 hours at 65°C in a 3.4× SSC, 0.15% SDS, 1 mM DTT hybridization buffer. Male D. melanogaster samples were used in competitive hybridizations with two male D. sechellia samples, two male D. simulans samples and two male D. yakuba samples, incorporating dye swaps to account for dye bias. These aCGH hybridizations were analyzed for the ability to detect significantly diverged genes. An additional six D. melanogaster versus D. yakuba aCGH hybridizations were available in order to asses the effect of increased technical replication. For this power analysis, only genes located on the autosomes were used because a subset of the hybridizations involved D. yakuba female genomic DNA of the same strain.
Microarray Data Analysis
Hybridized arrays were scanned with an Axon 4000B scanner (Axon Instruments) using Genepix 5.0 software (Axon Instruments). All raw array data have been submitted to GEO database (dataseries number GSE18416 sample number GSM459056-67). Features of poor quality (signal intensity < 2 standard deviations above background) and those of potentially erroneous sequences (mismatch between initial PCR product sequence prediction and current D. melanogaster database; refseq_rna 12/2008) were excluded. Features were only considered in the analysis if they survived these technical filters on multiple arrays for a given species comparison. Raw data from Genepix was imported into R, and LIMMA  was used to apply a background correction ("minimum") and within-array intensity normalization ("loess"). Because we expect the normalization of cross-species arrays to be affected by a substantial number of diverged genes in the non-platform species , we performed the within-array normalization using a set of ~1000 genes highly conserved (greater than ~95% sequence identity; determined with NCBI BLAST to Genbank) among D. melanogaster D. simulans and D. yakuba. A linear model was fitted to the data using "lmFit", and "eBayes" provided error shrinkage towards a pooled estimate of variation . Array features were tested for hybridization ratios that were significantly different from equal as assessed after a FDR multiple testing correction at P < 0.1 . The Gene Ontology analysis was conducted in GOstats , using features with FlyBase IDs, as listed in the GEO platform.
Genomic Sequence Divergence
The sequences of the D. melanogaster probes were predicted by blasting primers from the D. melanogaster probe (GEO Profiles accession: GPL6056) against the D. melanogaster Release 5 assembly and searching for unique and proximal (within 600 base pairs of each other) targets. We queried the resulting sequences against the D. simulans and D. yakuba NCBI genomes (chromosome) using "megablast" (2009) and against the full chromosome sequence assemblies for D. sechellia downloaded from flybase.org (release 1.3). From each heterologous genome, the top BLAST hit to each array feature (threshold e-value 10-14) was used to obtain the percent similarity between the two sequences and the length of the alignment.
We are grateful for the support of Dan Hartl (Harvard University) and members of his lab for assistance with microarray production and annotation. This work was supported by the Murdock Life Trust Foundation (S.C.P.R.) and by the Bauer Center for Genomics Research (H.A.H.).
- Rokas A, Abbot P: Harnessing genomics for evolutionary insights. Trends Ecol Evol. 2009, 24 (4): 192-200. 10.1016/j.tree.2008.11.004.PubMedView ArticleGoogle Scholar
- Turner DJ, Keane TM, Sudbery I, Adams DJ: Next-generation sequencing of vertebrate experimental organisms. Mamm Genome. 2009, 20 (6): 327-338. 10.1007/s00335-009-9187-4.PubMed CentralPubMedView ArticleGoogle Scholar
- Sibley CG, Ahlquist JE: The phylogeny of the hominoid primates, as indicated by DNA-DNA hybridization. J Mol Evol. 1984, 20 (1): 2-15. 10.1007/BF02101980.PubMedView ArticleGoogle Scholar
- Templeton AR: The Phylogeny of the hominoid primates - a statistical-analysis of the DNA-DNA hybridization data. Mol Biol Evol. 1985, 2 (5): 420-433.PubMedGoogle Scholar
- Pinkel D, Albertson DG: Comparative genomic hybridization. Annu Rev Genomics Hum Genet. 2005, 6: 331-354. 10.1146/annurev.genom.6.080604.162140.PubMedView ArticleGoogle Scholar
- West MAL, van Leeuwen H, Kozik A, Kliebenstein DJ, Doerge RW, St Clair DA, Michelmore RW: High-density haplotyping with microarray-based expression and single feature polymorphism markers in Arabidopsis. Genome Res. 2006, 16 (6): 787-795. 10.1101/gr.5011206.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller MR, Dunham JP, Amores A, Cresko WA, Johnson EA: Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res. 2007, 17 (2): 240-248. 10.1101/gr.5681207.PubMed CentralPubMedView ArticleGoogle Scholar
- Fan JB, Chee MS, Gunderson KL: Highly parallel genomic assays. Nat Rev Genet. 2006, 7 (8): 632-644. 10.1038/nrg1901.PubMedView ArticleGoogle Scholar
- Clark RM, Schweikert G, Toomajian C, Ossowski S, Zeller G, Shinn P, Warthmann N, Hu TT, Fu G, Hinds DA, Chen H, Frazer KA, Huson DH, Schölkopf B, Nordborg M, Rätsch G, Ecker JR, Weigel D: Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana. Science. 2007, 317 (5836): 338-342. 10.1126/science.1138632.PubMedView ArticleGoogle Scholar
- Hinds DA, Stuve LL, Nilsen GB, Halperin E, Eskin E, Ballinger DG, Frazer KA, Cox DR: Whole-genome patterns of common DNA variation in three human populations. Science. 2005, 307 (5712): 1072-1079. 10.1126/science.1105436.PubMedView ArticleGoogle Scholar
- Brunelle BW, Nicholson TL, Stephens RS: Microarray-based genomic surveying of gene polymorphisms in Chlamydia trachomatis. Genome Biol. 2004, 5 (6): 10.1186/gb-2004-5-6-r42.Google Scholar
- Giuntini E, Mengoni A, De Filippo C, Cavalieri D, Aubin-Horth N, Landry CR, Becker A, Bazzicalupo M: Large-scale genetic variation of the symbiosis-required megaplasmid pSymA revealed by comparative genomic analysis of Sinorhizobium meliloti natural strains. BMC Genomics. 2005, 6: 158-10.1186/1471-2164-6-158.PubMed CentralPubMedView ArticleGoogle Scholar
- Janvilisri T, Scaria J, Thompson AD, Nicholson A, Limbago BM, Arroyo LG, Songer JG, Grohn YT, Chang YF: Microarray identification of Clostridium difficile core components and divergent regions associated with host origin. J Bacteriol. 2009, 191 (12): 3881-3891. 10.1128/JB.00222-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhou D, Han Y, Song Y, Tong Z, Wang J, Guo Z, Pei D, Pang X, Zhai J, Li M, Cui B, Qi Z, Jin L, Dai R, Du Z, Bao J, Zhang X, Yu J, Wang J, Huang P, Yang R: DNA microarmy analysis of genome dynamics in Yersinia pestis: Insights into bacterial genome microevolution and niche adaptation. J Bacteriol. 2004, 186 (15): 5138-5146. 10.1128/JB.186.15.5138-5146.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Hinchliffe SJ, Isherwood KE, Stabler RA, Prentice MB, Rakin A, Nichols RA, Oyston PCF, Hinds J, Titball RW, Wren BW: Application of DNA mcroarrays to study the evolutionary genomics of Yersinia pestis and Yersinia pseudotuberculosis. Genome Res. 2003, 13 (9): 2018-2029. 10.1101/gr.1507303.PubMed CentralPubMedView ArticleGoogle Scholar
- Kato-Maeda M, Rhee JT, Gingeras TR, Salamon H, Drenkow J, Smittipat N, Small PM: Comparing genomes within the species Mycobacterium tuberculosis. Genome Res. 2001, 11 (4): 547-554. 10.1101/gr.166401.PubMed CentralPubMedView ArticleGoogle Scholar
- Dziejman M, Balon E, Boyd D, Fraser CM, Heidelberg JF, Mekalanos JJ: Comparative genomic analysis of Vibrio cholerae: Genes that correlate with cholera endemic and pandemic disease. Proc Natl Acad Sci USA. 2002, 99 (3): 1556-1561. 10.1073/pnas.042667999.PubMed CentralPubMedView ArticleGoogle Scholar
- Turner TL, Hahn MW, Nuzhdin SV: Genomic islands of speciation in Anopheles gambiae. PLoS Biol. 2005, 3 (9): 10.1371/journal.pbio.0030285.Google Scholar
- Riehle MM, Markianos K, Niaré O, Xu J, Li J, Touré AM, Podiougou B, Oduol F, Diawara S, Diallo M, Coulibaly B, Ouatara A, Kruglyak L, Traoré SF, Vernick KD: Natural malaria infection in Anopheles gambiae is regulated by a single genomic control region. Science. 2006, 312 (5773): 577-579. 10.1126/science.1124153.PubMedView ArticleGoogle Scholar
- Locke DP, Segraves R, Carbone L, Archidiacono N, Albertson DG, Pinkel D, Eichler EE: Large-scale variation among human and great ape genomes determined by array comparative genomic hybridization. Genome Res. 2003, 13 (3): 347-357. 10.1101/gr.1003303.PubMed CentralPubMedView ArticleGoogle Scholar
- Fortna A, Kim Y, MacLaren E, Marshall K, Hahn G, Meltesen L, Brenton M, Hink R, Burgers S, Hernandez-Boussard T, Karimpour-Fard A, Glueck D, McGavran L, Berry R, Pollack J, Sikela JM: Lineage-specific gene duplication and loss in human and great ape evolution. PLoS Biol. 2004, 2 (7): 937-954. 10.1371/journal.pbio.0020207.View ArticleGoogle Scholar
- Taboada EN, Acedillo RR, Luebbert CC, Findlay WA, Nash JHE: A new approach for the analysis of bacterial microarray-based comparative genomic hybridization: insights from an empirical study. BMC Genomics. 2005, 6: 1-10. 10.1186/1471-2164-6-78.View ArticleGoogle Scholar
- Le Quere A, Eriksen KA, Rajashekar B, Schutzenbubel A, Canback B, Johansson T, Tunlid A: Screening for rapidly evolving genes in the ectomycorrhizal fungus Paxillus involutus using cDNA microarrays. Mol Ecol. 2006, 15 (2): 535-550. 10.1111/j.1365-294X.2005.02796.x.PubMedView ArticleGoogle Scholar
- Murray AE, Lies D, Li G, Nealson K, Zhou J, Tiedje JM: DNA/DNA hybridization to microarrays reveals gene-specific differences between closely related microbial genomes. Proc Natl Acad Sci USA. 2001, 98 (17): 9853-9858. 10.1073/pnas.171178898.PubMed CentralPubMedView ArticleGoogle Scholar
- Porwollik S, Wong RMY, McClelland M: Evolutionary genomics of Salmonella: Gene acquisitions revealed by microarray analysis. Proc Natl Acad Sci USA. 2002, 99 (13): 8956-8961. 10.1073/pnas.122153699.PubMed CentralPubMedView ArticleGoogle Scholar
- Edwards-Ingram LC, Gent ME, Hoyle DC, Hayes A, Stateva LI, Oliver SG: Comparative genomic hybridization provides new insights into the molecular taxonomy of the Saccharomyces sensu stricto complex. Genome Res. 2004, 14 (6): 1043-1051. 10.1101/gr.2114704.PubMed CentralPubMedView ArticleGoogle Scholar
- Drosophila 12 Genomes Consortium, Clark AG, Eisen MB, Smith DR, Bergman CM, Oliver B, Markow TA, Kaufman TC, Kellis M, Gelbart W, Iyer VN, Pollard DA, Sackton TB, Larracuente AM, Singh ND, Abad JP, Abt DN, Adryan B, Aguade M, Akashi H, Anderson WW, Aquadro CF, Ardell DH, Arguello R, Artieri CG, Barbash DA, Barker D, Barsanti P, Batterham P, Batzoglou S, Begun D: Evolution of genes and genomes on the Drosophila phylogeny. Nature. 2007, 450 (7167): 203-218. 10.1038/nature06341.View ArticleGoogle Scholar
- O'Grady PM, Markow TA: Phylogenetic taxonomy in Drosophila. Fly (Austin). 2009, 3 (1): 1-5.View ArticleGoogle Scholar
- Fitzpatrick DA, Logue ME, Stajich JE, Butler G: A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis. BMC Evol Biol. 2006, 6: 99-10.1186/1471-2148-6-99.PubMed CentralPubMedView ArticleGoogle Scholar
- McQuiston JR, Herrera-Leon S, Wertheim BC, Doyle J, Fields PI, Tauxe RV, Logsdon JM: Molecular phylogeny of the Salmonellae: Relationships among Salmonella species and subspecies determined from four housekeeping genes and evidence of lateral gene transfer events. J Bacteriol. 2008, 190 (21): 7060-7067. 10.1128/JB.01552-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Townsend JP: Resolution of large and small differences in gene expression using models for the Bayesian analysis of gene expression levels and spotted DNA microarrays. BMC Bioinformatics. 2004, 5: 13-10.1186/1471-2105-5-54.View ArticleGoogle Scholar
- Hatfield GW, Baldi P: DNA microarrays and gene expression: from experiments to data analysis and modeling. 2002, Cambridge: Cambridge University PressGoogle Scholar
- Dopman EB, Hartl DL: A portrait of copy-number polymorphism in Drosophila melanogaster. Proc Natl Acad Sci USA. 2007, 104 (50): 19920-19925. 10.1073/pnas.0709888104.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim CC, Joyce EA, Chan K, Falkow S: Improved analytical methods for microarray-based genome-composition analysis. Genome Biol. 2002, 3 (11): RESEARCH0065-10.1186/gb-2002-3-11-research0065.PubMed CentralPubMedView ArticleGoogle Scholar
- Pin C, Reuter M, Pearson B, Friis L, Overweg K, Baranyi J, Wells J: Comparison of different approaches for comparative genetic analysis using microarray hybridization. Appl Microbiol Biotechnol. 2006, 72 (4): 852-859. 10.1007/s00253-006-0536-x.PubMedView ArticleGoogle Scholar
- Lindroos HL, Mira A, Repsilber D, Vinnere O, Naslund K, Dehio M, Dehio C, Andersson SGE: Characterization of the genome composition of Bartonella koehlerae by microarray comparative genomic hybridization profiling. J Bacteriol. 2005, 187 (17): 6155-6165. 10.1128/JB.187.17.6155-6165.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Gilad Y, Rifkin SA, Bertone P, Gerstein M, White KP: Multi-species microarrays reveal the effect of sequence divergence on gene expression profiles. Genome Res. 2005, 15 (5): 674-680. 10.1101/gr.3335705.PubMed CentralPubMedView ArticleGoogle Scholar
- Renn SCP, Aubin-Horth N, Hofmann HA: Fish and chips: functional genomics of social plasticity in an African cichlid fish. J Exp Biol. 2008, 211 (18): 3041-3056. 10.1242/jeb.018242.PubMed CentralPubMedView ArticleGoogle Scholar
- Page RB, Voss SR, Samuels AK, Smith JJ, Putta S, Beachy CK: Effect of thyroid hormone concentration on the transcriptional response underlying induced metamorphosis in the Mexican axolotl (Ambystoma). BMC Genomics. 2008, 9: 78-10.1186/1471-2164-9-78.PubMed CentralPubMedView ArticleGoogle Scholar
- Mancia A, Lundqvist ML, Romano TA, Peden-Adams MM, Fair PA, Kindy MS, Ellis BC, Gattoni-Celli S, McKillen DJ, Trent HF, Chen YA, Almeida JS, Gross PS, Chapman RW, Warr GW: A dolphin peripheral blood leukocyte cDNA microarray for studies of immune function and stress reactions. Dev Comp Immunol. 2007, 31 (5): 520-529. 10.1016/j.dci.2006.07.011.PubMedView ArticleGoogle Scholar
- Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, Marden JH: Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Mol Ecol. 2008, 17 (7): 1636-1647. 10.1111/j.1365-294X.2008.03666.x.PubMedView ArticleGoogle Scholar
- Zhu H, Gegear RJ, Casselman A, Kanginakudru S, Reppert SM: Defining behavioral and molecular differences between summer and migratory monarch butterflies. BMC Biol. 2009, 7: 14-10.1186/1741-7007-7-14.PubMed CentralPubMedView ArticleGoogle Scholar
- Edge SE, Morgan MB, Gleason DF, Snell TW: Development of a coral cDNA array to examine gene expression profiles in Montastraea faveolata exposed to environmental stress. Marine Pollution Bulletin. 2005, 51 (5-7): 507-523. 10.1016/j.marpolbul.2005.07.007.PubMedView ArticleGoogle Scholar
- Kassahn KS: Microarrays for comparative and ecological genomics: beyond single-species applications of array technologies. Journal of Fish Biology. 2008, 72 (9): 2407-2434. 10.1111/j.1095-8649.2008.01890.x.View ArticleGoogle Scholar
- Dhar AK, Dettori A, Roux MM, Klimpel KR, Read B: Identification of differentially expressed genes in shrimp (Penaeus stylirostris) infected with White spot syndrome virus by cDNA microarrays. Arch Virol. 2003, 148 (12): 2381-2396. 10.1007/s00705-003-0172-z.PubMedView ArticleGoogle Scholar
- Stillman JH, Teranishi KS, Tagmount A, Lindquist EA, Brokstein PB: Construction and characterization of EST libraries from the porcelain crab, Petrolisthes cinctipes. Integr Comp Biol. 2006, 46 (6): 919-930. 10.1093/icb/icl007.PubMedView ArticleGoogle Scholar
- Ki JS, Raisuddin S, Lee KW, Hwang DS, Han J, Rhee JS, Kim IC, Park HG, Ryu JC, Lee JS: Gene expression profiling of copper-induced responses in the intertidal copepod Tigriopus japonicus using a 6 K oligochip microarray. Aquat Toxicol. 2009, 93 (4): 177-187. 10.1016/j.aquatox.2009.04.004.PubMedView ArticleGoogle Scholar
- Shendure J, Ji HL: Next-generation DNA sequencing. Nat Biotechnol. 2008, 26 (10): 1135-1145. 10.1038/nbt1486.PubMedView ArticleGoogle Scholar
- Hegde P, Qi R, Abernathy K, Gay C, Dharap S, Gaspard R, Hughes JE, Snesrud E, Lee N, Quackenbush J: A concise guide to cDNA microarray analysis. Biotechniques. 2000, 29 (3): 548-PubMedGoogle Scholar
- Smyth GK: Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420. full_text.View ArticleGoogle Scholar
- van Hijum S, Baerends RJS, Zomer AL, Karsens HA, Martin-Requena V, Trelles O, Kok J, Kuipers OP: Supervised Lowess normalization of comparative genome hybridization data - application to lactococcal strain comparisons. BMC Bioinformatics. 2008, 9: 93-10.1186/1471-2105-9-93.PubMed CentralPubMedView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: 1-26.Google Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995, 57 (1): 289-300.Google Scholar
- Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20 (9): 1464-1465. 10.1093/bioinformatics/bth088.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.