- Research article
- Open Access
Sex-dependent gene expression in early brain development of chicken embryos
BMC Neuroscience volume 7, Article number: 12 (2006)
Differentiation of the brain during development leads to sexually dimorphic adult reproductive behavior and other neural sex dimorphisms. Genetic mechanisms independent of steroid hormones produced by the gonads have recently been suggested to partly explain these dimorphisms.
Using cDNA microarrays and real-time PCR we found gene expression differences between the male and female embryonic brain (or whole head) that may be independent of morphological differentiation of the gonads. Genes located on the sex chromosomes (ZZ in males and ZW in females) were common among the differentially expressed genes, several of which (WPKCI-8, HINT, MHM non-coding RNA) have previously been implicated in avian sex determination. A majority of the identified genes were more highly expressed in males. Three of these genes (CDK7, CCNH and BTF2-P44) encode subunits of the transcription factor IIH complex, indicating a role for this complex in neuronal differentiation.
In conclusion, this study provides novel insights into sexually dimorphic gene expression in the embryonic chicken brain and its possible involvement in sex differentiation of the nervous system in birds.
Sexual differences that arise in the brain during embryonic development underlie sex-specific reproductive behavior in adults of vertebrate species. Better knowledge of the mechanisms behind sexual differentiation of the nervous system can contribute to the understanding of the brain's functions and its susceptibility to disease , and is critically needed in the context of developmental neurotoxicity and chemically induced disruption of the neuroendocrine system . For a long time the idea prevailed that somatic tissues in the embryo are gender neutral until the gonads become hormonally active following their sex-specific differentiation, a process presumed to be governed by cell-autonomous programs genetically determined by the sex chromosomes [3–5]. However, several recent studies indicate that sex-specific neuronal characteristics may be influenced by genetic mechanisms preceding or occurring in parallel with the effects produced by the gonadal hormones . That the genetic sex of neurons contributes to the process of sexual differentiation is illustrated by studies of the neural song circuit in a gynandromorphic zebra finch, where genetically male and female brain cells develop differently in the presence of the same levels of circulating gonadal hormones .
Birds, in which males are homogametic (ZZ) and females heterogametic (ZW), provide attractive models for the study of normal and disturbed sex-related neurobehavioral development [7, 8]. The recent advent of chicken expressed sequence tag (EST) clone collections [9, 10] and the release of the first draft sequence of the chicken genome  has paved the way for genome-wide studies of sexual differentiation in the chicken. In the chicken, the indifferent embryonic gonad rudiments differentiate morphologically around embryonic day (ed) 6.5 [12, 13], whereas the expression of enzymes involved in steroidogenesis can be detected as early as ed2 , and estrogen receptor (ER) mRNA is detectable in the male and female urogenital system from ed4.5 . The earliest signs of sexually dimorphic expression of genes related to gonadal hormones and steroidogenesis, however, have been detected at ed5 when mRNA for anti-Müllerian-hormone is more abundant in males , and aromatase  in females.
In the present study, we used cDNA microarrays to investigate if there are gene expression differences between the male and female embryonic brain (or whole head) that may be independent of morphological differentiation of the gonads. We report evidence for sex-related gene expression differences from ed4 (whole head) through ed10 (brain) for several genes located on sex chromosomes. We propose that intrinsic genetic mechanisms are involved in sexual differentiation of the chicken brain, and that sex-linked genes may play key roles.
Identification of gene expression differences between male and female embryonic chicken brains
To study sex-specific gene expression in the developing brain of chicken embryos during stages around morphological differentiation of the gonads, total RNA extracted from male and female heads on ed4 and ed6, and brains on ed8 and ed10, was subjected to replicated microarray analysis using cDNA microarrays containing approximately 14,000 EST clones. For each sex, we generated a total of 12 biologically independent samples, comprising either individuals (ed10) or pools of four individuals (ed4, 6, 8 and 10). These 24 samples were analyzed on 20 array slides, with technical replication as shown in Figure 1. Using an empirical Bayes moderated t-test , and adjusted p-values according to Benjamini and Hochberg , we found 146 clones with significant differential expression (adjusted p-value < 0.0005, corresponding to a B score threshold of 3.413) across all investigated embryonic stages (Figure 2). For most of these clones, the log2 female/male ratios were relatively similar across the studied developmental stages (Figure 3, see Additional File 1). Importantly, for a majority of these clones the expression differences between males and females were present at the earliest embryonic stage (Figure 3, see Additional File 1), i.e. well before morphological gonadal differentiation.
Characterization of differentially expressed genes
The 146 identified clones (see Additional File 2) were found to represent 133 genes, 19 of which were more highly expressed in females, and 114 in males (Table 1). Among the 'female genes', three are located on the female-specific W chromosome (Table 2), one on the Z chromosome (Table 2) and 11 on autosomes, whereas for four genes the location is unknown (Table 1). Among the 114 'male genes', 61 are located on the Z chromosome (Table 2) and 29 on autosomes, and for another 29 genes the location is unknown (Table 1). Hence, given the fraction of Z-linked clones present in the microarray (257 clones out of ~14,000), genes annotated as being located on the Z chromosome appeared overrepresented among the differentially expressed genes (Table 1). As the chicken W chromosome is still poorly characterized, with many gaps in the sequence, the same conclusion is less straightforward to draw about W-linked genes (for further details see the UCSC Genome Browser Gateway ).
The differential expression of selected genes on ed4 (and, for some genes, ed8) was investigated by real-time PCR of individual samples (see Methods) and the results (see Additional File 3) are presented below in the context of the chromosomal locations of these genes.
Since the W chromosome is present only in female chickens, the expression of W-linked genes was reflected in the microarray analysis by high log2 female/male ratios (Figure 3), and in the real-time PCR analysis by even higher (technically 'infinite') ratios (Table 3). An exception from this was SPINW, for which the exclusive expression in females was not indicated by similarly high log2 female/male ratios in the microarray analysis (see Additional File 2). However, unlike what is known for the other W-linked genes reported here, SPINW has a counterpart (SPINZ) on the Z chromosome (see below). SPINW was found by real-time PCR to be expressed in the female head on ed4 (Table 3, see Additional File 3). Also the expression of WPKCI-8, for which female expression at all stages was detected by seven clones in the microarray analysis (see Additional File 2), was verified in the ed4 female head by real-time PCR.
The sequence of the EST clone WLA084D05 shows no similarity to any known gene or putative protein. It is part of the UniGene entry Gga.16155 and the TIGR Gene Index entry TC216880 (Release 10.0), located about 300 nucleotides from a chicken sequence similar (>90% sequence alignment match at the amino acid sequence level) to the Ube2r2 mouse sequence (see Additional Files 4 and 5). Using primers designed to amplify this putative Ube2r2 homolog, we did not detect any significant differential expression in embryonic chicken heads. However, using primers amplifying the WLA084D05 sequence, we found expression in heads from 4-day-old female embryos (Table 3). It remains to determine whether this EST may represent a novel W-linked gene (Figure 4, see Additional File 4), which we provisionally suggest calling 'avian brain W-linked transcript' (ABWT).
Whereas half of the genes (56 of 114) suggested from the microarray study to be more highly expressed in males are located on the Z chromosome, only one clone with female-enhanced expression is known to be Z-linked (Tables 1 and 2). We found 141 nucleotides of this clone (RJA001B07) to align with high similarity (98% nucleotide sequence identity) to a partial sequence of the GenBank entry with accession AB046698 (Figure 5), representing a Z-linked male hypermethylated (MHM) region. This sequence is thought to be specifically transcribed into a non-coding nuclear RNA in females . Using primers designed to discriminate between the RJA001B07 and AB046698 sequences (Figure 5), we found female-specific transcription for both of these MHM genomic regions (Table 4).
Real-time PCR analysis confirmed the male-enhanced microarray expression for all Z-linked genes examined (Table 4). We also found significantly higher expression in males (p < 0.05) for two Z-linked genes (CDK7 and the Z-linked WPKCI-8 homolog HINT), but not for a third (DMRT1), that were not represented in the microarray but are functionally related (see Discussion) to other genes identified as differentially expressed (Table 4).
A close to two-fold higher expression in males was observed for most of the Z-linked genes in the real-time PCR analysis (Table 4), suggestive of gene dosage. Exceptions were BTF2-P44, CDK7 and HMGCOA. No differential expression was detected by real-time PCR for three genes (ANXA1, ZO-2 and ZOV3) that were randomly selected among the non-differentially expressed Z-linked genes in the microarray experiment (Table 4), supporting that only certain Z-linked genes are more highly expressed in males. A fourth Z-linked gene (TXN) randomly selected among genes that were not differentially expressed in the microarray analysis was revealed by real-time PCR to be more highly expressed in males (Table 4). However, examination of the corresponding microarray clone (RJA044E12) by sequencing showed that the spotted cDNA contained vector sequence only.
Genes located on autosomes or having unknown chromosomal location
Differential expression could be confirmed by real-time PCR for only three genes out of 14 analyzed genes with autosomal (P311) or unknown location (ACAA2 and CETN3), all with higher expression in males (Table 5). Errors in the microarray were uncovered for 6 of 11 genes that could not be confirmed by real-time PCR (Tables 4 and 5). For example, the microarray spots thought to represent the male-enhanced autosomal genes AATF and ACAD8 were found to contain EST clones representing the Z-linked genes RPS6 and PLAA, respectively (Table 5). Two of the microarray clones showing female-enhanced expression that could not be confirmed by real-time PCR were found to be erroneous and contain repetitive sequences (data not shown), whereas two such clones (representing the autosomal genes CEZANNE2 and NAT5) appeared to be correctly annotated (Table 5). The discrepancy between microarray and real-time PCR data for non-sex chromosomal gene expression remains unexplained. Unspecific cross-hybridization of sex chromosome-linked target cDNA to non-sex chromosomal array probes may be one possibility. Because of this uncertainty we decided not to further follow up on non-sex chromosomal candidates. Only microarray data on sex-chromosomal genes are reported in Table 2 (for the whole data set of 146 clones, see Additional File 2). We confirmed the absence of differential expression for the autosomal gene MAT1, a gene selected for PCR analysis because of its functional association with other genes identified as differentially expressed (see Discussion).
Relatively little is known about the genes that are involved in neuronal sex differentiation in birds, whether influenced by gonadal steroid hormones or not. In the present study we used cDNA microarrays to identify genes that are differentially expressed between males and females in the developing chicken brain (or whole head) before (ed4) and during (ed6, 8 and 10) morphological differentiation of the gonads. We found statistical significance (adjusted p-value < 0.0005) for differential expression of 146 microarray clones across the embryonic stages and types of tissue sample (whole head for ed4 and 6, and brain for ed8 and 10). These clones correspond to ~1.4% of the ~10,700 clones included in the analysis. Using real-time PCR we were able to confirm female expression for all three W-linked genes identified, and the differential expression of all ten Z-linked genes tested and three genes with autosomal or unknown location. The real-time PCR analysis furthermore revealed one false negative among Z-linked genes (TXN), and some false positives among genes located on autosomal or unknown chromosomes. In many cases, however, this could be traced to errors in the microarrays.
A major finding was that sex-specific differential expression could be detected for numerous genes already at the earliest studied embryonic stage (ed4), before the stage (~ed5) at which the earliest sex differences in gonadal gene expression with regard to steroidogenesis have been found [15, 16, 21]. Expression of several genes for steroidogenic enzymes (P450scc, P450c17, 3-βHSD, 17βHSD, and aromatase) in chicken gonads already on ed4 has been documented, but no sex-dependent expression differences . Even if we cannot rule out the possibility that hormonal differences between males and females may exist on ed4, the early sex-specific expression of the genes reported in the present study may be primarily regulated by other mechanisms than hormonal control. In favor of this notion, we found no apparent relationship between the expected differences in hormonal levels associated with gonadal differentiation  and the magnitude of the expression differences between the sexes for these genes (Figure 3). Moreover, our preliminary data from chicken embryos exposed to ethinyl estradiol in ovo on ed4 show that the expression in the brain on ed10 for the majority of these genes is virtually unaffected by treatment in either sex (work in progress).
Gene expression differences independent of gonadal hormones have previously been suggested to play critical roles in the sexual differentiation of the brain in birds  and mice . Why an apparently different set of genes should be involved in the chicken compared to the mouse  is unclear. That such a large proportion of genes located on the sex chromosomes (Z and W) were identified as differentially expressed in the chicken is therefore intriguing considering the independent evolution and different gene content of the sex chromosomes in birds and mammals .
A majority of the genes indicated to be differentially expressed from the microarray analysis were more highly expressed in males than in females, and about half of these genes are located on the Z chromosome. This raises the issue of dosage compensation versus biallelic expression of Z-linked genes in males . Although biallelic expression has previously been demonstrated [26, 27], and might be inferred from the close to two-fold higher expression in males of many Z-linked genes in our study, the expression levels of some Z-linked genes were found both by microarray and PCR analysis not to differ significantly between the sexes. Such a lack of sex differences in the expression of certain genes located on the Z chromosome would be in line with previous evidence for dosage compensation . Taken together, our results suggest that if dosage compensation does occur in chicken it is unlikely to involve a widespread inactivation of the Z chromosome similar to that of the mammalian X chromosome. Studies by McQueen and coworkers  in chicken embryos have shown that at least six genes on the Z chromosome are dosage compensated. However, Kuroda and coworkers  found that transcription of five genes, including two of the genes studied by McQueen and coworkers , is taking place on both Z chromosomes of male chicken. This provides further support that dosage compensation in birds does not involve inactivation of a large majority of genes on the Z chromosome.
That epigenetic mechanisms such as methylation and acetylation may be involved in the sex-specific expression of genes located on the Z chromosome has previously been implicated by the female-specific expression of the Z-linked MHM region [20, 29]. In agreement with this, our microarray results showed high female/male ratios suggestive of expression in females only of two distinctive transcripts containing an MHM region. Whether these transcripts, represented by microarray clone RJA001B0 (Table 2) and GenBank entry AB046698, respectively, are expressed from the same or independent genomic regions is unclear. Their highly similar MHM regions (141 nucleotides with 98% identity) and their female-specific expression indicate some common function, possibly the repression of the adjacent gene DMRT1, suggested as a conserved sex determining gene [30–35], through the accumulation of a non-coding RNA . However, unlike in male gonads before and during gonadal differentiation [32, 36, 37], we found no evidence of male-enhanced expression of DMRT1 in the embryonic chicken brain.
We found high female/male expression ratios on ed4 for three genes/clones (WPKCI-8, SPINW and ABWT) located on the female-specific W chromosome, indicating significant early expression in the female embryonic brain. Because of its early expression (ed4.5) in the developing female gonads [38, 39], WPKCI-8 (also known as ASW) has been suggested to be involved in avian sex determination . Early WPKCI-8 expression (ed5) has also been detected in the spinal cord, spinal ganglion, and myotomes . The WPKCI-8 gene, which is reiterated approximately 40 times on the W chromosome , and its single-copy homolog on the Z chromosome HINT  both encode proteins belonging to the Hint family of dimeric nucleotide hydrolases . The Wpkci-8 protein, however, lacks the histidine triad (HIT) motif  and appears to act as a dominant suppressor of Hint activity, possibly by conferring signals for mislocalization or degradation of the Hint/Wpkci-8 heterodimers  (Figure 6). It is an interesting possibility that suppressed levels of active Hint in the brain of female chicken embryos, resulting from both lower expression of the HINT gene itself and inhibition by expressed Wpkci-8, could contribute to sexual differentiation of the nervous system. Female expression of WPKCI-8 and male-enhanced expression of HINT has previously been detected in the telencephalon of juvenile zebra finches . The Hint protein (also known as protein kinase C inhibitor/interacting protein) is well conserved between species [38, 39], but whether this protein is involved in sex determination or sexual differentiation in various species is unclear.
We found male-enhanced expression for two Z-linked genes (CCNH and CDK7) known to encode subunits of the cyclin-dependent kinase activating kinase (CAK) complex, and a third one (BTF2-P44) encoding a subunit of the transcription factor II H (TFIIH) core complex. Together these complexes constitute the basal transcription complex TFIIH, which is involved in transcription, DNA repair, and cell-cycle control . Among the genes encoding the nine proteins (CDK7, CCNH, MAT1, ERCC2, ERCC3, BTF2-P34, BTF2-44, BTF2-52, BTF2-62) in this complex , two (ERCC2 and BTF2-52) appear to be less conserved in the chicken genome (see Additional Files 6 and 7). The implications of the male-enhanced expression of some but not all components of CAK (we found no differential expression for its third component MAT1) and the TFIIH core complex remains to be elucidated. However, it is intriguing that both CAK and TFIIH are involved in estrogen receptor α (ERα) and androgen receptor (AR) transactivation in mammals [42–45], and that a functional CAK is important for both the mitotic and meiotic cell cycle [46–48]. The proteins encoded by the Z-linked genes HINT and APTX, are both members of the HIT protein superfamily . APTX, which in the microarrays is more highly expressed in males, encodes aprataxin, a protein involved in single-strand break repair [50, 51] and Hint seems to function as a positive regulator of components of TFIIH , including the CCNH homolog Ccl1 . Considering the role of TFIIH in DNA-repair, a sex specific control of DNA repair seems possible. In addition to TFIIH components, three other Z-linked genes associated with DNA repair (RAD1 , MSH3 , and DHFR ) were found in the microarray result to be more highly expressed in males.
Analogous to the sex-dependent expression of WPKCI-8 and HINT, the W-linked gene SPINW was expressed in females while the Z-linked homolog SPINZ showed male-enhanced expression. SPINW, which is prominently expressed in ovarian and somatic tissues, has been suggested to be involved in gonadal differentiation , but the early (ed4) sex-specific expression of these genes in the brain may indicate a more widespread involvement in sexual differentiation. We also identified a W-linked transcript, ABWT, which in addition to being expressed in the early female embryonic brain has also been found in cDNA libraries from chicken ovaries and heart (UniGene entry Gga.16155). ABWT maps to the genome in the vicinity of the homolog to mouse Ube2r2 (i.e. chicken gene UBE2R2), a gene belonging to the E2 family of ubiquitin-conjugating enzymes . As we were unable to detect any differential expression of UBE2R2 transcripts, however, the detected ABWT transcripts appear to be distinct from these. It remains to be determined whether ABWT is a splice variant of a UBE2R2 transcript, a novel gene encoding a yet unidentified protein or a non-coding RNA.
In conclusion, this study provides hints to mechanisms behind sexual differentiation of the nervous system in birds and raises several new questions. The expression of W-linked genes supports the presence of 'sex-control genes' similar to the mammalian SRY gene, in addition to the possible gene dosage effects from genes on the Z chromosomes. Future studies will reveal if the W- and Z-linked genes are differentially expressed in the early embryonic brain of different bird species and if differential expression of these genes is present also in other tissues. It should be noted, however, that WPKCI-8, an apparent key player in sexual differentiation [13, 38], is not reiterated on , and may even not be linked to , the W chromosome of the primitive ratites (emus and ostriches). The crucial importance of sex steroids for organization of the brain and behavior in birds has been shown in several studies (as reviewed by Balthazart et al  and Panzica et al ), but the results presented in this paper point toward the possible presence also of a genetic component in sex-specific neuronal differentiation. Consequently, it seems likely that genetic and hormonal control interact to organize the avian brain dimorphism during differentiation.
Embryos and sample collection
Fertilized eggs from White Leghorn fowl were purchased from OVA Production (Morgongåva, Sweden). The eggs were incubated at 37.5°C and 60% relative humidity, and were turned every 3 hours. After 8 or 10 days of incubation (ed8 and ed10), whole brains were dissected out. As brains could not be excised in a reproducible way from earlier embryos (ed4 and ed6), whole heads were collected. Heads were excised along a line from the mandibular and maxillary processes of the first pharyngeal arch to the caudal boundary of the myelencephalon (ed4), or from immediately in front of the maxillary process to the caudal boundary of the myelencephalon (ed6). Samples were immediately frozen in liquid nitrogen, and stored at -70°C. A tissue sample from each embryo was also collected for DNA isolation and genetic sexing according to a PCR-based method  in which intron sequences of different lengths in the W-linked gene CHD1W (females) and Z-linked gene CHD1Z (both sexes) are amplified.
Total RNA was isolated using the Micro-to-Midi Total RNA Purification System (Invitrogen, Carlsbad, CA). Any contaminating genomic DNA was digested by DNase treatment (DNA-free, Ambion, Austin, USA), according to the manufacturer's recommendations. RNA quality was checked using the Agilent 2100 Bioanalyzer and the RNA 6000 LabChip (Agilent Technologies, Palo Alto, CA, USA). Only high quality RNA, with no signs of degradation, was used for further experiments.
Spotted cDNA microarrays, containing 1,136 expressed-sequence-tag (EST) clones from a cDNA library of developing chicken brain  and 12,771 EST clones from four cDNA libraries from brain and testis of White Leghorn and Red Jungle Fowl , spotted in duplicate, were purchased from the Royal Institute of Technology, Stockholm, Sweden . Annotation and chromosome localization of the spotted clones can be found at the Stockholm Bioinformatics Center (SBC) .
Four microarray hybridizations were done for each embryonic stage (ed4, 6, 8 and 10), addressing biological and technical variation through a pooling and dye reversal strategy (Figure 1A). Equal amounts of total RNA extracted from eight male and eight female individual embryos were mixed in four pools with RNA from 4 embryos of the same sex in each pool.
This design, which we adopted from Churchill , is a trade-off between the need for biological and technical replication while keeping down animal consumption and microarray use. Although pooling is intended to average between individuals there is a risk that single individuals may bias the sample. Individual variation among the embryos included in the pools was addressed in the PCR analysis (see below), and was found to be relatively modest for the randomly selected individuals assayed (see Additional File 3). Moreover, four additional microarray hybridizations were done for ed10, using an alternative design (Figure 1B) in which all samples were derived from individual embryos and using other embryos than in the pooling strategy.
cDNA synthesis and microarray hybridization
Labeling of cDNA for microarray hybridization was done using the 3DNA Array 350 Expression Array Detection Kit (Genisphere Inc., Hatfield, PA) according to the manufacturer's protocol and recommendations. Briefly, total RNA (20 μg) was oligo-dT-primed with Cy3- or Cy5-capture-sequence primer (see Figure 1 for dye reversal scheme), and reverse transcribed using SuperScript II Reverse Transcriptase (Invitrogen, Carlsbad, CA). For cDNA hybridization, 5 μg cDNA was mixed with 1 μl mouse COT-1 DNA and either regular (ed8 and ed10S) or Enhanced (ed4, ed6, ed10) hybridization buffer. After hybridization at 60°C for 16 hours, microarray slides were washed on a rocking platform, first with 2 × SSC containing 0.2% SDS for 10 min at 55°C, followed by 2 × SSC and 0.2 × SSC for 10 min each at room temperature. 3DNA hybridization was performed at 60°C for 2.5 hours, followed by washes as for the cDNA hybridization.
Microarray data analysis
The microarrays were scanned with a GenePix 4000B scanner (Axon, Foster City, CA) at 10 μm resolution. The photomultiplier tube voltage settings were varied to obtain maximum signal intensities while saturating less than 0.1% of the spots. Images were analyzed with the GenePix Pro 5.0 (Axon) software, utilizing the option to find irregular features. Spots with visible artifacts or containing fewer than 35 pixels were manually flagged as bad. As each array contained the same clone spotted in duplicate in each half of the array, the average of the two duplicate spots was calculated for each array. If one of the duplicate spots was missing or flagged bad, the value was based only on the remaining spot. The base 2 logarithm (log2) ratio of the median spot intensity for each channel was used to quantify the fold difference in relative gene expression levels. Since all arrays had low background, no background subtraction was done. To remove systematic sources of variation, within-print-group loess normalization  was done, in which the relative weight of 0.1 was given to spots flagged as missing or bad. To identify clones differentially expressed between females and males independently of embryonic stage (ed4, 6, 8 or 10), tissue origin (brain or whole head), or experimental design (pools or individuals; see above), the arrays were collapsed across these parameters, resulting in 12 biologically independent samples per gender analyzed on 20 arrays.
Clones for which values were obtained for less than 17 (out of 20) arrays were excluded from further analysis, resulting in 10,792 clones across all arrays. In addition, 90 of the microarray clones known from previous RepeatMasker analysis to mainly consist of repetitive sequence were excluded from further analysis. The moderated t-statistics and log-odds of differential expression was calculated using the Limma package version 2.0.4  in the statistical software R version 2.2.1  (both freely available). To address the problem with multiple testing, the p-values were adjusted using the method of Benjamini and Hochberg . Because of the technical replication of the pooled samples (Figure 1A), only 12 arrays (not all 20) were independent measurements. We only considered clones to be differentially expressed if the adjusted p-value < 0.0005. Microarray data are available in the ArrayExpress database  at the accession number E-MEXP-266.
Annotation of EST clones detecting differential expression
We reevaluated gene identities for all differentially expressed EST clones (146 clones; see Results) by searching The Institute of Genomic Research (TIGR) Chicken Gene Indices (Release 10.0) . All sequences can be accessed by using the clone identity in a search against GenBank database . Functional annotations were reassessed by RepeatMasker  and searching the National Center for Biotechnology Information (NCBI) non-redundant and the Universal Protein Resource (UniProt) databases [70, 71]. The RepeatMasker program screens nucleotide sequences against a library of repetitive elements such as interspersed repeats and low complexity DNA sequences. We discarded sequences with more than 10% repetitive sequence.
Chromosomal localizations in the released chicken genome  were reinvestigated using the BLAST-like Alignment Tool (BLAT)  and the University of California Santa Clara (UCSC) Genome Browser . Designating functionality to chicken gene products is complicated by the limited knowledge of the avian genome or proteome. Assuming that well conserved protein regions and domains reflect similar protein functionality, cross-species comparisons of several gene products (see Additional Files 5, 6, 7) were done using data from the USCS genome browser, the Ensembl genome browser , TIGR chicken gene indices , the Pfam  and UniProt databases in combination with the BLAST2 (default parameters) . The figures in Additional Files 4 and 5 show sequence similarity and Pfam domains for several of the gene products mentioned in the discussion. Pfam is a database of multiple alignments of protein domains or conserved protein regions . The thicker boxes represent Pfam A domains which are based on manually crafted accurate multiple protein alignments, whereas the multiple colored thinner boxes represent Pfam B domains based on an automatic clustering of a nonredundant protein database.
To confirm the identity as well as the differential expression of selected genes, samples derived from five individual embryos of each sex (randomly selected among the eight that had been pooled for microarray analysis) on ed4 (and, for some genes, ed8) were subjected to real-time PCR. Primers were designed with Primer Express software (Applied Biosystems, Palo Alto, CA, USA), using default settings for the TaqMan mode, and ordered from DNA Technology A/S (Aarhus, Denmark). Primers were designed for amplicons to span exon-exon borders when possible, as determined by the alignment of EST sequences against the chicken genome (see above). Primer sequences are given in Table 6. Primers amplifying the endogenous reference 18S ribosomal RNA were from TaqMan® Ribosomal RNA Control Reagents (Applied Biosystems). For PCR, 2 μg total RNA was reverse transcribed in a final volume of 100 μl using TaqMan Reverse Transcription Reagents (Applied Biosystems) with random hexamer primers according to manufacturer instructions. Reactions excluding MultiScribe Reverse Transcriptase (Applied Biosystems) were performed as negative controls. cDNA targets at a 100-fold (for target genes) or 10,000-fold (for endogenous reference) final dilution were amplified in replicate wells (four for target genes and six for endogenous reference), using primer concentration 250 nM (reverse and forward primer) for target genes and 50 nM for 18S, respectively, in 1× qPCR Mastermix Plus for SYBR Green I (Eurogentec, Seraing, Belgium) in an ABI Prism 7000 Sequence Detector System (Applied Biosystems) with the following thermal profile: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 15 sec at 95°C and 1 min at 60°C in a volume of 25 μl. PCR products were checked by monitoring melting curves. Standard curves for each gene were obtained by amplifying (in quadruplicate) a four-fold dilution series of 1:50 through 1:12,800 final volume of a reference mixture containing equal amounts of cDNA from all ten individual samples. For each gene, after removing one outlying value (leaving three values for target genes and five for endogenous reference), a mean normalized gene expression (MNE) was calculated according to Muller et al :
Etarget and Eref denote the efficiency, and CTtarget and CTref the threshold cycle (CT) of the target and reference gene, respectively, in the PCR amplification. PCR data for each gene are reported in Tables 3, 4, 5 (also see Additional File 3). An un-paired Student's t-test, using the MNE values, was used to determine a statistically significant difference between the genders (p < 0.05) for each gene. The comparative Ct method  was used to calculate the average relative log2 fold difference between males and females for each gene.
Avian brain W-linked transcript
BLAST-like alignment tool
ensemble genome browser
Expressed sequence tag
- Female genes:
genes more highly expressed in females
Cyclin-dependent kinase activating kinase
histidine triad motif
- Male genes:
genes more highly expressed in males
Mean normalized expression
Transcription factor II H
University of California Santa Cruz
Morris JA, Jordan CL, Breedlove SM: Sexual differentiation of the vertebrate nervous system. Nat Neurosci. 2004, 7: 1034-1039. 10.1038/nn1325.
Schantz SL, Widholm JJ: Cognitive effects of endocrine-disrupting chemicals in animals. Environ Health Perspect. 2001, 109: 1197-1206.
Arnold AP, Gorski RA: Gonadal steroid induction of structural sex differences in the central nervous system. Annu Rev Neurosci. 1984, 7: 413-442. 10.1146/annurev.ne.07.030184.002213.
MacLusky NJ, Naftolin F: Sexual differentiation of the central nervous system. Science. 1981, 211: 1294-1302.
Dennis C: Brain development: the most important sexual organ. Nature. 2004, 427: 390-392. 10.1038/427390a.
Arnold AP: Sex chromosomes and brain gender. Nat Rev Neurosci. 2004, 5: 701-708. 10.1038/nrn1494.
Agate RJ, Grisham W, Wade J, Mann S, Wingfield J, Schanen C, Palotie A, Arnold AP: Neural, not gonadal, origin of brain sex differences in a gynandromorphic finch. Proc Natl Acad Sci U S A. 2003, 100: 4873-4878. 10.1073/pnas.0636925100.
Halldin K, Axelsson J, Brunstrom B: Effects of endocrine modulators on sexual differentiation and reproductive function in male Japanese quail. Brain Res Bull. 2005, 65: 211-218. 10.1016/j.brainresbull.2004.11.020.
Boardman PE, Sanz-Ezquerro J, Overton IM, Burt DW, Bosch E, Fong WT, Tickle C, Brown WR, Wilson SA, Hubbard SJ: A comprehensive collection of chicken cDNAs. Curr Biol. 2002, 12: 1965-1969. 10.1016/S0960-9822(02)01296-4.
Savolainen P, Fitzsimmons C, Arvestad L, Andersson L, Lundeberg J: ESTs from brain and testis of White Leghorn and red junglefowl: annotation, bioinformatic classification of unknown transcripts and analysis of expression levels. To appear in Cytogenetic and Genome Research.
Hillier LW, Miller W, Birney E, Warren W, Hardison RC, Ponting CP, Bork P, Burt DW, Groenen MA, Delany ME, Dodgson JB, Chinwalla AT, Cliften PF, Clifton SW, Delehaunty KD, Fronick C, Fulton RS, Graves TA, Kremitzki C, Layman D, Magrini V, McPherson JD, Miner TL, Minx P, Nash WE, Nhan MN, Nelson JO, Oddy LG, Pohl CS, Randall-Maher J, Smith SM, Wallis JW, Yang SP, Romanov MN, Rondelli CM, Paton B, Smith J, Morrice D, Daniels L, Tempest HG, Robertson L, Masabanda JS, Griffin DK, Vignal A, Fillon V, Jacobbson L, Kerje S, Andersson L, Crooijmans RP, Aerts J, van der Poel JJ, Ellegren H, Caldwell RB, Hubbard SJ, Grafham DV, Kierzek AM, McLaren SR, Overton IM, Arakawa H, Beattie KJ, Bezzubov Y, Boardman PE, Bonfield JK, Croning MD, Davies RM, Francis MD, Humphray SJ, Scott CE, Taylor RG, Tickle C, Brown WR, Rogers J, Buerstedde JM, Wilson SA, Stubbs L, Ovcharenko I, Gordon L, Lucas S, Miller MM, Inoko H, Shiina T, Kaufman J, Salomonsen J, Skjoedt K, Wong GK, Wang J, Liu B, Yu J, Yang H, Nefedov M, Koriabine M, Dejong PJ, Goodstadt L, Webber C, Dickens NJ, Letunic I, Suyama M, Torrents D, von Mering C, Zdobnov EM, Makova K, Nekrutenko A, Elnitski L, Eswara P, King DC, Yang S, Tyekucheva S, Radakrishnan A, Harris RS, Chiaromonte F, Taylor J, He J, Rijnkels M, Griffiths-Jones S, Ureta-Vidal A, Hoffman MM, Severin J, Searle SM, Law AS, Speed D, Waddington D, Cheng Z, Tuzun E, Eichler E, Bao Z, Flicek P, Shteynberg DD, Brent MR, Bye JM, Huckle EJ, Chatterji S, Dewey C, Pachter L, Kouranov A, Mourelatos Z, Hatzigeorgiou AG, Paterson AH, Ivarie R, Brandstrom M, Axelsson E, Backstrom N, Berlin S, Webster MT, Pourquie O, Reymond A, Ucla C, Antonarakis SE, Long M, Emerson JJ, Betran E, Dupanloup I, Kaessmann H, Hinrichs AS, Bejerano G, Furey TS, Harte RA, Raney B, Siepel A, Kent WJ, Haussler D, Eyras E, Castelo R, Abril JF, Castellano S, Camara F, Parra G, Guigo R, Bourque G, Tesler G, Pevzner PA, Smit A, Fulton LA, Mardis ER, Wilson RK: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.
Bruggeman V, Van As P, Decuypere E: Developmental endocrinology of the reproductive axis in the chicken embryo. Comp Biochem Physiol A Mol Integr Physiol. 2002, 131: 839-846. 10.1016/S1095-6433(02)00022-3.
Smith CA, Sinclair AH: Sex determination: insights from the chicken. Bioessays. 2004, 26: 120-132. 10.1002/bies.10400.
Andrews JE, Smith CA, Sinclair AH: Sites of Estrogen Receptor and Aromatase Expression in the Chicken Embryo. General and Comparative Endocrinology. 1997, 108: 182-190. 10.1006/gcen.1997.6978.
Nishikimi H, Kansaku N, Saito N, Usami M, Ohno Y, Shimada K: Sex differentiation and mRNA expression of P450c17, P450arom and AMH in gonads of the chicken. Mol Reprod Dev. 2000, 55: 20-30. 10.1002/(SICI)1098-2795(200001)55:1<20::AID-MRD4>3.0.CO;2-E.
Nomura O, Nakabayashi O, Nishimori K, Yasue H, Mizuno S: Expression of five steroidogenic genes including aromatase gene at early developmental stages of chicken male and female embryos. J Steroid Biochem Mol Biol. 1999, 71: 103-109. 10.1016/S0960-0760(99)00127-2.
Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology. 2004, 3: Article 3.-10.2202/1544-6115.1027.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B. 1995, 57: 289-300.
UCSC (University of California Santa Cruz) genome browser. http://genome.ucsc.edu/cgi-bin/hgGateway/.
Teranishi M, Shimada Y, Hori T, Nakabayashi O, Kikuchi T, Macleod T, Pym R, Sheldon B, Solovei I, Macgregor H, Mizuno S: Transcripts of the MHM region on the chicken Z chromosome accumulate as non-coding RNA in the nucleus of female cells adjacent to the DMRT1 locus. Chromosome Res. 2001, 9: 147-165. 10.1023/A:1009235120741.
Kamata R, Takahashi S, Morita M: Gene expression of sex-determining factors and steroidogenic enzymes in the chicken embryo: influence of xenoestrogens. Gen Comp Endocrinol. 2004, 138: 148-156. 10.1016/j.ygcen.2004.05.011.
Arnold AP: The gender of the voice within: the neural origin of sex differences in the brain. Current Opinion in Neurobiology. 2003, 13: 759-764. 10.1016/j.conb.2003.10.005.
Dewing P, Shi T, Horvath S, Vilain E: Sexually dimorphic gene expression in mouse brain precedes gonadal differentiation. Brain Res Mol Brain Res. 2003, 118: 82-90. 10.1016/S0169-328X(03)00339-5.
Fridolfsson AK, Cheng H, Copeland NG, Jenkins NA, Liu HC, Raudsepp T, Woodage T, Chowdhary B, Halverson J, Ellegren H: Evolution of the avian sex chromosomes from an ancestral pair of autosomes. PNAS. 1998, 95: 8147-8152. 10.1073/pnas.95.14.8147.
Ellegren H: Dosage compensation: do birds do it as well?. Trends in Genetics. 2002, 18: 25-28. 10.1016/S0168-9525(01)02553-7.
Kuroda Y, Arai N, Arita M, Teranishi M, Hori T, Harata M, Mizuno S: Absence of Z-chromosome inactivation for five genes in male chickens. Chromosome Res. 2001, 9: 457-468. 10.1023/A:1011672227256.
Kuroiwa A, Yokomine T, Sasaki H, Tsudzuki M, Tanaka K, Namikawa T, Matsuda Y: Biallelic expression of Z-linked genes in male chickens. Cytogenet Genome Res. 2002, 99: 310-314. 10.1159/000071609.
McQueen HA, McBride D, Miele G, Bird AP, Clinton M: Dosage compensation in birds. Curr Biol. 2001, 11: 253-257. 10.1016/S0960-9822(01)00070-7.
Bisoni L, Batlle-Morera L, Bird AP, Suzuki M, McQueen HA: Female-specific hyperacetylation of histone H4 in the chicken Z chromosome. Chromosome Res. 2005, 13: 205-214. 10.1007/s10577-005-1505-4.
Raymond CS, Shamu CE, Shen MM, Seifert KJ, Hirsch B, Hodgkin J, Zarkower D: Evidence for evolutionary conservation of sex-determining genes. Nature. 1998, 391: 691-695. 10.1038/35618.
Raymond CS, Parker ED, Kettlewell JR, Brown LG, Page DC, Kusz K, Jaruzelska J, Reinberg Y, Flejter WL, Bardwell VJ, Hirsch B, Zarkower D: A region of human chromosome 9p required for testis development contains two genes related to known sexual regulators. Hum Mol Genet. 1999, 8: 989-996. 10.1093/hmg/8.6.989.
Smith CA, McClive PJ, Western PS, Reed KJ, Sinclair AH: Conservation of a sex-determining gene. Nature. 1999, 402: 601-602. 10.1038/45130.
Kettlewell JR, Raymond CS, Zarkower D: Temperature-dependent expression of turtle Dmrt1 prior to sexual differentiation. Genesis. 2000, 26: 174-178. 10.1002/(SICI)1526-968X(200003)26:3<174::AID-GENE2>3.0.CO;2-J.
Guan G, Kobayashi T, Nagahama Y: Sexually dimorphic expression of two types of DM (Doublesex/Mab-3)-domain genes in a teleost fish, the Tilapia (Oreochromis niloticus). Biochem Biophys Res Commun. 2000, 272: 662-666. 10.1006/bbrc.2000.2840.
Shibata K, Takase M, Nakamura M: The Dmrt1 expression in sex-reversed gonads of amphibians. Gen Comp Endocrinol. 2002, 127: 232-241. 10.1016/S0016-6480(02)00039-4.
Raymond CS, Kettlewell JR, Hirsch B, Bardwell VJ, Zarkower D: Expression of Dmrt1 in the genital ridge of mouse and chicken embryos suggests a role in vertebrate sexual development. Dev Biol. 1999, 215: 208-220. 10.1006/dbio.1999.9461.
Shan Z, Nanda I, Wang Y, Schmid M, Vortkamp A, Haaf T: Sex-specific expression of an evolutionarily conserved male regulatory gene, DMRT1, in birds. Cytogenet Cell Genet. 2000, 89: 252-257. 10.1159/000015626.
Hori T, Asakawa S, Itoh Y, Shimizu N, Mizuno S: Wpkci, encoding an altered form of PKCI, is conserved widely on the avian W chromosome and expressed in early female embryos: implication of its role in female sex determination. Mol Biol Cell. 2000, 11: 3645-3660.
O'Neill M, Binder M, Smith C, Andrews J, Reed K, Smith M, Millar C, Lambert D, Sinclair A: ASW: a gene with conserved avian W-linkage and female specific expression in chick embryonic gonad. Dev Genes Evol. 2000, 210: 243-249. 10.1007/s004270050310.
Pace HC, Brenner C: Feminizing chicks: a model for avian sex determination based on titration of Hint enzyme activity and the predicted structure of an Asw-Hint heterodimer. Genome Biol. 2003, 4: R18-10.1186/gb-2003-4-3-r18.
Wade J, Peabody C, Coussens P, Tempelman RJ, Clayton DF, Liu L, Arnold AP, Agate R: A cDNA microarray from the telencephalon of juvenile male and female zebra finches. J Neurosci Methods. 2004, 138: 199-206. 10.1016/j.jneumeth.2004.04.007.
Chen D, Riedl T, Washbrook E, Pace PE, Coombes RC, Egly JM, Ali S: Activation of estrogen receptor alpha by S118 phosphorylation involves a ligand-dependent interaction with TFIIH and participation of CDK7. Mol Cell. 2000, 6: 127-137. 10.1016/S1097-2765(00)00014-9.
Talukder AH, Mishra SK, Mandal M, Balasenthil S, Mehta S, Sahin AA, Barnes CJ, Kumar R: MTA1 Interacts with MAT1, a Cyclin-dependent Kinase-activating Kinase Complex Ring Finger Factor, and Regulates Estrogen Receptor Transactivation Functions. J Biol Chem. 2003, 278: 11676-11685. 10.1074/jbc.M209570200.
Lee DK, Duan HO, Chang C: From androgen receptor to the general transcription factor TFIIH. Identification of cdk activating kinase (CAK) as an androgen receptor NH(2)-terminal associated coactivator. J Biol Chem. 2000, 275: 9308-9313. 10.1074/jbc.275.13.9308.
Ito S, Takeyama K, Yamamoto A, Sawatsubashi S, Shirode Y, Kouzmenko A, Tabata T, Kato S: In vivo potentiation of human oestrogen receptor alpha by Cdk7-mediated phosphorylation. Genes Cells. 2004, 9: 983-992. 10.1111/j.1365-2443.2004.00777.x.
Larochelle S, Pandur J, Fisher RP, Salz HK, Suter B: Cdk7 is essential for mitosis and for in vivo Cdk-activating kinase activity. Genes Dev. 1998, 12: 370-381.
Matt Kim J, McGaughy JT, Kent Bogle R, Ravnik SE: Meiotic Expression of the Cyclin H/Cdk7 Complex in Male Germ Cells of the Mouse. Biol Reprod. 2001, 64: 1400-1408. 10.1095/biolreprod64.5.1400.
Wallenfang MR, Seydoux G: cdk-7 Is required for mRNA transcription and cell cycle progression in Caenorhabditis elegans embryos. Proc Natl Acad Sci U S A. 2002, 99: 5527-5532. 10.1073/pnas.082618399.
Brenner C: Hint, Fhit, and GalT: function, structure, evolution, and mechanism of three branches of the histidine triad superfamily of nucleotide hydrolases and transferases. Biochemistry. 2002, 41: 9003-9014. 10.1021/bi025942q.
Mosesso P, Piane M, Palitti F, Pepe G, Penna S, Chessa L: The novel human gene aprataxin is directly involved in DNA single-strand-break repair. Cell Mol Life Sci. 2005, 62: 485-491. 10.1007/s00018-004-4441-0.
Gueven N, Becherel OJ, Kijas AW, Chen P, Howe O, Rudolph JH, Gatti R, Date H, Onodera O, Taucher-Scholz G, Lavin MF: Aprataxin, a novel protein that protects against genotoxic stress. Hum Mol Genet. 2004, 13: 1081-1093. 10.1093/hmg/ddh122.
Bieganowski P, Garrison PN, Hodawadekar SC, Faye G, Barnes LD, Brenner C: Adenosine monophosphoramidase activity of Hint and Hnt1 supports function of Kin28, Ccl1, and Tfb3. J Biol Chem. 2002, 277: 10852-10860. 10.1074/jbc.M111480200.
Melo J, Toczyski D: A unified view of the DNA-damage checkpoint. Current Opinion in Cell Biology. 2002, 14: 237-245. 10.1016/S0955-0674(02)00312-5.
Peltomaki P: DNA mismatch repair and cancer. Mutation Research/Reviews in Mutation Research. 2001, 488: 77-85. 10.1016/S1383-5742(00)00058-2.
Trimble JJ, Murthy SC, Bakker A, Grassmann R, Desrosiers RC: A gene for dihydrofolate reductase in a herpesvirus. Science. 1988, 239: 1145-1147.
Itoh Y, Hori T, Saitoh H, Mizuno S: Chicken spindlin genes on W and Z chromosomes: transcriptional expression of both genes and dynamic behavior of spindlin in interphase and mitotic cells. Chromosome Res. 2001, 9: 283-299. 10.1023/A:1016694513051.
Semplici F, Meggio F, Pinna LA, Oliviero S: CK2-dependent phosphorylation of the E2 ubiquitin conjugating enzyme UBC3B induces its interaction with beta-TrCP and enhances beta-catenin degradation. Oncogene. 2002, 21: 3978-3987. 10.1038/sj.onc.1205574.
Balthazart J, Tlemcani O, Ball GF: Do sex differences in the brain explain sex differences in the hormonal induction of reproductive behavior? What 25 years of research on the Japanese quail tells us. Horm Behav. 1996, 30: 627-661. 10.1006/hbeh.1996.0066.
Panzica GC, Viglietti-Panzica C, Balthazart J: The sexually dimorphic medial preoptic nucleus of quail: a key brain area mediating steroid action on male sexual behavior. Front Neuroendocrinol. 1996, 17: 51-125. 10.1006/frne.1996.0002.
Fridolfsson AK, Ellegren H: A simple and universal method for molecular sexing of non-ratite birds. Journal of Avian Biology. 1999, 30: 116-121.
Royal Institute of Technology (KTH) Microarray Center. http://www.biotech.kth.se/molbio/microarray/index.html/.
Stockholm Bioinformatics Center (SBC) Chicken EST identities. http://www.sbc.su.se/~arve/chicken/.
Churchill GA: Fundamentals of experimental design for cDNA microarrays. Nat Genet. 2002, 32 Suppl: 490-495. 10.1038/ng1031.
Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.
Ihaka RGR: R: A Language for Data Analysis and Graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299--314.
ArrayExpress database. http://www.ebi.ac.uk/arrayexpress/.
Quackenbush J, Cho J, Lee D, Liang F, Holt I, Karamycheva S, Parvizi B, Pertea G, Sultana R, White J: The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species. Nucl Acids Res. 2001, 29: 159-164. 10.1093/nar/29.1.159.
NCBI (National Center for Biotechnology Information) GenBank. http://www.ncbi.nlm.nih.gov/Genbank/index.html/.
A.F.A. Smit & P. Green: Institute for Systems Biology. http://www.repeatmasker.org/cgi-bin/WEBRepeatMasker.
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2005, 33 Database Issue: D34-8.
Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LSL: UniProt: the Universal Protein knowledgebase. Nucl Acids Res. 2004, 32: D115-119. 10.1093/nar/gkh131.
Kent WJ: BLAT---The BLAST-Like Alignment Tool. Genome Res. 2002, 12: 656-664. 10.1101/gr.229202. Article published online before March 2002.
Karolchik D, Baertsch R, Diekhans M, Furey TS, Hinrichs A, Lu YT, Roskin KM, Schwartz M, Sugnet CW, Thomas DJ, Weber RJ, Haussler D, Kent WJ: The UCSC Genome Browser Database. Nucl Acids Res. 2003, 31: 51-54. 10.1093/nar/gkg129.
Hubbard T, Andrews D, Caccamo M, Cameron G, Chen Y, Clamp M, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T, Down T, Durbin R, Fernandez-Suarez XM, Gilbert J, Hammond M, Herrero J, Hotz H, Howe K, Iyer V, Jekosch K, Kahari A, Kasprzyk A, Keefe D, Keenan S, Kokocinsci F, London D, Longden I, McVicker G, Melsopp C, Meidl P, Potter S, Proctor G, Rae M, Rios D, Schuster M, Searle S, Severin J, Slater G, Smedley D, Smith J, Spooner W, Stabenau A, Stalker J, Storey R, Trevanion S, Ureta-Vidal A, Vogel J, White S, Woodwark C, Birney E: Ensembl 2005. Nucleic Acids Res. 2005, 33 Database Issue: D447-53.
Lee Y, Tsai J, Sunkara S, Karamycheva S, Pertea G, Sultana R, Antonescu V, Chan A, Cheung F, Quackenbush J: The TIGR Gene Indices: clustering and assembling EST and known genes and integration with eukaryotic genomes. Nucleic Acids Res. 2005, 33 Database Issue: D71-4.
Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR: The Pfam protein families database. Nucleic Acids Res. 2004, 32: D138-41. 10.1093/nar/gkh121.
Tatusova TA, Madden TL: BLAST 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol Lett. 1999, 174: 247-250. 10.1111/j.1574-6968.1999.tb13575.x.
Muller YP, Janovjak H, Miserez AR, Dobbie S: Processing of Gene Expression Data generated by Quantitative Real-Time RT-PCR. Biotechniques. 2002, 32: 1372-4, 1376, 1378-9..
Wong ML, Medrano JF: Real-time PCR for mRNA quantitation. Biotechniques. 2005, 39: 75-85.
We want to thank Professor Leif Andersson and the KTH Microarray Resource Center for providing the chicken cDNA arrays. Dr Peter Savolainen at the KTH Microarray Resource Center is thanked for valuable discussions concerning the quality of the spotted microarrays and for sequencing cDNA clones. Raili Engdahl and Lena Norgren are acknowledged for excellent technical assistance. This work was supported by the Swedish EPA (the ReproSafe program) and the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning (Formas).
All authors conceived and designed the experiments, but KK designed the microarray study. AM performed the experiment related to the animal material. KK, BS and non-authors (see Acknowledgments) carried out the microarray and PCR analysis. BS and KK analyzed the data. BB and LD contributed reagents, materials, and analysis tools. BS, KK and MS drafted the manuscript with significant contributions from all authors. All authors read and approved the final manuscript.
Birger Scholz, Kim Kultima contributed equally to this work.
Electronic supplementary material
Additional File 3: Graphs of PCR results. The average mean normalized expression (MNE) for females (red) and males (blue) of selected gene (listed alphabetically) on ed4 and/or ed8, analyzed by real-time PCR (n = 5). Error bars indicate standard error of the mean. For calculations, see Methods. (DOC 5 MB)
Additional File 4: Alignment of ABWT. A. UCSC genome browser window showing the alignment of ABWT and three additional annotation tracks. 'Non-Chicken RefSeq Genes' shows the alignment of non-chicken RefSeq sequences for Ube2r2 from mouse (Mus musculus) and zebrafish (Danio rerio). 'Chicken mRNAs from GenBank' describes the alignment of additional chicken mRNAs in this region. 'Repeating Elements by RepeatMasker ' shows the parts of the genome region masked for repetitive sequences. This UCSC Genome Browser window can be accessed by searching with 'chrW_random:175,450-180,581' on the UCSC chicken genome browser site . B. Alignment data for ABWT against the chicken genome. (JPEG 63 KB)
Additional File 5: Species comparisons of the UBE2R2 protein. Species comparison between chicken (Gg), human (Hs), mouse (Mm), Xenopus leavis/tropicalis (Xl/Xt) and Saccharomyces cerevisiae (Sc) of the UBE2R2 protein using the protein sequences listed in additional file 7. A comparison between Pfam domains (representing protein domains or conserved protein regions) and the number of positive amino acid matches against the chicken sequence is shown (see Methods). (DOC 40 KB)
Additional File 6: Species comparisons of proteins associated with TFIIH complex. Species comparison between chicken (Gg), human (Hs), mouse (Mm), Xenopus leavis/tropicalis (Xl/Xt) and Saccharomyces cerevisiae (Sc) of TFIIH proteins using the protein sequences listed in additional file 7. A comparison between Pfam domains (representing protein domains or conserved protein regions) and the number of positive amino acid matches against the chicken sequence is shown (see Methods). (DOC 368 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.