Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Sterility and Gene Expression in Hybrid Males of Xenopus laevis and X. muelleri

  • John H. Malone,

    Affiliation Department of Biology, The University of Texas at Arlington, Arlington, Texas, United States of America

  • Thomas H. Chrzanowski,

    Affiliation Department of Biology, The University of Texas at Arlington, Arlington, Texas, United States of America

  • Pawel Michalak

    To whom correspondence should be addressed. E-mail: michalak@uta.edu

    Affiliation Department of Biology, The University of Texas at Arlington, Arlington, Texas, United States of America

Abstract

Background

Reproductive isolation is a defining characteristic of populations that represent unique biological species, yet we know very little about the gene expression basis for reproductive isolation. The advent of powerful molecular biology tools provides the ability to identify genes involved in reproductive isolation and focuses attention on the molecular mechanisms that separate biological species. Herein we quantify the sterility pattern of hybrid males in African Clawed Frogs (Xenopus) and apply microarray analysis of the expression pattern found in testes to identify genes that are misexpressed in hybrid males relative to their two parental species (Xenopus laevis and X. muelleri).

Methodology/Principal Findings

Phenotypic characteristics of spermatogenesis in sterile male hybrids (X. laevis x X. muelleri) were examined using a novel sperm assay that allowed quantification of live, dead, and undifferentiated sperm cells, the number of motile vs. immotile sperm, and sperm morphology. Hybrids exhibited a dramatically lower abundance of mature sperm relative to the parental species. Hybrid spermatozoa were larger in size and accompanied by numerous undifferentiated sperm cells. Microarray analysis of gene expression in testes was combined with a correction for sequence divergence derived from genomic hybridizations to identify candidate genes involved in the sterility phenotype. Analysis of the transcriptome revealed a striking asymmetric pattern of misexpression. There were only about 140 genes misexpressed in hybrids compared to X. laevis but nearly 4,000 genes misexpressed in hybrids compared to X. muelleri.

Conclusions/Significance

Our results provide an important correlation between phenotypic characteristics of sperm and gene expression in sterile hybrid males. The broad pattern of gene misexpression suggests intriguing mechanisms creating the dominance pattern of the X. laevis genome in hybrids. These findings significantly contribute to growing evidence for allelic dominance in hybrids and have implications for the mechanism of species differentiation at the transcriptome level.

Introduction

Biological species remain cohesive by a lack of gene flow between interspecific populations and the mechanisms that maintain this pattern are manifest by various forms of reproductive isolation. Postzygotic reproductive isolation is characterized by dysfunctional phenotypes observed in F1 interspecific hybrids including but not limited to inviability and/or sterility resulting in decreased gene flow between species. Haldane (1922) [1] observed that the heterogametic sex (XY or ZW) typically suffers the most dysfunctional effects of hybridization. Haldane's rule has been shown to be broadly applicable across diverse groups of animals suggesting that understanding the basis for Haldane's rule provides a key towards understanding the mechanisms for how species become and/or remain reproductively isolated.

Much work, using mainly forward genetic approaches, has focused on the genetic basis of postzygotic reproductive isolation and Haldane's rule. These studies have provided support for both dominance effects and faster male evolution as the main mechanisms generating Haldane's rule; however only recently have studies focused on gene expression as related to sterile hybrids and this opens a new avenue towards understanding the proximate causes of reproductive isolation [2][12]. By starting at the level of phenotype, analyzing the transcriptome found in both species and compared to the dysfunctional hybrid, this reverse genetics approach provides the candidate loci that contribute to reproductive isolation, identifies targets that have evolved disparately between the two species, and allows a test of what evolutionary forces create gene misexpression in dysfunctional hybrids. Although correlational in nature, this approach advances our understanding of reproductive isolation by generating testable hypotheses for which future functional experiments can be designed.

Frogs of the genus Xenopus offer an exciting new system to explore the expression basis of reproductive isolation. Xenopus are characterized genomically by allopolyploidization and range from diploid (n = 20) to dodecaploid (n = 108) numbers of chromosomes [13]. Most species can be crossed to produce viable progeny and males from lab produced and wild caught hybrids are sterile whereas hybrid females are fertile [14][21]. Sex-reversal experiments and sex ratios from backcross progeny have established that a dominant allele in females of X. laevis, X. muelleri, and other species determines sex and therefore Xenopus have ZW sex determination [22][27]. Additionally, no morphologically distinct Z or W sex chromosomes have been identified in Xenopus and therefore Xenopus have homomorphic sex chromosomes [28], [29]. As Xenopus females are heterogametic and males are homogametic, yet in interspecies crosses the males are consistently sterile and females are fully or partially fertile, Xenopus provide an exception to Haldane's rule. Given the genetic peculiarities as well as the variety of genetic and genomic tools available for Xenopus as a model system, Xenopus offer an excellent system to address questions related to gene expression, reproductive isolation, and speciation.

We explored gene expression and reproductive isolation in hybrids of Xenopus laevis and X. muelleri, both tetraploid species (2n = 36) which hybridize in the wild [21]. To do so, we focused on the phenotypic characteristics of spermatogenesis in sterile male hybrids (X. laevis x X. muelleri). Affymetrix Xenopus laevis Genome Arrays were used to assay the transcriptome in testes combined with a correction for sequence divergence from genomic hybridizations that allowed discovery of the broad pattern of misexpression as well as the identification of candidate genes involved in the sterility phenotype. This approach importantly allows a correlation of the broad gene expression pattern to the phenotypic characteristics observed for sperm in sterile hybrids and identifies the loci misexpressed in hybrids relative to the two parental species.

Results

Sperm Abundance

One explanation for why male Xenopus laevis x X. muelleri are sterile could be due to phenotypic defects associated with the process of spermatogenesis. The production of amphibian sperm is a complex physiological process involving six key stages of differentiation that include primary spermatogonia, secondary spermatogonia, primary spermatocytes, secondary spermatocytes, and spermatids [30]. Additionally, it is widely accepted that the process of spermatogenesis is under hormonal control and that the injection of gonadotropins stimulates sperm production [31], [32]. If hybrid male Xenopus are sterile due to phenotypic defects associated with spermatogenesis, we predict that characteristics of sperm quality (i.e., abundance, motility, and morphology) should be different in hybrids compared to the parental species. We tested this hypothesis by injecting males with human chorionic gonadotropin hormone (hCG) to assay the effect of hormone induced stimulation on spermatogenesis in hybrids and the two parental species and then compared sperm characteristics of these injected males with uninjected sexually mature males.

Sperm abundance was quantified using a novel sperm assay that allowed detection of live, dead, and undifferentiated sperm cells. Testes were homogenized and then incubated with fluorescent dyes that intercalate with DNA of sperm cells based on whether the cell is live or dead. Counts of live, dead, and undifferentiated sperm cells were made using epifluorescence microscopy and we tested the null hypothesis that sperm abundance was the same in hybrids compared to the two parental species.

There was a dramatic difference in the abundance of sperm in Xenopus laevis compared to hybrids (Fig. 1a–c). About 40 times more sperm cells were found in X. laevis compared to hybrids (F1,9 = 135.4; P = 0.000) and there was no effect of the hCG treatment on the total number of sperm cells (F1,9 = 0.001; P = 0.978; Fig. 2). Only one X. muelleri was available for analysis and the abundance of sperm for this injected X. muelleri male (22,900 sperm/microliter) was comparable to the abundance of sperm for injected X. laevis (Mean sperm = 26,833 sperm/microliter) but not comparable to injected hybrids (Mean = 738 sperm/microliter).

thumbnail
Figure 1. Visualization of live (green) and dead (red) sperm in sperm density assay using dual emission filter for SYBR14 and propidium iodide.

Xenopus laevis (A); hybrid (B); X. muelleri (C); and brightfield combined with fluorescent image of undifferentiated cells of hybrids (D).

https://doi.org/10.1371/journal.pone.0000781.g001

thumbnail
Figure 2. The number of total sperm cells in Xenopus laevis (black circles) and hybrids (striped squares) in uninjected sexually mature males and hCG injected sexually mature males.

Error bars represent±1 standard error.

https://doi.org/10.1371/journal.pone.0000781.g002

Proportions of live and dead sperm cells revealed an effect of hCG treatment on Xenopus laevis only. The proportion of live sperm cells was about 40% greater in normal X. laevis compared to hybrids (F1,8 = 7.13; P = 0.028) and there was an affect of hCG treatment on the proportion of live sperm cells (F1,8 = 15.75; P = 0.004). This hCG treatment effect can be explained mainly by the X. laevis specific response to hCG treatment (Taxa x Treatment interaction F1,8 = 6.31; P = 0.036) as the proportion of live cells decreased by about 50% (because the proportion of dead cells increased from 30% to 50% from the hCG treatment) in X. laevis but did not change in hybrids (Fig. 3). The proportion of live and dead cells for X. muelleri (Live = 0.52; Dead = 0.48) was higher compared to injected X. laevis (Mean proportion Live X. laevis = 0.22; Mean proportion Dead X. laevis = 0.78) but was closer to the mean proportion of live and dead cells for normal X. laevis (Mean proportion Live = 0.69; Mean proportion Dead = 0.30)

thumbnail
Figure 3. Proportion of live (green circles), dead (red squares), and undifferentiated (black triangles) sperm cells in Xenopus laevis (A) and hybrids (B) in normal/uninjected compared to hCG injected males.

Error bars represent±1 standard error.

https://doi.org/10.1371/journal.pone.0000781.g003

There was a striking difference in the proportion of undifferentiated sperm cells in Xenopus laevis compared to hybrids (Fig. 1b, d). In X. laevis there were practically no undifferentiated sperm cells but in hybrids the majority (about 50%) of the sperm suspension was comprised of undifferentiated sperm (F1,8 = 46.073; P = 0.000). The hCG treatment had no effect on the number of undifferentiated sperm cells (F1,8 = 0.013; P = 0.911; Fig. 3). There was no significant number of undifferentiated sperm cells for the sample of X. muelleri.

Sperm Motility

Motility of sperm is important to successful fertilization in amphibian species. Xenopus, like many amphibian species, have external fertilization and sperm become motile and swim to fertilize eggs only after entering a lower osmolality environment [33]. Hybrids produce a lower abundance of mature sperm and while sterility may be explained in part by this lower abundance, sterility could be complete if sperm in hybrids failed to become motile. We tested the null hypothesis that sperm motility was the same in hybrids compared to the two parental species by counting the number of motile and immotile sperm following activation.

There were about 70 times more motile sperm in X. laevis compared to hybrids (F1,5 = 188.3; P = 0.000) and there was no effect of hCG treatment on the number of motile sperm (F1,5 = 0.258; P = 0.633). Despite the large difference in number of motile sperm cells, proportions of motile to immotile sperm were the same between X. laevis and hybrids (F1,5 = 0.570; P = 0.484; Fig. 4). There was 12% motile and 88% immotile sperm out of a total of 238 sperm cells observed in X. muelleri. This was close to the mean percentage of sperm for X. laevis (Mean = 18% motile and Mean = 82% immotile).

thumbnail
Figure 4. Numbers of motile (green circles) and immotile (red x's) sperm in Xenopus laevis and hybrids in normal/uninjected compared to hCG injected males.

Error bars represent±1 standard error.

https://doi.org/10.1371/journal.pone.0000781.g004

Sperm Morphology

Abnormal sperm morphology is frequently associated with infertility [34] and misshapen or abnormal sperm morphology in hybrid males would suggest an explanation for sterility in hybrids. We measured the area of sperm cells as a proxy for size and tested the null hypothesis that the size of mature sperm was the same in hybrids compared to the parental species.

Hybrid sperm were dramatically larger compared to the sperm of Xenopus laevis and X. muelleri (F1,51 = 74.0; P = 0.000). Hybrid sperm were about 28 µm2 larger than Xenopus laevis (Bonferroni corrected P = 0.000) and 25 µm2 larger than X. muelleri sperm (Bonferroni corrected P = 0.000). The size of sperm did not differ between X. laevis and X. muelleri (Bonferroni corrected P = 0.513; Fig. 5).

thumbnail
Figure 5. Comparison of sperm area (µm2) in Xenopus laevis, X. muelleri, and hybrids and representative brightfield images of sperm.

Error bars represent±1 standard error and * denotes significance from a Bonferroni multiple comparison test. Hybrids have larger sperm compared to X. laevis and X. muelleri but sperm area does not differ between X. laevis and X. muelleri.

https://doi.org/10.1371/journal.pone.0000781.g005

Gene Expression

Our analyses of sperm characteristics in hybrids compared to X. laevis and X. muelleri revealed several phenotypic differences that may contribute to the sterility pattern. We next analyzed the transcriptional pattern in hybrids and the two parental species to explore the proximate mechanism of gene expression and its contribution to the phenotype of hybrids. This approach allowed us to identify loci that may be involved in generating the sterility phenotype and simultaneously examine the broad pattern of gene expression in hybrids compared to the parental species.

We used Affymetrix Xenopus laevis Genome Arrays to generate a transcriptional profile of gene expression in testes of hybrids, X. laevis, and X. muelleri. The Affymetrix microarray is designed for X. laevis and it is widely accepted that hybridizing RNA from a heterospecific species to a microarray designed for a related species can have a dramatic impact on the signal recovered from microarrays [35][37]. To control for this effect we directly assayed sequence divergence by hybridizing genomic DNA from Xenopus laevis and X. muelleri each separately to the Affymetrix Xenopus laevis GeneChip® Genome Array. By taking the ratio of hybridization intensity of X. muelleri/X. laevis for each probe on the array we then screened out probes that did not hybridize properly due to sequence divergence in X. muelleri. To be conservative we set the lower ratio of hybridization intensity to 0.99 and explored variation in the number of probes eliminated at a variety of upper thresholds with the idea that higher intensity for X. muelleri is not as damaging compared to probes that have a hybridization signal lower for X. muelleri when compared to X. laevis (Fig. 6). In examining the threshold variation on the number of probesets remaining, we chose to analyze the datasets generated at the 1.01 and 1.10 ratio level. Both datasets provided a similar general pattern but the effect was more prominent for the less stringent threshold. We report here the results generated using an upper threshold of 1.10 in subsequent analyses of gene expression which resulted in the removal of 226,841 individual probes and provided 11,485 probesets for further analysis. We tested the null hypothesis that the expression level for each of the 11,485 probesets was the same in three separate contrasts (hybrids vs. X. laevis; hybrids vs. X. muelleri; and X. laevis vs. X. muelleri).

thumbnail
Figure 6. Number of probesets remaining at various thresholds from comparing the hybridization intensity of Xenopus muelleri vs. X. laevis.

https://doi.org/10.1371/journal.pone.0000781.g006

Microarray analysis revealed an asymmetrical pattern of gene misexpression between hybrids and the two parental species. Only 1.2% of genes (142/11,485) in our sample of the transcriptome were misexpressed between Xenopus laevis and hybrids whereas about 35% of genes (3,995/11,485) were misexpressed between X. muelleri and hybrids (Fig. 7). There were more genes upregulated in hybrids relative to X. laevis (92 vs. 50; G = 12.61; df = 1; P<0.001) but there were more genes upregulated in X. muelleri compared to hybrids (2,236 vs. 1759; G = 57.1; df = 1; P<0.001). Complete results for each of the three main contrasts can be found in Tables S1, S2 and S3. The top thirty most misexpressed transcripts in the two contrast tests were mainly dominated by EST sequences with little functional information available but our analysis of gene expression in the testis would suggest that these EST sequences have testis related function. Among the annotated transcripts, many are known to have functions in the process of spermatogenesis, spermiogenesis, or testis related functions or are involved in regulating polymerase II transcription (Table S4, S5, S6 and S7). There were 56 transcripts recovered as differentially expressed in both contrasts (Table 1) and these transcripts may have a more crucial role in the sterility pattern in hybrid males of Xenopus.

thumbnail
Figure 7. Volcano plots from FDR corrected t-tests of statistical significance (vertical axis) against magnitude of expression change (horizontal axis), where each point corresponds to a gene/transcript.

Expression change (fold-change) is defined as a log2-transformed ratio of mean nonhybrid to mean hybrid expression level. (A) Xenopus laevis vs. Hybrids; (B) Xenpous muelleri vs. Hybrids. The red line denotes FDR adjusted alpha 0.05. The horizontal deviation from 0 towards the right or left reflects hybrid underexpression or overexpression, respectively.

https://doi.org/10.1371/journal.pone.0000781.g007

thumbnail
Table 1. Candidate loci (56 total) common to both contrasts (X. laevis vs hybrids; X. muelleri vs. hybrids) divided according to corresponding expression behavior (hybrid>(laevis = muelleri); hybrid<(laevis = muelleri); laevis>hybrid>muelleri; laevis<hybrid<muelleri).

https://doi.org/10.1371/journal.pone.0000781.t001

There was also a dramatic difference in expression between the two parental species. About 60% of genes were differentially expressed between X. laevis and X. muelleri (6956/11485). Of these differentially expressed transcripts, about 40% (2824/6956) were upregulated in X. laevis relative to X. muelleri and 60% (4132/6956) were upregulated in X. muelleri relative to X. laevis.

Comparing the overlap in genes misexpressed in individual hybrid contrasts to the three classes of expression behavior between X. laevis and X. muelleri (X. laevis>X. muelleri; X. laevis<X. muelleri; X. laevis = X. muelleri) revealed that in general hybrids have an intermediate level of expression compared to the two parental species (Table 2). For example, of the 92 genes upregulated in hybrids relative to X. laevis; 92% were also upregulated in X. muelleri relative to X. laevis. Similarly, of the 1,759 genes that were upregulated in hybrids relative to X. muelleri, 92% were upregulated in X. laevis compared to X. muelleri. These results suggest a general pattern of intermediate expression in hybrids and are consistent with a semi-dominant model of expression difference even despite the strong asymmetrical pattern of misexpression in hybrids compared to the two parental species.

thumbnail
Table 2. Comparison in the overlap of transcripts recovered as differentially expressed from the two contrasts with hybrids (Xenopus laevis vs. hybrids and X. muelleri vs. hybrids) and the interspecies contrast (Xenopus laevis vs. X. muelleri).

https://doi.org/10.1371/journal.pone.0000781.t002

We further asked whether the observed expression pattern is unbiased by sequence divergence still unaccounted for by the masking procedure or distorted by the procedure itself. To address this question, we used genomic sequence information available for Xenopus tropicalis (JGI, v3.0) as an outgroup to help identify highly conserved sequences among distantly related species. X. tropicalis belongs to a clade (Silurana) that have diverged from X. laevis and X. muelleri for more than 70 million years, compared to a >20-million year divergence between X. laevis and X. muelleri [38], [39]. We blasted (BLASTN 2.2.12 [40]) all 495,232 X. laevis probe sequences provided by Affymetrix against the X. tropicalis genome and selected 20 probesets with the lowest E-values (10−6 to 0.04 averaged across the 32 probes per probeset) corresponding to the most conserved sequences (100-97.2% mean identity) between X. laevis and X. tropicalis (Table S8). As these two species are more distantly related than X. laevis and X. muelleri are to each other, it is reasonable to assume that these sequences will tend to be conserved in X. muelleri as well (no X. muelleri genome information is available), and thus their expression pattern will not be confounded by sequence divergence. We also reasoned that if the masking procedure was valid and unbiased, the expression patterns from the mask and the most conserved probesets would be similar. Indeed, the mean differences between normalized expression levels from these two methods were almost identical: 0.503 and 0.427 for the X. laevis-hybrids contrast and 1.040 and 1.091 for the X. muelleri-hybrids contrast (Kolmogorov-Smirnov tests, P>0.05). This approximately twofold expression difference between the X. laevis-hybrid and the X. muelleri-hybrid comparisons from the conserved probesets was significant (t-test, P = 0.029), confirming the asymmetric pattern of expression differences. This asymmetry persisted across all masking criteria, in the absence of the mask (X. laevis vs. hybrids: 180/15491 = 1.2%, X. muelleri vs. hybrids: 9345/15491 = 60.3%), and was robust to changes in the method of normalization (RMA vs. MAS 5.0 scaling – results not shown), providing additional support for the expression pattern.

Discussion

Our analysis of sperm abundance, motility, and morphology in hybrids suggests that the process and control of spermatogenesis is severely perturbed in hybrids. Hybrids have about 40 times less total sperm compared to parental species. Mature sperm in hybrids are capable of motility and the highly reduced abundance of sperm suggests that these numbers may be insufficient to fertilize eggs in vivo. Additionally, differential live/dead staining of sperm cells in hybrids showed that the majority of cells consisted of undifferentiated sperm suggesting an overall pattern of arrested development of spermatogenesis and/or spermiogenesis. Cells that have the characteristic shape of sperm in hybrids are larger compared to X. laevis and X. muelleri raising the possibility that these sperm cells have abnormal numbers of chromosomes. Kobel (1996) [18] suggested that sperm cells of hybrids were aneuploid and our results may be consistent with these observations. Future work focusing on chromosome content of hybrid sperm would aid in determining the end result of meiosis in hybrids.

Sperm production in hybrids did not respond to hCG hormone treatment whereas the proportion of dead cells increased dramatically in X. laevis. The dosage of hCG used in the experiment is typical for stimulating reproduction in Xenopus but our results suggest the possibility that the dose may be too high. hCG triggers the production of LH and FSH which are involved in the differentiation of sperm [31], [32]. Overstimulation by excessive hCG could result in sperm death and may explain the increased numbers of dead sperm in X. laevis but why hybrids fail to respond to hCG treatment remains an intriguing question. If receptors for hCG are more sensitive in hybrids, an abnormally high dose of hCG may cause these receptors to become inactivated leading to a lack of response by the process of spermatogenesis. Alternatively, hormone receptors in hybrids may fail to function properly due to the hybrid genetic background which could lead to a lack of response to hCG treatment. This later hypothesis would implicate a possible mechanism for the overall different pattern of spermatogenesis observed in hybrids. Spermatogenesis, a process tightly controlled by hormonal interactions with the hypothalamic-pituitary-gonadal axis [30], may be misregulated due to improper interactions with hormones and receptors. This hypothesis has received no investigation in studies of reproductive isolation despite considerable attention directed to the process of spermatogenesis as an explanation for Haldane's rule [41], [42]. The lack of response to hCG in hybrids suggests that hormonal regulation of spermatogenesis may contribute to reproductive isolation in hybrids and support for this hypothesis could be generated by assaying receptor function in hybrids compared to parental species.

The genes misexpressed in hybrids may offer clues to the loci involved in the sterility phenotype characteristic of hybrid males and therefore potential loci that contribute to reproductive isolation between X. laevis and X. muelleri. Many of the most misexpressed loci consist of EST targets and therefore currently there is little functional information available. One EST of interest though has been identified to be the transcription factor TFIIE complex which is a transcription factor important to the function of RNA polymerase II transcription. This EST was downregulated 4 times in hybrids relative to X. laevis suggesting that Pol II transcription may be impacted by the lack of TFIIE in hybrids. Interestingly, this gene was not found to be misexpressed compared to X. muelleri.

Many interesting genes related to spermatogenesis in other organisms appear in the candidate gene lists. For example, type 2 retinaldehyde dehydrogenase (RALDH2) was found to be 3.2 times lower in hybrids compared to X. laevis. This gene catalyzes the important developmental modulator retinoic acid and is exclusively expressed in mouse testis [43]. Caesin kinase I (CKIe), a gene involved in protein amino acid phosphoryalation through the utilization of ATP but not GTP was downregulated 4.3 times in hybrids relative to X. laevis and a unique form of CKle is expressed in rat testis [44].

The most misexpressed gene in hybrids relative to X. laevis was neuropeptide Y (NPY). NPY is one of the most abundant and widespread nueropeptides in mammals and has been suggested to play a role in controlling reproductive function. In particular, NPY potentiates the release of LH and FSH in response to GnRH [45], [46] and is modulated by testosterone. NPY is predominantly synthesized in the hypothalamic arcuate nucleus but also expressed in the testes (mostly in Leydig cells) of mouse and rat [47], [48]. Importantly, there is a strong association between testicular development and the expression of NPY [48]. NPY was upregulated 16 times in hybrids relative to X. laevis and also was overexpressed in X. muelleri relative to X. laevis. Other upregulated genes in hybrids relative to X. laevis and that have documented expression patterns in the testes include PAK5 (p21 activated kinase 2 plays a role in Sertoli-germ cell anchoring dynamics of rats; [49] and dynll2 (dynein light chain-1 was highly expressed in mouse and rat testis; [50]).

Examining genes misexpressed in hybrids compared to X. muelleri reveals that there are also many interesting genes related to the function of the testis in other organisms. Occludin which was upregulated about 9 times more in X. muelleri compared to hybrids has an important role in the formation of tight junctions surrounding Sertoli cells and seminiferous tubules in mammals [51] and salamanders [52]. Goosecoid is an important transcription factor during early development of Xenopus and was found to be about 8 times overexpressed in X. muelleri relative to hybrids. A related gene (goosecoid-like) is expressed primarily in the brain and primoridial germ cells in mouse [53]. B1, a X. borealis specific gene included on the X. laevis microarray, was upregulated about 7 times in X. muelleri compared to hybrids and plays an important role in the assembly of TFIIIA and expression of 5S RNA by RNA polymerase III [54]. Two tubulin related genes (similar to alpha-tubulin and beta-tubulin) were overexpressed 140 and 85 times respectively in hybrids relative to X. muelleri. Beta tubulin 2 has a testis restricted expression profile in Aedes aegypti and Drosophila [55], [56]. Another gene of interest is cyclin B2, upregulated 70 times in hybrids compared to X. muelleri. Cyclins have a demonstrated role in spermatogenesis in the eel Anguilla japonica [57] and mouse [58] and play a role in cell cycle regulation particularly in response to cancerous aberrations in germ cells [59]. Ferritin was found to be upregulated 68 times in hybrids compared to X. muelleri. Recently mitochondrial ferritin has been discovered to have a testis specific expression profile in mammals and flies [60] leading to the possibility that we recovered a signal of this ferritin product in Xenopus. Finally, HMG box protein, a transcriptional repressor that was upregulated about 62 times in hybrids relative to X. muelleri, has been found to be expressed strongly in mouse testis [61].

The gene expression data suggest two major patterns for expression behavior in hybrids relative to the two parental species. First, for misexpressed genes in hybrids the overwhelming preponderance of misexpressed genes follow a semi-dominant model of expression behavior because hybrid expression is intermediate and/or additive compared to expression differences between the two species (Table 2). Hybrids of these two species of Xenopus then appear to follow a very different overall pattern of misexpression compared to hybrids in other organisms. For example, hybrids of Drosophila have distinctly nonadditive expression behavior because hybrids have more misexpression compared to the genes misexpressed between the two parental species [3][4]. The same can be said for hybrids of maize [62], Arabidopsis [8], [63][65], and Senecio [7]. In the case of Xenopus, more genes are misexpressed between species than for hybrids compared to each parental species.

The second major pattern observed from the microarray results is that there was a difference in the expression pattern of hybrids compared to the two parental species. Surprisingly, a relatively few number of genes were misexpressed between X. laevis and hybrids (142/11,485) while a substantial number of genes were misexpressed between X. muelleri and hybrids (3995/11485). In effect, adult testes of hybrids have an expression pattern very similar to X. laevis but substantially different compared to X. muelleri. These results suggest that most genes in the hybrid genetic background follow a dominant pattern of expression because the overall level of expression in hybrids is equal to that of X. laevis only. This pattern of expression in hybrid testes differs from analyses of proteins in hybrids of X. laevis and X. borealis (a related species to X. muelleri). Both species specific proteins of X. laevis and X. borealis are expressed in hybrid ovary [66] and species-specific copies of several allozyme loci appear to both be expressed in various developmental stages of tadpoles, adult heart, and adult liver [67]. However, LDH isozyme in early embryos resembled that of the maternal species (X. laevis victorianus) and this persisted through metamorphosis [68]. These studies suggest that the allelic dominance observed in the expression profile of hybrid testes is not distributed in other tissue types but rather may be confined to the testes.

One interpretation of these results is that the collective molecular processes that occur in adult frog testes (i.e. spermio- and spermatogenesis) are under the influence of massive genomic imprinting and specifically, the maternal but not the paternal copies of genes are expressed. This would provide an explanation for the similarity of hybrid expression profiles to X. laevis and the dissimilarity compared to X. muelleri because the mother of these hybrids is X. laevis and not X. muelleri. No evidence for imprinting has been found for Xenopus but this hypothesis would require that imprinting may occur in Xenopus and that maternal imprinting is the controller of the expression pattern in whole frog testes [69]. One fruitful avenue to explore this hypothesis would be to analyze the expression profiles from testes of reciprocal cross hybrids (X. muelleri x X. laevis). If the expression profile becomes more similar to X. muelleri then this would be evidence of strong maternal effects, possible related to genomic imprinting, controlling the expression pattern in frog testes.

An alternative mechanism that may explain the dissimilar expression profiles is that in hybrids the paternal copy of chromosomes is eliminated or quiescent at the expression level. Hybrids then would express only the maternally contributed X. laevis portion of the genome but this would require that the maternally contributed portion be upregulated two fold to match the expression profile of X. laevis. The mechanisms involved here are difficult to explain but one possibility may be that paralogous gene copies, normally repressed, become expressed thereby providing expressed products to match that of X. laevis. Genomic elimination of paternal copies does occur in insects and frogs [70][74]. For example in the treefrog, gene silencing and monoallelic dominance of one genome occurs in crosses of tetraploid Hyla versicolor and diploid H. arborea [72]. In F1 triploid hybrids, several allozyme loci had unexpected frequencies suggesting that in some cases the paternal gene was silenced (i.e. Mpi-2) or that the maternal allele was silenced and the single paternal allele had double expression (i.e. Got-1). The implication here is that the foreign copy can be recognized and eliminated, and given the polyploid nature of Xenopus a similar mechanism may be operating. Indeed, amphibians have remarkably examples of genomic exclusion. One well documented example consists of the hybridogenetic system of the Rana esculenta complex in Europe where the entire paternal genome is eliminated during oogenesis [75].

Bringing together two divergent genomes into a common genetic background can induce a variety of genomic changes including genome instability, changes in chromatin, and transcriptome shock [7], [76][79]. In particular, gene silencing appears to be one major result of bringing together two divergent genomes into a common genomic environment via hybridization. For example, interesting patterns regarding rRNA gene expression have been established in interspecific hybrids of Xenopus. Both Xenopus laevis and X. borealis (a species closely related to X. muelleri) display two nucleoli per cell during embryogenesis. Hybrids between these two species have one nucleolus in the majority of cells suggesting that the nucleolar organizer is inactivated due to genomic incompatibility between the two species [80]. More specifically, in hybrids between Xenopus laevis and X. borealis (in the cited papers mistakenly referred to as X. muelleri; [81], the rRNA of X. laevis is preferentially transcribed and X. borealis is repressed [82]. Repression of X. borealis rRNA is basically complete until the swimming tadpole stage during which a low level of X. borealis rRNA is detectable. In adults, the level of expression of the X. borealis copy ranges from zero to significant amounts indicating the presence of much variation between individual hybrids. This phenomenon, first called “differential amphiplasty” and currently termed “nucleolar dominance”, is an epigenetic phenomena that is controlled at the level of transcription and occurs in a wide variety of interspecific hybrids including the plant taxa Crepis [83], [84], Salix [85], Ribes [86], [87], Solanum [88], Hordeum [89][92], Triticum [93][96], Agropyron [97], Brassica [98], [99], and Arabidopsis [100] as well as hybrids of Drosophila [101], [102], and mammal somatic cell hybrids of mouse and human [103][109]. Nucleolar dominance operates on tens of millions of base pairs of chromosomal DNA and this large scale inactivation mechanism is only rivaled by X chromosome inactivation in mammalian cells [110].

Can the mechanisms of nucleolar dominance, a pattern specifically related to the RNA polymerase I transcription machinery of rRNA and not with RNA polymerase II transcription of protein coding genes provide explanations for the asymmetrical transcriptomic expression pattern observed in hybrids of Xenopus? The mechanisms that discriminate maternal or paternal rRNA genes remain unclear but three hypotheses have been proposed to account for the differential silencing of species specific alleles: the species-specific transcription factor hypothesis, the enhancer-imbalance hypothesis, and the chromatin imprinting hypothesis [111][113].

Species-Specific Transcription Factor Hypothesis

The DNA sequences coding for ribosomal proteins (18S, 5.8S, and 25/26S) are highly conserved but intergenic rDNA sequences are highly divergent amongst eukaryotes [112][115]. According to the species-specific transcription factor hypothesis, because of this divergent evolution rRNA gene transcription can occur only if co-evolved species-specific transcription factors bind to corresponding regulatory sequences. When genomes from two different species are brought together, there is a failure of rDNA transcription because of a lack of appropriate specific transcription factors. Cell-free transcription systems of mouse-human hybrid cell lines revealed that the promoter for mouse or the promoter for human rRNA would not function in a cell-free extract made from another species [116][119]. In plants, the tomato rRNA gene promoter does not work when transfected into Arabidopsis protoplasts and the tobacco rRNA gene promoter does not work when placed into a bean cell-free extract [120], [121]. These findings suggest that transcription factors had co-evolved with their corresponding regulatory DNA sequences because transcription function was erased when forced into a heterospecific environment.

One prediction of this co-evolution model is that the amount of sequence divergence would be inversely related to how effective the mechanisms of the species-specific transcription hypothesis operate and this appears to be the case. When more closely related species are brought together into a common genomic environment, rRNA gene promoters are fully functional. For example, in Brassica napus, a hybrid derived from B. rapa and B. oleracea, only B. rapa rRNA is transcribed reflecting a pattern of nucleolar dominance but both rRNA genes were transcribed when transfected into protoplasts of the other species [122]. The same pattern occurs in Arabidopsis suecica, a hybrid of A. thaliana and A. arenosa; A. thaliana and A. arenosa rRNA genes are expressed in both species and transfected A. thaliana rRNA is active when introduced into the hybrid [100]. These congeneric interspecific hybrids are more closely related compared to the previous examples and here species specific transcription factors work fine in a heterospecific background. In the case of Xenopus, it is possible at least for some loci that transcription factors may fail to bind to their corresponding regulatory targets causing a silencing of underdominant genes because the two species examined in this study are substantially divergent. The estimated divergence time for the last common ancestor of X. laevis and X. muelleri is 21–35 mya [38], [39], [123]. However, this mechanism could operate on a genome wide scale in adult testes only if a key transcription factor, or factors, derived from maternally contributed X. laevis genes and involved in the regulation of X. muelleri chromatin was sufficiently divergent to misregulate the proper expression of X. muelleri alleles. Allele-specific imprinting at the level of chromatin has been recently reported in Arabidopsis [124] and in the case of Xenopus evidence for this hypothesis would be obtained by identifying chromatin remolding factors that operate from the maternal genome to regulate paternal copy chromatin and/or paternal copy transcription.

Enhancer Imbalance Hypothesis

The enhancer imbalance hypothesis proposes that differences in the amount of intergenic sequences provide a competitive environment for the transcription of rRNA genes. If particular intergenic sequences are more effective at controlling transcription then these sequences will outcompete sequences contributed by the other species creating an imbalance in the way that enhancers operate to sequester transcription factors. In Xenopus, the repetitive DNA elements located upstream of the rRNA promoter act as orientation and position independent enhancers of transcription [125][128]. In X. laevis, each intergenic spacer contains a 60/81 bp repeat with a 42 bp core element that shares 90% sequence similarity with the gene promoter. In X. borealis, complete 60/81 bp repeats are absent but the 42 bp core element is present suggesting that the 42 bp core element controls the enhancer effect [111], [129][131]. In hybrids, only the rRNA genes of X. laevis are transcribed [80], [81], [132] and this appears to be due to the more numerous enhancer elements located in X. laevis sequence compared to X. borealis which sequester a critical transcription factor necessary for transcription [126], [133], [134]. Co-injection of rRNA minigenes with complete spacers into oocytes of X. laevis and X. borealis showed that X. laevis minigenes were preferentially transcribed and this occurred even when injecting recombinant constructs that had spacer and promoter regions swapped [134]. These experiments support the notion that species specific differences in intergenic regions, and not gene promoters, are responsible for the failure of the X. borealis allele to be transcribed in hybrids.

The enhancer imbalance hypothesis stems directly from biochemical interactions of transcription factors and the enhancer elements that regulate gene expression and therefore follows models of cis-trans evolution that have recently been proposed [135], [136]. Differences in genomic content between species then may be consistent with this model because differing amounts of intergenic DNA and/or enhancers may cause problems with transcribing one allele from one of the two species. Examining published C-values for Xenopus (available on the Animal Genome Database-http://www.genomesize.com/), the genomic content of Xenopus laevis and X. muelleri differs by about 430 Mb (Mean X. laevis 3231 Mb; Mean X. muelleri 3660 Mb; t = −2.76; df = 1; P = 0.02). By combining these two genomes into the same genomic environment, different enhancer content resulting from differences in genome size, may contribute to the patterns observed in gene expression of X. laevis, X. muelleri, and their F1 hybrid.

Chromatin Imprinting Hypothesis

DNA methylation and covalent modifications of histone tails within nucleosomes (acetylation, phosphorylation, sumoylation, ubiquination, and methylation) have an impact on chromatin structure, gene transcription, and epigenetic information [137]. Indeed modifications to chromatin are applied and removed in a highly specific fashion leading to the proposition of a histone code that is read by chromatin-associated factors [138]. The chromatin imprinting hypothesis proposes that changes in the regulation of chromatin state explain the dominance seen in the preferential transcription of rRNA genes [111]. DNA hypomethylation and histone hyperacetylation correlate with transcriptional activity and DNA hypermethylation and histone hypoacetylation correlate with transcriptional silencing and therefore allele specific regulation of chromatin state may render some gene copies to be transcribed while rendering others to be silent [139].

There is a striking correlation between genome size and the direction of silencing in interspecific hybrids. On average, the parental species with the smaller genome is always dominant whereas under-dominant loci came from the smaller genome of the two parental species [111]. The argument here is that larger genomes have more heterochromatin and when combined into the same genomic environment via hybridization this increased heterochromatin content changes the spatial structure of the genome within the nucleus leading to preferential silencing of one parental copy. In effect, the positional configuration of chromatin in the hybrid nucleus is imbalanced and the genome with the larger content will be subjected to heterochromatinization via chromatin remolding mechanisms [140]. As stated above, the genomic content differs between X. laevis and X. muelleri and the smaller of the two genomes (X. laevis) is the species for which hybrids are most similar at the expression level. Additionally, several genes involved in chromatin maintenance are misexpressed in hybrids compared to X. muelleri but not in hybrids compared to X. laevis (i.e. Ube2e2 that has ubiquitin-protein ligase activity and SET which is involved in chromatin remolding).

Our analysis of the phenotype of hybrids revealed a dramatically lower abundance of sperm in hybrids, increased numbers of undifferentiated sperm cells, a lack of response to hCG treatment, larger sperm compared to both parental species and even despite these phenotypic abnormalities, the presence of a few motile sperm. Coupling these phenotypic data with microarray analyses provides loci that may lead to the sterile phenotypic condition of hybrids. These data are important for identifying genes germane to postzygotic reproductive isolation between species and we present the first lists of these loci for Xenopus. At the broad gene expression level, more genes are misexpressed between species compared to hybrids. This is different compared to gene expression analyses of hybrids in other organisms where in general more misexpression occurs in hybrids. Genes that are misexpressed in hybrids compared to each species follow a semi-dominant model of expression behavior because most genes have an intermediate level of expression in hybrids. Despite the semi-dominance of misexpressed genes, there was strong pattern of dominance for the X. laevis genome because hybrids had fewer genes misexpressed compared to X. laevis than compared to X. muelleri. This pattern implies a silencing mechanism for the paternally inherited X. muelleri alleles and has implications for the evolution of species differentiation at the expression level [141]. This silencing mechanism may be a natural phenomenon related to amphibian testis function or alternatively could be a genuine response to the hybrid genetic background. Most interesting is that this pattern of silencing, reminiscent of the widespread phenomenon of nucleloar dominance, has also been documented in Arabidopsis allotetraploid hybrids [65]. We propose three hypotheses derived from large scale silencing of Xenopus rRNA genes which could explain the species-level dominance pattern. By doing so we focus much needed attention on the molecular mechanisms that promote gene transcription and in particular how these mechanisms relate to the expression basis of reproductive isolation in Xenopus.

Future research will be required to address the intriguing question of why Xenopus do not conform to Haldane's rule, this longstanding generalization of evolutionary biology. There is a common misconception that Haldane's rule does not apply to taxa that lack heteromorphic sex chromosomes, which is an example of confusing the phenomenon (heterozygous sex is afflicted) with its cause (dominance patterns due to heteromorphic sex chromosomes). In fact, sterility patterns (but not inviability patterns) in various taxa lacking heteromorphic sex chromosomes conform well to Haldane's rule [142][144], and suggest a role of other evolutionary mechanisms, such as faster male evolution [42].

Materials and Methods

Experimental Crosses

Four biological replicate hybrid clutches were produced at The University of Texas at Arlington. The female parent for each of these crosses was Xenopus l. laevis and these originate from the Cape region of South Africa or from the laboratory of Jacques Robert at The University of Rochester. The male parents were X. muelleri and originate from Nkambeni area of Swaziland and were kindly provided by R. C. Tinsley. For the expression analysis, we deliberately used X. laevis from genetically diverse sources (Xenopus Express, J. Robert's Lab, and R. C. Tinsley's Lab) to increase intragroup heterogeneity resulting in the lower rate of the null hypothesis rejections and more conservative statistical inference.

Crosses were conducted in the following manner: each female was injected with 500 IU human chorionic gonadotropin (hCG; Choluron or Carolina Biological Supply) and placed into a container with 10 L of dechlorinated water. Males were injected with 200 IU hCG and then paired with individual females. Each hybrid clutch was raised to sexual maturity and clutches were maintained in 12L:12D photoperiod at 25°C air temperature. Animals were fed Xenopus pellets from Xenopus Express Inc. every two days.

Sperm Assays

We quantified characteristics of the sperm in hybrids and the two parental species. From each of three biological replicate clutches we selected two sexually mature hybrids. Sexual maturity was determined by the presence of dark colored nuptial pads indicative of circulating androgens and characteristic of sexual maturity [145]. One hybrid was injected with 150 IU hCG and the other was not injected with hCG in order to test the effect of hCG treatment on sperm abundance. For X. laevis, we used F1 males from Xenopus Express and whose parents originate from the Cape of South Africa and/or the lab of J. Roberts at the University of Rochester and selected animals based on sexual maturity. Again, one of the two animals was injected with 150 IU hCG and then left overnight in 10 L of dechlorinated water. For X. muelleri we selected one F1 male X. muelleri that was produced at the University of Texas at Arlington and whose parents originate from the Nkambeni area of Swaziland. This male was injected with 150 IU hCG.

Males were killed by immersion in MS-222 and testes were harvested and placed in 500 µl of DeBoer's Solution (110 mM NaCl; 1.3 mM CaCl2; 1.3 mM KCl). Testes were homogenized using a handheld pestle to form a stock sperm solution. Ten µl of the stock sperm solution was aliquoted and further diluted 3∶1 with distilled water. The 3∶1 dilution ratio was found to provide the maximum number of motile sperm in X. laevis [33] and was used to assay for sperm motility. Motile sperm and nonmotile sperm were enumerated during a 2 minute interval. This procedure was repeated four subsequent times/individual to provide a metric of sperm motility. Two individuals per treatment (injected and uninjected) and per taxa (X. laevis and hybrids) were used for sperm motility assays.

We assayed the abundance of live and dead sperm cells using the LIVE/DEAD Sperm Viability Kit from Molecular Probes (Invitrogen, Inc.). The LIVE/DEAD Sperm Kit is designed specifically for assaying live and dead sperm cells. Five µl of a 50-fold dilution of the stock SYBR 14 dye and 5 µl of 2.4 mM propidium iodide were added to the stock sperm suspension followed by a 10 minute incubation at 36°C for each dye. Stained cells were collected on black polycarbonate filters (1.0 µm pore-size, GE Osmonics) and enumerated at a magnification of 165× using epifluorescence microscopy (Olympus BH-2). Briefly, 10 µl of stained sperm were added to 2 ml of particle-free distilled-water (0.2 µm pore-size filtered) contained in the filter tower of a 10- place filtration unit (Hoeffer). Following filtration, filters were placed on slides, a drop of immersion oil added to the surface, and a cover slip placed on top. Live and dead sperm cells were enumerated from fifteen separate, randomly-selected visual fields. Images were captured using a digital camera (Olympus DP70).

We measured the area of sperm acrosomes in square micrometers as a proxy for the size for spermatozoa using images captured from placing 10 ul of sperm stock solution onto a microscope slide. Thirteen sperm cells from hybrids, 21 sperm cells from X. laevis, and 20 sperm cells from X. muelleri were selected at random from the images and the area in square micrometers of each sperm cell was measured using ImageJ software (NIH).

We constructed general linear models in SYSTAT v. 8.0 (SAS, 1998) to test the null hypothesis that there was no difference in the total number of sperm cells, the proportion of live, dead, and undifferentiated sperm cells, and the proportion of motile and nonmotile sperm between taxa and between the hCG treatment. An ANOVA was used to test the hypothesis that there was no difference in the size of spermatozoa between treatments.

Microarray Experiments

RNA was extracted from adult testis in Xenopus laevis (n = 4), hybrids of X. laevis x. X. muelleri (n = 4) and X. muelleri (n = 3). Adults were euthanized with MS-222 and testes were immediately removed, placed in RNA extraction solution, and homogenized. RNA was recovered using GeneHunter and Ambion RiboPure total RNA kits. Samples of RNA were checked for purity by examination of the 28S and 18S ribosomal RNA bands from denaturing gel electrophoresis, by 260/280 ratios from scans with a Nanodrop ND 1000 spectrophotometer, and by readouts of the Agilent Bioanalyzer. Total RNA samples were prepared and hybridized to Affymetrix Xenopus laevis GeneChip Genome Arrays at the University of Texas Southwestern Medical Center Microarray Array Core Facility following standard Affymetrix protocols. Affymetrix Microarray Analysis Suite (MAS) v.5.0 was used to scan and process each microarray chip. The signals of quality control and poly(A) transcripts revealed that hybridizations were of high quality in all chips. Quality control probesets (i.e., spike in and housekeeping genes) were removed in subsequent statistical analyses. Non-unique probesets (i.e. interrogating different transcript variants from same genes) that represent <4% of the entire array were not masked out.

Hybridizing RNA from a heterospecific species to a microarray designed for a related species can have a dramatic impact on the signal recovered from microarrays [35][37]. Consequently, we directly assayed sequence divergence by hybridizing genomic DNA from Xenopus laevis and X. muelleri each separately to the Affymetrix Xenopus laevis GeneChip Genome Array. Genomic DNA was extracted from liver from eight individuals/species (four male and four female) using a QIAGEN DNeasy kit. Genomic DNA was quantified using a Nanodrop ND 1000 spectrophotometer and 2.5 µg/µl from each of the eight individuals was pooled to produce 20 µg genomic DNA for subsequent fragmentation. The 20 µg pooled genomic DNA was fragmented by a DNase I digestion using 2.5 µl 10× One-Phor-All Buffer (Amersham Biosciences), 0.015 U/µg DNase I (Amersham Biosciences), and 14.125 µl H2O to form a 25 µl reaction. Reactions were incubated at 37°C for 10 minutes and the DNase I was inactivated by incubating the reaction at 100°C for 10 minutes. Reactions were visualized on 2.5% agarose gels stained with SYBR Green to confirm that fragments were in the 50–200 bp size range necessary for hybridization to Affymetrix microarrays. Fragments were labeled on the 3′ termini and hybridized to the X. laevis Affymetrix Genome Array following standard protocols recommended by Affymetrix.

We used an electronic mask to eliminate probes that behaved poorly due to sequence divergence. To do so, we modified the Xspecies perl script (Version 1.1) in [146] to incorporate information gained by comparing the ratio of genomic hybridization intensity of X. muelleri to X. laevis for each probe pair. The Xspecies perl script selects a probe-set when one or more PM probe-pair (s) meets user-specified criteria of hybridization intensity and creates a new probe-definition file (.cdf) that contains only those probe-pairs that meet the user-specified criteria. This, in theory, should create a hybridization signal that interrogates gene expression in a less biased manner compared to not eliminating these probes. To conservatively guard against sequence divergence we set the lower ratio of hybridization intensity to 0.99 and explored variation in the number of probes eliminated at a variety of upper thresholds with the idea that higher intensity for X. muelleri is not as damaging compared to probes that have a hybridization signal lower for X. muelleri when compared to X. laevis (Fig. 6). In examining the threshold variation on the number of probesets remaining, we chose to analyze the .cdf generated at the 1.01 and 1.10 level. Both sets provided a similar general pattern but the effect was more prominent for the less stringent threshold. We report here the results generated using an upper threshold of 1.10 in subsequent analyses of gene expression which resulted in the removal of 226,841 individual probes and provided 11,485 probesets for further analysis.

Data Analysis

Using the chip definition files created from the Xspecies perl script analyses, we conducted two separate comparisons to uncover patterns of differential expression between Xenopus laevis and hybrids, X. muelleri compared to hybrids, and X. laevis compared to X. muelleri. First, the Xenopus laevis and hybrid data (filtered with the Xspecies mask) were normalized using RMAexpress [147] with default parameters for background correction and quantile normalization. These RMA normalized data were then imported into the R statistical environment and tested for differences in expression between X. laevis and hybrids for each of the 11,485 genes using a moderated t-statistic based on an empirical Bayes method in the limma package found in Bioconductor [148]. The TopTable function was then used to output the FDR-adjusted P-values for differential expression. We normalized X. muelleri and hybrid chips together using RMA and repeated the analyses to uncover differential expression between X. muelleri and hybrids. Finally, we normalized X. laevis and X. muelleri chips together using RMA and repeated the analysis to uncover genes misexpressed between the two species. This later test was used to discover patterns of expression behavior in hybrids compared to interspecies expression behavior.

Supporting Information

Table S1.

Genes differentially expressed between Xenopus laevis and hybrids. Probe ID is the Affymetrix reference number for a particular probeset. Columns B–E are expression values after RMA normalization for Xenopus laevis and columns F–I are hybrids. MeanLaevis is the mean expression value of X. laevis and MeanHybrid is the mean expression value for hybrids. L-H is the fold change in expression (MeanLaevis - MeanHybrid). t is the moderated t-statistic. P.Value is the unadjusted P value obtained from the empirical Bayes function. adj.P.Val is the FDR corrected P value according to Benjamini and Hochberg (1995). B is the B-statistic which is the log-odds that a gene is differentially expressed. Description is the annotation information for a probeset as given by Affymetrix.

https://doi.org/10.1371/journal.pone.0000781.s001

(0.09 MB XLS)

Table S2.

Genes differentially expressed between Xenopus muelleri and hybrids. Probe ID is the Affymetrix reference number for a particular probeset. Columns B–D are expression values after RMA normalization for Xenopus laevis and columns E–H are hybrids. MeanMuelleri is the mean expression value of X. muelleri and MeanHybrid is the mean expression value for hybrids. M-H is the fold change in expression (MeanMuelleri - MeanHybrid). t is the moderated t-statistic. P.Value is the unadjusted P value obtained from the empirical Bayes function. adj.P.Val is the FDR corrected P value according to Benjamini and Hochberg (1995). B is the B-statistic which is the log-odds that a gene is differentially expressed. Description is the annotation information for a probeset as given by Affymetrix.

https://doi.org/10.1371/journal.pone.0000781.s002

(2.19 MB XLS)

Table S3.

Genes differentially expressed between Xenopus muelleri and X. laevis. Probe ID is the Affymetrix reference number for a particular probeset. Columns B–D are expression values after RMA normalization for Xenopus muelleri and columns E–H are X. laevis. MeanMuell is the mean expression value of X. muelleri and MeanLaev is the mean expression value for hybrids. M-L is the fold change in expression (MeanMuell - MeanLaev). t is the moderated t-statistic. P.Value is the unadjusted P value obtained from the empirical Bayes function. adj.P.Val is the FDR corrected P value according to Benjamini and Hochberg (1995). B is the B-statistic which is the log-odds that a gene is differentially expressed. Description is the annotation information for a probeset as given by Affymetrix.

https://doi.org/10.1371/journal.pone.0000781.s003

(3.78 MB XLS)

Table S4.

Top 30 candidate transcripts upregulated in X. laevis and differentially expressed between X. laevis and hybrid. Expression values are in log2 scale; SD = standard deviation of expression values. P values are adjusted according to FDR moderated t-tests.

https://doi.org/10.1371/journal.pone.0000781.s004

(0.09 MB DOC)

Table S5.

Top 30 candidate transcripts upregulated in hybrids and differentially expressed between X. laevis and hybrid. Expression values are in log2 scale; SD = standard deviation of expression values. P values are adjusted according to FDR moderated t-tests.

https://doi.org/10.1371/journal.pone.0000781.s005

(0.09 MB DOC)

Table S6.

Top 30 candidate transcripts upregulated in X. muelleri and differentially expressed between X. muelleri and hybrid. Expression values are in log2 scale; SD = standard deviation of expression values. P values are adjusted according to FDR moderated t-tests.

https://doi.org/10.1371/journal.pone.0000781.s006

(0.09 MB DOC)

Table S7.

Top 30 candidate transcripts upregulated in hybrids and differentially expressed between X. muelleri and hybrid. Expression values are in log2 scale; SD = standard deviation of expression values. P values are adjusted according to FDR moderated t-tests.

https://doi.org/10.1371/journal.pone.0000781.s007

(0.09 MB DOC)

Table S8.

Top 20 transcripts with the highest sequence conservation between Affymetrix probesets and Xenopus tropicalis genome. Probe ID is the Affymetrix reference number for a particular probeset. Columns B–C are mean differences between RMA-normalized expression values of Xenopus laevis and hybrids (B) and Xenopus muelleri and hybrids (C). Columns D and E are absolute values from columns B and C, respectively. Column F is mean % sequence identity of sequences between Affymetrix probesets and Xenopus tropicalis genome generated by BLAST. Column G is average alignment length from BLAST. Column H is average mismatch number of the alignments from F and G. Column I is mean E-value generated by BLAST. J column contains Affymetrix descriptions of transcripts.

https://doi.org/10.1371/journal.pone.0000781.s008

(0.02 MB XLS)

Acknowledgments

This work was conducted following the protocols of the University of Texas-Arlington Animal Care Committee (Protocol No. A05.001). We thank Professor R. C. Tinsley for donating specimens of Xenopus muelleri for use in this study, members of the Genome Biology Group at The University of Texas-Arlington for insightful discussions, and editorial comments from Jacob Malcom, all of which improved the quality of the manuscript.

Author Contributions

Conceived and designed the experiments: PM JM. Performed the experiments: JM. Analyzed the data: JM. Contributed reagents/materials/analysis tools: PM TC. Wrote the paper: JM.

References

  1. 1. Haldane JBS (1922) Sex ratio and unisexual sterility in animal hybrids. J Genet 12: 101–109.
  2. 2. Reiland J, Noor MAF (2002) Little qualitative RNA misexpression in sterile male F1 hybrids of Drosophila pseudoobscura and D. persimilis. BMC Evol Biol 2: 16.
  3. 3. Michalak P, Noor MAF (2003) Genome-wide patterns of expression in Drosophila pure species and hybrid males. Mol Biol Evol 20: 1070–1076.
  4. 4. Ranz JM, Namgyal K, Gibson G, Hartl DL (2004) Anomalies in the expression profile of interspecific hybrids of Drosophila melanogaster and Drosophila simulans. Genome Research 14: 373–379.
  5. 5. Barbash DA, Lorigan JG (2006) Lethality in Drosophila melanogaster/Drosophila simulans species hybrids is not associated with substantial transcriptional misregulation. J Exp Zool (Mol Dev Evol) 306B:
  6. 6. Haerty W, Singh RS (2006) Gene regulation divergence is a major contributor to the evolution of Dobzhansky-Muller incompatibilities between species of Drosophila. Mol Biol Evol 23: 1707–1714.
  7. 7. Hegarty MJ, Barker GL, Wilson ID, Abbott RJ, Edwards KJ, et al. (2006) Transcriptome shock after interspecific hybridization in Senecio is ameliorated by genome duplication. Curr Biol 16: 1652–1659.
  8. 8. Josefsson C, Dilkes B, Comai L (2006) Parent-dependent loss of gene silencing during interspecies hybridization. Curr Biol 16: 1322–1328.
  9. 9. Lai Z, Gross BL, Zou Y, Andrews J, Rieseberg LH (2006) Microarray analysis reveals differential gene expression in hybrid sunflower species. Mol Ecol 15: 1213–1227.
  10. 10. Michalak P (2006) RNA world – the dark matter of evolutionary genomics. J Evol Biol 19: 1768–1774.
  11. 11. Moehring AJ, Teeter KC, Noor MAF (2007) Genome-wide patterns of expression in Drosophila pure species and hybrid males. II. Examination of multiple-species hybridizations, platforms, and life cycle stages. Mol Biol Evol 24: 137–145.
  12. 12. Ortíz-Barrientos D, Counterman BA, Noor MAF (2007) Gene expression divergence and the origin of hybrid dysfunctions. Genetica 129: 71–81.
  13. 13. Kobel HR, Du Pasquier L (1986) Genetics of polyploidy Xenopus. Trends Genet 2: 310–315.
  14. 14. Kobel HR, Du Pasquier L (1975) Production of large clones of histocompatible, fully identical Clawed Toads (Xenopus). Immunogenetics 2: 87–91.
  15. 15. Kobel HR, Du Pasquier L, Tinsley RC (1981) Natural hybridization and gene introgression between Xenopus gilli and Xenopus laevis laevis (Anura: Pipidae). J Zool Lond 194: 317–322.
  16. 16. Picker MD (1985) Hybridization and habitat selection in Xenopus gilli and Xenopus laevis in the southwestern Cape Province. Copeia 1985: 574–580.
  17. 17. Picker MD, de Villiers AL (1989) The distribution and conservation status of Xenopus gilli (Anura: Pipidae). Biol Conserv 49: 169–183.
  18. 18. Kobel HR (1996) Allopolyploid speciation. In: Tinsley RC, Kobel HR, editors. The biology of Xenopus. Oxford: Clarendon Press. pp. 391–401.
  19. 19. Kobel HR, Loumont C, Tinsley RC (1996) The extant species. In: Tinsley RC, Kobel HR, editors. The biology of Xenopus. Oxford: Clarendon Press. pp. 9–33.
  20. 20. Picker MD, Harrison JA, Wallace D (1996) Natural hybridization between Xenopus laevis and X. gilli in the south-western Cape Province, South Africa. In: Tinsley RC, Kobel HR, editors. The biology of Xenopus. Oxford: Clarendon Press. pp. 61–71.
  21. 21. Fischer WJ, Koch WA, Elepfandt A (2000) Sympatry and hybridization between the clawed frogs Xenopus laevis laevis and Xenopus muelleri (Pipidae). J Zool Lond 252: 99–107.
  22. 22. Mikamo K, Witschi E (1963) Functional sex-reversal in genetic females of Xenopus laevis, induced by implanted testes. Genetics 48: 1411–1421.
  23. 23. Mikamo K, Witschi E (1964) Masculinization and breeding of the WW Xenopus. Experientia 20: 622–623.
  24. 24. Witschi E (1971) Mechanisms of sexual differentiation. In: Hamburg M, Barrington E, editors. Hormones in development. New York: Appleton Century Crofts. pp. 601–618.
  25. 25. Kobel HR (1985) Sex determination in polyploid Xenopus. S Afr J Sci 81: 205–206.
  26. 26. Hayes TB (1998) Sex determination and primary sex differentiation in amphibians: genetic and developmental mechanisms. J Exp Zool 281: 373–399.
  27. 27. Miyata S, Koike S, Kubo T (1999) Hormonal reversal and the genetic control of sex differentiation in Xenopus. Zool Sci 16: 335–340.
  28. 28. Tymowska J, Kobel HR (1972) Karyotype analysis of Xenopus muelleri (Peters) and Xenopus laevis (Daudin), Pipidae. Cytogenetics 11: 270–278.
  29. 29. Schmid M, Steinlein C (1991) Chromosome banding in Amphibia, XVI. High-resolution replication banding patterns in Xenoups laevis. Chromosoma 101: 123–132.
  30. 30. Scheltinga DM, Jamieson BGM (2003) Spermatogenesis and the mature spermatozoon: form, function and phylogenetic implications. In: Jamieson BGM, editor. Reproductive biology and phylogeny of anura. Plymouth: Science Publishers, Inc.. pp. 119–251.
  31. 31. Pozzi AG, Rosemblit C, Ceballos NR (2006) Effect of human gonadotropins on spermiation and androgen biosynthesis in the testis of the toad Bufo arenarum (Amphibia, Anura). J Exp Zool 305A: 96–102.
  32. 32. Rosemblit C, Pozzi AG, Ceballos NR (2006) Relationship between steroidogenesis and spermiation in Rana catesbeiana and Leptodactylus ocellatus. J Comp Physiol B 176: 559–566.
  33. 33. Christensen JR, Pauli BD, Richardson JS, Bishop CA, Elliott J (2004) Effects of pH and dilution on African clawed frog (Xenopus laevis) sperm motility. Can J Zool 82: 555–563.
  34. 34. Pilder SH, Olds-Clarke P, Orth JM, Jester WF, Dugan L (1997) Hst7: a male sterility mutation perturbing sperm motility, flagellar assembly, and mitochondrial sheath differentiation. J Androl 18: 663–667.
  35. 35. Nuzhdin SV, Wayne ML, Harmon KL, McIntyre LM (2004) Common pattern of evolution of gene expression level and protein sequence in Drosophila. Mol Biol Evol 21: 1308–1317.
  36. 36. Gilad Y, Rifkin SA, Bertone P, Gerstein M, White KP (2005) Multi-species microarrays reveal the effect of sequence divergence on gene expression profiles. Genome Research 15: 674–680.
  37. 37. Gilad Y, Borevitz J (2006) Using DNA microarrays to study natural variation. Curr Opin Genet Dev 16: 553–558.
  38. 38. Evans BJ, Kelley DB, Tinsley RC, Melnick DJ, Cannatella DC (2004) A mitochondrial DNA phylogeny of African clawed frogs: phylogeography and implications for polyploidy evolution. Mol Phylo Evol 33: 197–213.
  39. 39. Evans BJ, Kelley DB, Melnick DJ, Cannatella DC (2005) Evolution of RAG-1 in polyploidy clawed frogs. Mol Biol Evol 22: 1193–1207.
  40. 40. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.
  41. 41. Wu C-I, Davis AW (1993) Evolution of postmating reproductive isolation: the composite nature of Haldane's rule and its genetic basis. Am Nat 142: 187–212.
  42. 42. Wu C-I, Johnson NA, Palopoli MF (1996) Haldane's rule and its legacy: why are there so many sterile males? Trends Ecol Evol 11: 281–284.
  43. 43. Hsu LC, Chang WC, Yoshida A (2000) Mouse type-2 retinaldehyde dehydrogenase (RALDH2): genomic organization, tissue-dependent expression, chromosome assignment and comparison to other types. Biochim Biophys Acta 1492: 289–293.
  44. 44. Graves PR, Haas DW, Hagedorn CH, DePaoli-Roach AA, Roach PJ (1993) Molecular cloning, expression, and characterization of a 49-kilodalton casein kinase I isoform from rat testis. J Biol Chem 268: 6394–6401.
  45. 45. Bauer-Dantoin AC, Mcdonald JK, Levine JE (1992) Neuropeptide Y potentiates luteinizing hormone (LH)-releasing hormone-induced LH secretion only under conditions leading to preovulatory LH surges. Endocrinology 131: 2946–2952.
  46. 46. Knox KL, Bauer-Dantoin AC, Levine JE, Schwartz NB (1995) Unmasking of neuropeptide-Y inhibitory effects on in vitro gonadotropin secretion from pituitaries of metestrous, but not proestrous rats. Endocrinology 136: 187–194.
  47. 47. Kanzaki M, Fujisawa M, Okuda Y, Okada H, Arakawa S, et al. (1996) Expresssion and regulation of neuropeptide Y messenger ribonucleic acid in cultured immature rat Leydig and Sertoli cells. Endocrinology 137: 1249–1257.
  48. 48. Terado M, Nomura M, Mineta K, Fujimoto N, Matsumoto T (2006) Expression of neuropeptide Y gene in mouse testes during testicular development. Asian J Androl 8: 443–449.
  49. 49. Siu MK, Wong CH, Lee WM, Cheng CY (2005) Sertoli-germ cell anchoring junction dynamics in the testis are regulated by an interplay of lipid and protein kinases. J Biol Chem 280: 25029–25047.
  50. 50. Wang RA, Zhao M, Meistrich ML, Kumar R (2005) Stage-specific expression of dynein light chain-1 and its interacting kinase, p21-activated kinase-1, in rodent testes: implications in spermiogenesis. J Histochem Cytochem 53: 1235–1243.
  51. 51. Gye MC (2004) Expression of occludin in canine testis and epididymis. Reprod Dom Anim 39: 43–47.
  52. 52. Jin Y, Uchida I, Eto K, Kitano T, Abe S-I (2007) Size-selective junctional barrier and Ca2+ - independent cell adhesion in the testis of Cynops pyrrhogaster: expression and function of occludin. Mol Reprod Dev DOI 10.1002:
  53. 53. Galili N, Nayak S, Epstein JA, Buck CA (2000) Rnf4, a RING protein expressed in the developing nervous and reproductive systems, interacts with Gscl, a gene within the DiGeorge critical region. Develop Dyn 218: 102–111.
  54. 54. Kaulen H, Pognonec P, Gregor PD, Roeder RG (1991) The Xenopus B1 factor is closely related to the mammalian activator USF and is implicated in the developmental regulation of TFIIIA gene expression. Mol Cell Biol 11: 412–424.
  55. 55. Nielsen MG, Caserta JM, Kidd SJ, Phillips CM (2006) Functional constraint underlies 60 million year stasis of Dipteran testis-specific β-tubulin. Evol Dev 8: 23–29.
  56. 56. Smith RC, Walter MF, Hice RH, O'Brochta DA, Atkinson PW (2007) Testis-specific expression of the β2 tubulin promoter of Aedes aegypti and its application as a genetic sex-separation marker. Insect Mol Biol 16: 61–71.
  57. 57. Kajiura-Kobayashi H, Kobayashi T, Nagahama Y (2004) The cloning of cyclin B3 and its gene expression during hormonally induced spermatogenesis in the teleost, Anguilla japonica. Biochem Biophys Res Comm 323: 288–292.
  58. 58. Brandeis M, Rosewell I, Carrington M, Crompton T, Jacobs MA, et al. (1998) Cyclin B2-null mice develop normally and are fertile whereas cyclin B1-null mice die in utero. Proc Natl Acad Sci USA 95: 4344–4349.
  59. 59. Diederichs S, Bäumer N, Schultz N, Hamra KF, Schrader MG, et al. (2005) Expression patterns of mitotic and meiotic cell cycle regulators in testicular cancer and development. Int J Cancer 116: 207–217.
  60. 60. Missirlis F, Holmberg S, Georgieva T, Dunkov BC, Rouault TA, et al. (2006) Characterization of mitochondrial ferritin in Drosophila. Proc Natl Acad Sci USA 103: 5893–5898.
  61. 61. Smith JM, Bowles J, Wilson M, Koopman P (2004) HMG box transcription factor gene Hbp1 is expressed in germ cells of the developing mouse testis. Develop Dyn 230: 366–370.
  62. 62. Auger DL, Gray AD, Ream TS, Kato A, Coe EH Jr. , et al. (2005) Nonadditive gene expression in diploid and triploid hybrids of maize. Genetics 169: 389–397.
  63. 63. Comai L, Madlung A, Josefsson C, Tyagi A (2003) Do different parental ‘heteromes’ cause genomic shock in newly formed allopolyploids? Philos Trans R Soc Lond B Biol Sci 358: 1149–1155.
  64. 64. Wang J, Tian L, Madlung A, Lee HS, Chen M, et al. (2004) Stochastic and epigenetic changes of gene expression in Arabidopsis polyploids. Genetics 167: 1961–1973.
  65. 65. Wang J, Tian L, Lee H-S, Wei NE, Jiang H, et al. (2006) Genomewide nonadditive gene regulation in Arabidopsis allotetraploids. Genetics 172: 507–517.
  66. 66. De Robertis EM, Black P (1979) Hybrids of Xenopus laevis and Xenopus borealis express proteins from both parents. Develop Biol 68: 334–339.
  67. 67. Wall DA, Blackler AW (1974) Enzyme patterns in two species of Xenopus and their hybrids. Develop Biol 36: 379–390.
  68. 68. Vonwyl E (1983) Expression of the lactate dehydrogenase genes of Xenopus species and interspecies hybrids during early development. Comp Biochem Physiol B 76: 17–21.
  69. 69. de la Casa-Esperón E, Sapienza C (2003) Natural selection and the evolution of genome imprinting. Annu Rev Genet 37: 349–370.
  70. 70. Crouse HV (1960) The controlling element in sex chromosome behavior in Sciara. Genetics 45: 1429–1443.
  71. 71. Brun LO, Stuart J, Gaudichon V, Aronstein K, French-Constant RH (1995) Functional haplodiploidy: a mechanism for the spread of insecticide resistance in an important international insect pest. Proc Natl Acad Sci USA 92: 9861–9865.
  72. 72. Mable BK, Bogart JP (1995) Hybridization between tetraploid and diploid species of treefrogs (Genus Hyla). J Heredity 86: 432–440.
  73. 73. Dobson SL, Tanouye MA (1998) Evidence for a genomic imprinting sex determination mechanism in Nasonia vitripennis (Hymenoptera: Chalcidoidea). Genetics 149: 233–242.
  74. 74. Goday C, Esteban MR (2001) Chromsome elimination in sciarid flies. BioEssays 23: 242–250.
  75. 75. Vinogradov AE, Borkin LJ, Gunther R, Rosanov JM (1990) Genome elimination in diploid and triploid Rana esculenta males: cytological evidence from DNA flow cytometry. Genome 33: 619–627.
  76. 76. Osborn TC, Pires CJ, Birchler JA, Auger DL, Chen ZJ, et al. (2003) Understanding mechanisms of novel gene expression in polyploids. Trends Genet 19: 141–147.
  77. 77. Adams KL, Wendel JF (2005) Novel patterns of gene expression in polyploidy plants. Trends Genet 21: 539–543.
  78. 78. Birchler JA, Riddle NC, Auger DL, Veitia RA (2005) Dosage balance in gene regulation: biological implications. Trends Genet 21: 219–226.
  79. 79. Chen ZJ, Ni Z (2006) Mechanisms of genomic rearrangements and gene expression changes in plant polyploids. BioEssays 28: 240–252.
  80. 80. Blackler AW, Gecking CA (1972) Transmission of sex cells of one species through the body of a second species in the genus Xenopus. II. Interspecific matings. Dev Biol 27: 385–394.
  81. 81. Brown DD, Dawid IB, Reeder RH (1977) Xenopus borealis misidentified as Xenopus mulleri. Dev Biol 59: 266–267.
  82. 82. Honjo T, Reeder RH (1973) Preferential transcription of Xenopus laevis ribosomal RNA in interspecies hybrids between Xenopus laevis and Xenopus mulleri. J Mol Biol 80: 217–228.
  83. 83. Navashin M (1934) Chromosomal alterations caused by hybridization and their bearing upon certain general genetic problems. Cytologia 5: 169–203.
  84. 84. Wallace H, Langridge WHR (1971) Differential amphiplasty and the control of ribosomal RNA synthesis. Heredity 27: 1–13.
  85. 85. Wilkinson J (1944) The cytology of Salix in relation to its taxonomy. Ann Bot 8: 269–284.
  86. 86. Keep E (1960) Amphiplasty in Ribes. Nature 188: 339.
  87. 87. Keep E (1962) Satellite and nucleolar number in hybrids between Ribes nigrum and R. grossulari and in their backcross. Can J Genet Cytol 4: 206–218.
  88. 88. Yeh BP, Peloquin SJ (1965) The nucleolus associated chromosome of Solanum species and hybrids. Am J Bot 52: 626.
  89. 89. Kasha KJ, Sadasivaiah RS (1971) Genome relationships between Hordeum vulgare L. and H. bulbosum L. Chromosoma 35: 264–287.
  90. 90. Lange W, Jochemsen G (1976) Karyotypes, nucleoli, and amphiplasty in hybrids between Hordeum vulgare L. and Hordeum bulbosum L. Genetica 46: 217–233.
  91. 91. Nicoloff H (1979) ‘Nucelolar dominance’ as observed in barley translocation lines with specifically reconstructed SAT chromosomes. Theor Appl Genet 55: 247–251.
  92. 92. Schubert I, Kunzel G (1990) Position-dependent NOR activity in barley. Chromosoma 99: 352–359.
  93. 93. Crosby AR (1957) Nucleolar activity of lagging chromosomes in wheat. Am J Bot 44: 813–822.
  94. 94. Flavell RB, O'Dell M, Thompson WF (1988) Regulation of cytosine methylation in ribosomal DNA and nucleolus organizer expression in wheat. J Mol Biol 204: 523–534.
  95. 95. Thomposon WF, Flavell RB (1988) DNase I sensitivity of ribosomal RNA genes in chromatin and nucleolar dominance in wheat. J Mol Biol 204: 535–548.
  96. 96. Sardana R, O'Dell M, Flavell RB (1993) Correlation between the size of the intergenic regulatory region, the status of cytosine methylation of rRNA genes and nucleoloar expression in wheat. Mol Gen Genet 236: 155–162.
  97. 97. Heneen WK (1962) Karyotype studies in Agropyron junceum, A. repens and their spontaneous hybrids. Hereditas 48: 471–502.
  98. 98. Chen ZJ, Pikaard CS (1997) Epigenetic silencing of RNA polymerase I transcription: a role for DNA methylation and histone modification in nucleolar dominance. Genes Dev 11: 2124–2136.
  99. 99. Chen ZJ, Pikaard CS (1997) Transcriptional analysis of nucleolar dominance in polyploidy plants: biased expression/silencing of progenitor rRNA genes is developmentally regulated in Brassica. Proc Natl Acad Sci USA 94: 3442–3447.
  100. 100. Chen ZJ, Comai L, Pikaard CS (1998) Gene dosage and stochastic effects determine the severity and direction of uniparental rRNA gene silencing (nucleolar dominance) in Arabidopsis allopolyploids. Proc Natl Acad Sci USA 95: 14891–14896.
  101. 101. Durica DS, Krider HM (1977) Studies on the ribosomal RNA cistrons in interspecific Drosophila hybrids. I. Nucleolar dominance. Dev Biol 59: 62–74.
  102. 102. Durica DS, Krider HM (1978) Studies on the ribosomal RNA cistrons in Drosophila hybrids. II. Heterochromatic regions mediating nucleolar dominance. Genetics 89: 37–64.
  103. 103. Elicieri GL, Green H (1969) Ribosomal RNA synthesis in human-mouse hybrid cells. J Mol Biol 41: 253–260.
  104. 104. Bramwell ME, Handmaker SD (1971) Ribosomal RNA synthesis in human-mouse hybrid cells. Biochim Biophys Acta 232: 580–593.
  105. 105. Miller OJ, Miller DA, Dev VG, Tantravahi R, Croce CM (1976) Expression of human and suppression of mouse nucleolus organizer activity in mouse-human somatic cell hybrids. Proc Natl Acad Sci USA 73: 4531–4535.
  106. 106. Croce CM, Talavera A, Basilico C, Miller OJ (1977) Suppression of production of mouse 28S rRNA in mouse-human hybrids segregating chromosomes. Proc Natl Acad Sci USA 74: 694–697.
  107. 107. Soprano KJ, Dev VG, Croce CM, Baserga R (1979) Reactivation of silent rRNA genes by simian virus 40 in mouse-human hybrid cells. Proc Natl Acad Sci USA 76: 3885–3889.
  108. 108. Soprano KJ, Baserga R (1980) Reactivation of ribosomal RNA genes in human-mouse hybrid cells by 12-O-tetradecanoylphorbol 13-acetate. Proc Natl Acad Sci USA 77: 1566–1569.
  109. 109. Onishi T, Berglund C, Reeder RH (1984) On the mechanism of nucleolar dominance in mouse-human somatic cell hybrids. Proc Natl Acad Sci USA 81: 484–487.
  110. 110. Heard E, Clerc P, Avner P (1997) X-chromosome inactivation in mammals. Annu Rev Genet 31: 571–610.
  111. 111. Reeder RH (1985) Mechanisms of nucleolar dominance in animals and plants. J Cell Biol 101: 2013–2016.
  112. 112. Pikaard CS (2000) Nucleolar dominance: uniparental gene silencing on a multi-megabase scale in genetic hybrids. Plant Molecular Biology 43: 163–177.
  113. 113. Viegas W, Neves N, Silva M, Caperta A, Morais-Cecílio L (2002) Nucleolar dominance: A ‘David and Goliath’ chromatin imprinting process. Curr Genomics 3: 563–576.
  114. 114. Appels R, Dvorak J (1982) The wheat ribosomal DNA spacer region: its structure and variation in populations and amongst species. Theor Appl Genet 63: 337–348.
  115. 115. Dover GA, Flavell RB (1984) Molecular co-evolution: rDNA divergence and the maintenance of function. Cell 38: 622–623.
  116. 116. Grummt I, Roth E, Paule MR (1982) Ribosomal RNA transcription is species specific. Nature 296: 173–174.
  117. 117. Miesfeld R, Arnheim N (1984) Species-specific rDNA transcription is due to promoter-specific binding factors. Mol Cell Biol 4: 221–227.
  118. 118. Miesfeld R, Solner-Webb B, Croce C, Arnheim N (1984) The absence of human-specific ribosomal DNA transcription factor leads to nucleolar dominance in mouse-human hybrid cells. Mol Cell Biol 4: 1306–1312.
  119. 119. Schnapp A, Rosenbauer H, Grumm I (1991) Trans-acting factors involved in species-specificity and control of mouse ribosomal gene transcription. Mol Cell Biochem 104: 137–147.
  120. 120. Fan H, Yakura K, Miyanishi M, Sugita M, Sugiura M (1995) Transcription of plant RNA polymerase I-dependent rRNA genes is species-specific. Plant J 8: 295–298.
  121. 121. Doelling JH, Pikaard CS (1996) Species-specificity of rRNA gene transcription in plants manifested as a switch in polymerase specificity. Nucleic Acids Res 24: 4725–4732.
  122. 122. Frieman M, Chen ZJ, Saez-Vasquez J, Shen LA, Pikaard CS (1999) RNA polymerase I transcription in a Brassica interspecific hybrid and its progenitors: tests of transcription factor involvement in nucleolar dominance. Genetics 152: 451–460.
  123. 123. Chain FJJ, Evans BJ (2006) Multiple mechanisms promote the retained expression of gene duplicates in the tetraploid frog Xenopus laevis. PLoS Genetics 2: e56.
  124. 124. Gehring M, Huh JH, Tzung-Fu H, Peneterman J, Yeonhee C, et al. (2006) DEMETER DNA glycosylase establisheds MEDEA polycomb gene self-imprinting by allele-specific demethylation. Cell 124: 495–506.
  125. 125. Busby SJ, Reeder RH (1983) Spacer sequences regulate transcription of ribosomal gene plasmids infected in Xenopus embryos. Cell 34: 989–996.
  126. 126. Reeder RH, Roan JG, Dunaway M (1983) Spacer regulation of Xenopus ribosomal gene transcription: competition in oocytes. Cell 35: 449–456.
  127. 127. Labhart P, Reeder RH (1984) Enhancer-like properties of the 60/81 bp elements in the ribosomal gene spacer of Xenopus laevis. Cell 37: 285–289.
  128. 128. Pikaard CS, Reeder RH (1988) Sequence elements essential for function of the Xenopus laevis ribosomal DNA enhancers. Mol Cell Biol 8: 4282–4288.
  129. 129. Boseley P, Moss T, Machler M, Portmann R, Birnstiel M (1979) Sequence of the spacer DNA in a ribosomal gene unit of X. laevis. Cell 17: 19–31.
  130. 130. Sollner-Webb B, Reeder RH (1979) The nucleotide sequence of the initiation and termination sites for ribosomal RNA transcription in X. laevis. Cell 18: 485–499.
  131. 131. Moss T (1983) A transcriptional function for the repetitive ribosomal spacer in Xenopus laevis. Nature 302: 223–228.
  132. 132. Cassidy DM, Blacker AW (1974) Repression of nucleloar organizer activity in an interspecific hybrid of the genus Xenopus. Dev Biol 41: 84–96.
  133. 133. Reeder RH (1984) Enhancers and ribosomal gene spacers. Cell 38: 349–351.
  134. 134. Reeder RH, Roan JG (1984) The mechanism of nucleolar dominance in Xenopus hybrids. Cell 38: 39–44.
  135. 135. Porter AH, Johnson NA (2002) Speciation despite gene flow when developmental pathways evolve. Evol Int J Org Evol 56: 2103–2111.
  136. 136. Landry CR, Wittkopp PJ, Taubes CH, Ranz JM, Clark AG, et al. (2005) Compensatory cis-trans evolution and the dysregulation of gene expression in interspecific hybrids of Drosophila. Genetics 171: 1813–1822.
  137. 137. Shilatifard A (2006) Chromatin modifications by methylation and ubiquitnation: implications in the regulation of gene expression. Annu Rev Biochem 75: 243–269.
  138. 138. Jenuwein T, Allis CD (2001) Translating the histone code. Science 293: 1074–1080.
  139. 139. McStay B (2006) Nucleolar dominance: a model for rRNA gene silencing. Genes and Dev 20: 1207–1214.
  140. 140. Richards EJ, Elgin SCR (2002) Epigenetic codes for heterochromatin formation and silencing, rounding up the usual suspects. Cell 108: 489–500.
  141. 141. Malone JH, Hawkins DL Jr, Michalak P (2006) Sex-biased gene expression in a ZW sex determination system. J Mol Evol 63: 427–436.
  142. 142. Presgraves DC (2002) Patterns of postzygotic isolation in Lepidoptera. Evolution 56: 1168–1183.
  143. 143. Presgraves DC, Orr HA (1998) Haldane's rule in taxa lacking a hemizygous X. Science 282: 952–954.
  144. 144. Presgraves DC, Balagopalan L, Abmayr SM, Orr HA (2003) Adaptive evolution drives divergence of a hybrid inviability gene between two species of Drosophila. Nature 423: 715–719.
  145. 145. Kelley DB (1996) Sexual differentiation in Xenopus laevis. In: Tinsley RC, Kobel HR, editors. The biology of Xenopus. Oxford: Clarendon Press. pp. 143–176.
  146. 146. Hammond JP, Broadley MR, Craigon DJ, Higgins J, Emmerson ZF, et al. (2005) Using genomic DNA-based pobe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species. Plant Meth 1: 10.
  147. 147. Bolstad B (2005) RMAExpress. RMAExpress.html.
  148. 148. Smyth G (2004) Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology 3: Article 3.