Differential DNA methylation associated with multiple sclerosis and disease modifying treatments in an underrepresented minority population.

Black and Hispanic American patients frequently develop earlier onset of multiple sclerosis (MS) and a more severe disease course that can be resistant to disease modifying treatments. The objectives were to identify differential methylation of genomic DNA (gDNA) associated with disease susceptibility and treatment responses in a cohort of MS patients from underrepresented minority populations. Patients with MS and controls with non-inflammatory neurologic conditions were consented and enrolled under an IRB-approved protocol. Approximately 64% of donors identified as Black or African American and 30% as White, Hispanic-Latino. Infinium MethylationEPIC bead arrays were utilized to measure epigenome-wide gDNA methylation of whole blood. Data were analyzed in the presence and absence of adjustments for unknown covariates in the dataset, some of which corresponded to disease modifying treatments. Global patterns of differential methylation associated with MS were strongest for those probes that showed relative demethylation of loci with lower M values. Pathway analysis revealed unexpected associations with shigellosis and amoebiasis. Enrichment analysis revealed an over-representation of probes in enhancer regions and an under-representation in promoters. In the presence of adjustments for covariates that included disease modifying treatments, analysis revealed 10 differentially methylated regions (DMR's) with an FDR <1E-77. Five of these genes (ARID5B, BAZ2B, RABGAP1, SFRP2, WBP1L) are associated with cancer risk and cellular differentiation and have not been previously identified in MS studies. Hierarchical cluster and multi-dimensional scaling analysis of differential DNA methylation at 147 loci within those DMR's was sufficient to differentiate MS donors from controls. In the absence of corrections for disease modifying treatments, differential methylation in patients treated with dimethyl fumarate was associated with immune regulatory pathways that regulate cytokine and chemokine signaling, axon guidance, and adherens junctions. These results demonstrate possible associations of gastrointestinal pathogens and regulation of cellular differentiation with MS susceptibility in our patient cohort. This work further suggests that analyses can be performed in the presence and absence of corrections for immune therapies. Because of their high representation in our patient cohort, these results may be of specific relevance in the regulation of disease susceptibility and treatment responses in Black and Hispanic Americans.

Construction of a microsatellites-based linkage map for the white grouper (Epinephelus aeneus).

The white grouper (Epinephelus aeneus) is a promising candidate for domestication and aquaculture due to its fast growth, excellent taste, and high market price. A linkage map is an essential framework for mapping quantitative trait loci for economic traits and the study of genome evolution. DNA of a single individual was deep-sequenced, and microsatellite markers were identified in 177 of the largest scaffolds of the sequence assembly. The success rate of developing polymorphic homologous markers was 94.9% compared with 63.1% of heterologous markers from other grouper species. Of the 12 adult mature fish present in the broodstock tank, two males and two females were identified as parents of the assigned offspring by parenthood analysis using 34 heterologous markers. A single full-sib family of 48 individuals was established for the construction of first-generation linkage maps based on genotyping data of 222 microsatellites. The markers were assigned to 24 linkage groups in accordance to the 24 chromosomal pairs. The female and male maps consisting of 203 and 202 markers spanned 1053 and 886 cM, with an average intermarker distance of 5.8 and 5.0 cM, respectively. Mapping of markers to linkage groups ends was enriched by using markers originating from scaffolds harboring telomeric repeat-containing RNA. Comparative mapping showed high synteny relationships among the white grouper, kelp grouper (E. bruneus), orange-spotted grouper (E. coioides), and Nile tilapia (Oreochromis niloticus). Thus, it would be useful to integrate the markers that were developed for different groupers, depending on sharing of sequence data, into a comprehensive consensus map.

Molecular heterochrony and the evolution of sociality in bumblebees (Bombus terrestris).

Sibling care is a hallmark of social insects, but its evolution remains challenging to explain at the molecular level. The hypothesis that sibling care evolved from ancestral maternal care in primitively eusocial insects has been elaborated to involve heterochronic changes in gene expression. This elaboration leads to the prediction that workers in these species will show patterns of gene expression more similar to foundress queens, who express maternal care behaviour, than to established queens engaged solely in reproductive behaviour. We tested this idea in bumblebees (Bombus terrestris) using a microarray platform with approximately 4500 genes. Unlike the wasp Polistes metricus, in which support for the above prediction has been obtained, we found that patterns of brain gene expression in foundress and queen bumblebees were more similar to each other than to workers. Comparisons of differentially expressed genes derived from this study and gene lists from microarray studies in Polistes and the honeybee Apis mellifera yielded a shared set of genes involved in the regulation of related social behaviours across independent eusocial lineages. Together, these results suggest that multiple independent evolutions of eusociality in the insects might have involved different evolutionary routes, but nevertheless involved some similarities at the molecular level.

Social regulation of maternal traits in nest-founding bumble bee (Bombus terrestris) queens.

During the nest-founding phase of the bumble bee colony cycle, queens undergo striking changes in maternal care behavior. Early in the founding phase, prior to the emergence of workers in the nest, queens are reproductive and also provision and feed their offspring. However, later in the founding phase, queens reduce their feeding of larvae and become specialized on reproduction. This transition is synchronized with the emergence of workers in the colony, who assume the task of feeding their siblings. Using a social manipulation experiment with the bumble bee Bombus terrestris, we tested the hypothesis that workers regulate the transition from feeding brood to specialization on reproduction in nest-founding bumble bee queens. Consistent with this hypothesis, we found that early-stage nest-founding queens with workers prematurely added to their nests reduce their brood-feeding behavior and increase egg laying, and likewise, late-stage nest-founding queens increase their brood-feeding behavior and decrease egg-laying when workers are removed from their nests. Further, brood-feeding and egg-laying behaviors were negatively correlated. We used Agilent microarrays designed from B. terrestris brain expressed sequenced tags (ESTs) to explore a second hypothesis, that workers alter brain gene expression in nest-founding queens. We found evidence that brain gene expression in nest-founding queens is altered by the presence of workers, with the effect being much stronger in late-stage founding queens. This study provides new insights into how the transition from feeding brood to specialization on reproduction in queen bumble bees is regulated during the nest initiation phase of the colony cycle.

Whole-genome resequencing of two elite sires for the detection of haplotypes under selection in dairy cattle.

Using a combination of whole-genome resequencing and high-density genotyping arrays, genome-wide haplotypes were reconstructed for two of the most important bulls in the history of the dairy cattle industry, Pawnee Farm Arlinda Chief ("Chief") and his son Walkway Chief Mark ("Mark"), each accounting for ∼7% of all current genomes. We aligned 20.5 Gbp (∼7.3× coverage) and 37.9 Gbp (∼13.5× coverage) of the Chief and Mark genomic sequences, respectively. More than 1.3 million high-quality SNPs were detected in Chief and Mark sequences. The genome-wide haplotypes inherited by Mark from Chief were reconstructed using ∼1 million informative SNPs. Comparison of a set of 15,826 SNPs that overlapped in the sequence-based and BovineSNP50 SNPs showed the accuracy of the sequence-based haplotype reconstruction to be as high as 97%. By using the BovineSNP50 genotypes, the frequencies of Chief alleles on his two haplotypes then were determined in 1,149 of his descendants, and the distribution was compared with the frequencies that would be expected assuming no selection. We identified 49 chromosomal segments in which Chief alleles showed strong evidence of selection. Candidate polymorphisms for traits that have been under selection in the dairy cattle population then were identified by referencing Chief's DNA sequence within these selected chromosome blocks. Eleven candidate genes were identified with functions related to milk-production, fertility, and disease-resistance traits. These data demonstrate that haplotype reconstruction of an ancestral proband by whole-genome resequencing in combination with high-density SNP genotyping of descendants can be used for rapid, genome-wide identification of the ancestor's alleles that have been subjected to artificial selection.

Transcriptome profiling of the small intestinal epithelium in germfree versus conventional piglets.

BACKGROUND: To gain insight into host-microbe interactions in a piglet model, a functional genomics approach was used to address the working hypothesis that transcriptionally regulated genes associated with promoting epithelial barrier function are activated as a defensive response to the intestinal microbiota. Cesarean-derived germfree (GF) newborn piglets were colonized with adult swine feces, and villus and crypt epithelial cell transcriptomes from colonized and GF neonatal piglets were compared using laser-capture microdissection and high-density porcine oligonucleotide microarray technology. RESULTS: Consistent with our hypothesis, resident microbiota induced the expression of genes contributing to intestinal epithelial cell turnover, mucus biosynthesis, and priming of the immune system. Furthermore, differential expression of genes associated with antigen presentation (pan SLA class I, B2M, TAP1 and TAPBP) demonstrated that microbiota induced immune responses using a distinct regulatory mechanism common for these genes. Specifically, gene network analysis revealed that microbial colonization activated both type I (IFNAR) and type II (IFNGR) interferon receptor mediated signaling cascades leading to enhanced expression of signal transducer and activator of transcription 1 (STAT1), STAT2 and IFN regulatory factor 7 (IRF7) transcription factors and the induction of IFN-inducible genes as a reflection of intestinal epithelial inflammation. In addition, activated RNA expression of NF-kappa-B inhibitor alpha (NFkappaBIA; a.k.a I-kappa-B-alpha, IKBalpha) and toll interacting protein (TOLLIP), both inhibitors of inflammation, along with downregulated expression of the immunoregulatory transcription factor GATA binding protein-1 (GATA1) is consistent with the maintenance of intestinal homeostasis. CONCLUSION: This study supports the concept that the intestinal epithelium has evolved to maintain a physiological state of inflammation with respect to continuous microbial exposure, which serves to sustain a tight intestinal barrier while preventing overt inflammatory responses that would compromise barrier function.

Discovery and characterization of 91 novel transcripts expressed in cattle placenta.

BACKGROUND: Among the eutherian mammals, placental architecture varies to a greater extent than any other tissue. The diversity of placental types, even within a single mammalian order suggests that genes expressed in placenta are under strong Darwinian selection. Thus, the ruminant placenta may be a rich source of genes to explore adaptive evolutionary responses in mammals. The aim of our study was to identify novel transcripts expressed in ruminant placenta, and to characterize them with respect to their expression patterns, organization of coding sequences in the genome, and potential functions. RESULTS: A combination of bioinformatics, comparative genomics and transcript profiling was used to identify and characterize 91 novel transcripts (NTs) represented in a cattle placenta cDNA library. These NTs have no significant similarity to any non-ferungulate DNA or RNA sequence. Proteins longer than 100 aa were predicted for 29 NTs, and 21 are candidate non-coding RNAs. Eighty-six NTs were found to be expressed in one or more of 18 different tissues, with 39 (42%) showing tissue-preference, including six that were expressed exclusively in placentome. The authenticity of the NTs was confirmed by their alignment to cattle genome sequence, 42 of which showed evidence of mRNA splicing. Analysis of the genomic context where NT genes reside revealed 61 to be in intergenic regions, whereas 30 are within introns of known genes. The genes encoding the NTs were found to be significantly associated with subtelomeric regions. CONCLUSION: The 91 lineage-specific transcripts are a useful resource for studying adaptive evolutionary responses of the ruminant placenta. The presence of so many genes encoding NTs in cattle but not primates or rodents suggests that gene loss and gain are important mechanisms of genome evolution in mammals. Furthermore, the clustering of NT genes within subtelomeric regions suggests that such regions are highly dynamic and may foster the birth of novel genes. The sequencing of additional vertebrate genomes with defined phylogenetic relationships will permit the search for lineage-specific genes to take on a more evolutionary context that is required to understand their origins and functions.

Dosage compensation is less effective in birds than in mammals.

BACKGROUND: In animals with heteromorphic sex chromosomes, dosage compensation of sex-chromosome genes is thought to be critical for species survival. Diverse molecular mechanisms have evolved to effectively balance the expressed dose of X-linked genes between XX and XY animals, and to balance expression of X and autosomal genes. Dosage compensation is not understood in birds, in which females (ZW) and males (ZZ) differ in the number of Z chromosomes. RESULTS: Using microarray analysis, we compared the male:female ratio of expression of sets of Z-linked and autosomal genes in two bird species, zebra finch and chicken, and in two mammalian species, mouse and human. Male:female ratios of expression were significantly higher for Z genes than for autosomal genes in several finch and chicken tissues. In contrast, in mouse and human the male:female ratio of expression of X-linked genes is quite similar to that of autosomal genes, indicating effective dosage compensation even in humans, in which a significant percentage of genes escape X-inactivation. CONCLUSION: Birds represent an unprecedented case in which genes on one sex chromosome are expressed on average at constitutively higher levels in one sex compared with the other. Sex-chromosome dosage compensation is surprisingly ineffective in birds, suggesting that some genomes can do without effective sex-specific sex-chromosome dosage compensation mechanisms.

Differential expression profiling of the blind subterranean mole rat Spalax ehrenbergi superspecies: bioprospecting for hypoxia tolerance.

The blind subterranean mole rat of the Spalax ehrenbergi superspecies, living underground and exposed to fluctuating oxygen and carbon dioxide levels, is an excellent model of hypoxic tolerance. Unique structural and functional adaptations of the cardiovascular and respiratory systems allow these underground mammals to survive at severely reduced oxygen tension. Elucidation of the natural variation and evolutionary changes under hypoxia within this superspecies may have biomedical applications in ischemic syndromes and cancer. In this study, we have compared expression profiles of muscle tissue at normoxic (21%) and hypoxic (3%) levels of oxygen concentration between two allospecies of the S. ehrenbergi superspecies exhibiting differential hypoxia tolerance in accordance with their ecological regimes. Profiling was performed by cross-species hybridization using a mouse cDNA array containing 15,000 gene elements. Results uncover species-specific responses to hypoxic stress among numerous genes involved in angiogenesis, apoptosis, and oxidative stress management. Among the most striking results are differential expressions of cardiac ankyrin repeat protein (Carp), activating transcription factor 3 (Atf3), LIM and cysteine-rich domains 1 (Lmcd1), cysteine and glycine-rich protein 2 (Csrp2), and ras homolog gene family, member B (RhoB). These findings support the hypothesis that allospecies of the S. ehrenbergi superspecies are variably adapted to fluctuating oxygen tension. Differences may involve specific metabolic pathways and functional adaptations at the structural and molecular levels.

Discovery of eight novel divergent homologs expressed in cattle placenta.

Ten divergent homologs were identified using a subtractive bioinformatic analysis of 12,614 cattle placenta expressed sequence tags followed by comparative, evolutionary, and gene expression studies. Among the 10 divergent homologs, 8 have not been identified previously. These were named as follows: cattle cerebrum and skeletal muscle-specific transcript 1 (CSSMST1), cattle intestine-specific transcript 1 (CIST1), hepatitis A virus cellular receptor 1 amino-terminal domain-containing protein (HAVCRNDP), prolactin-related proteins 8, 9, and 11 (PRP8, PRP9, and PRP11, respectively) and secreted and transmembrane protein 1A and 1B (SECTM1A and SECTM1B, respectively). In addition, two previously known divergent genes were identified, trophoblast Kunitz domain protein 1 (TKDP1) and a new splice variant of TKDP4. Nucleotide substitution analysis provided evidence for positive selection in members of the PRP gene family, SECTM1A and SECTM1B. Gene expression profiles, motif predictions, and annotations of homologous sequences indicate immunological and reproductive functions of the divergent homologs. The genes identified in this study are thus of evolutionary and physiological importance and may have a role in placental adaptations.

Accelerated expansion of group IID-like phospholipase A2 genes in Bos taurus.

Low-molecular-weight, calcium-dependent phospholipase A2 genes (PLA2s) that belong to the secreted type of PLA2s are clustered within a syntenic group on human 1p35-p36 and mouse 4qD3. We reassembled trace files available from the Whole Genome Sequencing (WGS) Project, obtaining an 86-kb contig with three tandem PLA2G2D duplications in the Hereford strain. We used mate-pair data to monitor the assembly and to exclude chimeric clones, demonstrating that the current WGS data may be assembled even in a highly repetitive region with a coverage exceeding fivefold. The genomic structure indicated that most of the PLA2G2D transcripts are formed by four exons. Two alternative first exons were present in all duplications. In two duplications insertions of satellite DNA in the third intron created a novel exon that gave rise to a two-exon product. Linkage and comparative mapping placed the bovine PLA2G2 locus on BTA2, indicating that it evolved from an ancestral PLA2G2D locus common to human, cattle, and rodents. Bovine PLA2G2D variants were capable of encoding 147-amino-acid polypeptides that consisted of putative signal peptide and metal-binding domains. Cysteine residues were conserved in positions analogous to those forming the seven disulfide bonds characteristic of PLA2G2 genes. Quantitative PCR analysis of bovine PLA2G2D transcripts indicated that their expression levels varied between the dry period and lactation in the mammary gland samples and that their expression was polymorphic in liver tissue. The recent burst of duplication and divergence of the bovine PLA2G2D genes and their polymorphic nature are typical of innate immune response genes.

Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle.

We previously localized a quantitative trait locus (QTL) on chromosome 6 affecting milk fat and protein concentration to a 4-cM confidence interval, centered on the microsatellite BM143. We characterized the genes and sequence variation in this region and identified common haplotypes spanning five polymorphic sites in the genes IBSP, SPP1, PKD2, and ABCG2 for two sires heterozygous for this QTL. Expression of SPP1 and ABCG2 in the bovine mammary gland increased from parturition through lactation. SPP1 and all the coding exons of ABCG2 and PKD2 were sequenced for these two sires. The single nucleotide change capable of encoding a substitution of tyrosine-581 to serine (Y581S) in the ABCG2 transporter was the only polymorphism corresponding to the segregation status of all 3 heterozygous and 15 homozygous sires for the QTL in the Israeli and U.S. Holstein populations. The allele substitution fixed effects on the genetic evaluations of 335 Israeli sires were -341 kg milk, +0.16% fat, and +0.13% protein (F-value = 200). No other polymorphism gave significant effect for fat and protein concentration in models that also included Y581S. The allele substitution effects on the genetic evaluations of 670 cows, daughters of two heterozygous sires, were -226 kg milk, 0.09% fat, and 0.08% protein (F-value = 394), with partial dominance towards the 581S homozygotes. We therefore propose that Y581S in ABCG2 is the causative site for this QTL.

A 7872 cDNA microarray and its use in bovine functional genomics.

The strategy used to create and annotate a 7872 cDNA microarray from cattle placenta and spleen cDNA sequences is described. This microarray contains approximately 6300 unique genes, as determined by BLASTN and TBLASTX similarity search against the human and mouse UniGene and draft human genome sequence databases (build 34). Sequences on the array were annotated with gene ontology (GO) terms, thereby facilitating data analysis and interpretation. A total of 3244 genes were annotated with GO terms. The array is rich in sequences encoding transcription factors, signal transducers and cell cycle regulators. Current research being conducted with this array is described, and an overview of planned improvements in our microarray platform for cattle functional genomics is presented.

Induction of overlapping genes by fasting and a peroxisome proliferator in pigs: evidence of functional PPARalpha in nonproliferating species.

Peroxisome proliferator-activated receptor alpha (PPARalpha), a key regulator of fatty acid oxidation, is essential for adaptation to fasting in rats and mice. However, physiological functions of PPARalpha in other species, including humans, are controversial. A group of PPARalpha ligands called peroxisome proliferators (PPs) causes peroxisome proliferation and hepatocarcinogenesis only in rats and mice. To elucidate the role of PPARalpha in adaptation to fasting in nonproliferating species, we compared gene expressions in pig liver from fasted and clofibric acid (a PP)-fed groups against a control diet-fed group. As in rats and mice, fasting induced genes involved with mitochondrial fatty acid oxidation and ketogenesis in pigs. Those genes were also induced by clofibric acid feeding, indicating that PPARalpha mediates the induction of these genes. In contrast to rats and mice, little or no induction of genes for peroxisomal or microsomal fatty acid oxidation was observed in clofibric acid-fed pigs. Histology showed no significant hyperplasia or hepatomegaly in the clofibric acid-fed pigs, whereas it showed a reduction of glycogen by clofibric acid, an effect of PPs also observed in rats. Copy number of PPARalpha mRNA was higher in pigs than in mice and rats, suggesting that peroxisomal proliferation and hyperresponse of several genes to PPs seen only in rats and mice are unrelated to the abundance of PPARalpha. In conclusion, PPARalpha is likely to play a central role in adaptation to fasting in pig liver as in rats and mice.

27-Hydroxycholesterol inhibits neutral sphingomyelinase in cultured human endothelial cells.

To study the effect of 27-hydroxycholesterol (27OHC) on the catabolism of sphingomyelin, we cultured endothelial cells (ECs) from human umbilical veins with 27OHC, then measured activities of acid sphingomyelinase (ASMase) and neutral sphingomyelinase (NSMase) and sphingomyelin consumption by using [14C]sphingomyelin, and determined NSMase mRNA expressions by RT-PCR method. The results indicated that [14C]sphingomyelin accumulated in cells treated with 27OHC, and that the activities of both NSMase and ASMase were inhibited in ECs cultured with 27OHC. To further study the effect of 27OHC on NSMase, we used desipramine, an inhibitor of ASMase, to exclude the possible interference of ASMase's residual activity at neutral condition. Also, we observed the significant inhibition of NSMase activity by using glutathione, an inhibitor of NSMase, but found no further impact when 27OHC was added later. To determine whether the inhibition of NSMase activity was directly due to the effect of 27OHC, we exposed cell homogenate to 27OHC, and found no inhibitive effect of 27OHC on the activity of NSMase. All of our data confirmed that 27OHC had only an indirect inhibitive effect on NSMase. Our finding that no change of the NSMase mRNA expression by 27OHC indicated that the inhibitive effect of 27OHC on NSMase activity occurred at a post-transcriptional level. We suggest that an altered membrane fluidity caused by 27OHC could be involved in the inhibited activity of NSMase.

A 1463 gene cattle-human comparative map with anchor points defined by human genome sequence coordinates.

A second-generation 5000 rad radiation hybrid (RH) map of the cattle genome was constructed primarily using cattle ESTs that were targeted to gaps in the existing cattle-human comparative map, as well as to sparsely populated map intervals. A total of 870 targeted markers were added, bringing the number of markers mapped on the RH(5000) panel to 1913. Of these, 1463 have significant BLASTN hits (E < e(-5)) against the human genome sequence. A cattle-human comparative map was created using human genome sequence coordinates of the paired orthologs. One-hundred and ninety-five conserved segments (defined by two or more genes) were identified between the cattle and human genomes, of which 31 are newly discovered and 34 were extended singletons on the first-generation map. The new map represents an improvement of 20% genome-wide comparative coverage compared with the first-generation map. Analysis of gene content within human genome regions where there are gaps in the comparative map revealed gaps with both significantly greater and significantly lower gene content. The new, more detailed cattle-human comparative map provides an improved resource for the analysis of mammalian chromosome evolution, the identification of candidate genes for economically important traits, and for proper alignment of sequence contigs on cattle chromosomes.

A 3800 gene microarray for cattle functional genomics: comparison of gene expression in spleen, placenta, and brain.

A cDNA microarray representing approximately 3800 cattle genes was created for functional genomic studies. The array elements were selected from > 7000 cDNA clones identified in a large-scale expressed sequence tag (EST) project that utilized spleen and normalized and subtracted placenta cDNA libraries. Sequence similarity searches of the 3820 ESTs represented on the array using BLASTN identified 3290 (86.1%) as putative human orthologs, with the remainder consisting of "novel" genes or highly divergent orthologs. Experiments were conducted with a prototype 768 gene microarray created from spleen cDNAs and with the 3800 gene array that included genes from spleen and placenta. The 768 gene array was used to profile RNA transcripts expressed by adult and fetal spleen. The 3800 gene array was used to profile transcripts expressed by adult brain and placenta. Microarray analysis of RNA extracted from fetal and adult spleen identified 29 genes that were differentially expressed two-fold or more. Transcriptional differences of two of these genes, IGJ and CTSS, were confirmed using TaqMan technology. The comparison of brain and placenta revealed 400 genes expressed at higher levels in brain and 72 genes expressed at higher levels in placenta. These results demonstrate the potential power of microarrays for understanding the molecular mechanisms of cattle development, disease resistance, nutrition, fertility and production traits.

Annotated expressed sequence tags and cDNA microarrays for studies of brain and behavior in the honey bee.

To accelerate the molecular analysis of behavior in the honey bee (Apis mellifera), we created expressed sequence tag (EST) and cDNA microarray resources for the bee brain. Over 20,000 cDNA clones were partially sequenced from a normalized (and subsequently subtracted) library generated from adult A. mellifera brains. These sequences were processed to identify 15,311 high-quality ESTs representing 8912 putative transcripts. Putative transcripts were functionally annotated (using the Gene Ontology classification system) based on matching gene sequences in Drosophila melanogaster. The brain ESTs represent a broad range of molecular functions and biological processes, with neurobiological classifications particularly well represented. Roughly half of Drosophila genes currently implicated in synaptic transmission and/or behavior are represented in the Apis EST set. Of Apis sequences with open reading frames of at least 450 bp, 24% are highly diverged with no matches to known protein sequences. Additionally, over 100 Apis transcript sequences conserved with other organisms appear to have been lost from the Drosophila genome. DNA microarrays were fabricated with over 7000 EST cDNA clones putatively representing different transcripts. Using probe derived from single bee brain mRNA, microarrays detected gene expression for 90% of Apis cDNAs two standard deviations greater than exogenous control cDNAs. [The sequence data described in this paper have been submitted to Genbank data library under accession nos. BI502708-BI517278. The sequences are also available at http://titan.biotec.uiuc.edu/bee/honeybee_project.htm.]