A new method to remove hybridization bias for interspecies comparison of global gene expression profiles uncovers an association between mRNA sequence divergence and differential gene expression in Xenopus

The recent sequencing of a large number of Xenopus tropicalis expressed sequences has allowed development of a high-throughput approach to study Xenopus global RNA gene expression. We examined the global gene expression similarities and differences between the historically significant Xenopus laevis model system and the increasingly used X.tropicalis model system and assessed whether an X.tropicalis microarray platform can be used for X.laevis. These closely related species were also used to investigate a more general question: is there an association between mRNA sequence divergence and differences in gene expression levels? We carried out a comprehensive comparison of global gene expression profiles using microarrays of different tissues and developmental stages of X.laevis and X.tropicalis. We (i) show that the X.tropicalis probes provide an efficacious microarray platform for X.laevis, (ii) describe methods to compare interspecies mRNA profiles that correct differences in hybridization efficiency and (iii) show independently of hybridization bias that as mRNA sequence divergence increases between X.laevis and X.tropicalis differences in mRNA expression levels also increase.     

Identification of novel microRNAs in Xenopus laevis metaphase II arrested eggs

Using a combination of deep sequencing and bioinformatics approach, we for the first time identify miRNAs and their relative abundance in mature, metaphase II arrested eggs in Xenopus laevis. We characterize 115 miRNAs that have been described either in Xenopus tropicalis (85), X. laevis (9), or other vertebrate species (21) that also map to known Xenopus pre-miRNAs and to the X. tropicalis genome. In addition, 72 new X. laevis putative candidate miRNAs are identified based on mapping to X. tropicalis genome within regions that have the propensity to form hairpin loops. These data expand on the availability of genetic information in X. laevis and identify target miRNAs for future functional studies.     

One of the duplicated matrix metalloproteinase-9 genes is expressed in regressing tail during anuran metamorphosis

The drastic morphological changes of the tadpole are induced during the climax of anuran metamorphosis, when the concentration of endogenous thyroid hormone is maximal. The tadpole tail, which is twice as long as the body, shortens rapidly and disappears completely in several days. We isolated a cDNA clone, designated as Xl MMP-9TH, similar to the previously reported Xenopus laevis MMP-9 gene, and showed that their Xenopus tropicalis counterparts are located tandemly about 9 kb apart from each other in the genome. The Xenopus MMP-9TH gene was expressed in the regressing tail and gills and the remodeling intestine and central nervous system, and induced in thyroid hormone-treated tail-derived myoblastic cultured cells, while MMP-9 mRNA was detected in embryos. Three thyroid hormone response elements in the distal promoter and the first intron were involved in the upregulation of the Xl MMP-9TH gene by thyroid hormone in transient expression assays, and their relative positions are conserved between X. laevis and X. tropicalis promoters. These data strongly suggest that the MMP-9 gene was duplicated, and differentiated into two genes, one of which was specialized in a common ancestor of X. laevis and X. tropicalis to be expressed in degenerating and remodeling organs as a response to thyroid hormone during metamorphosis.     

A Xenopus tropicalis oligonucleotide microarray works across species using RNA from Xenopus laevis

Microarrays have great potential for the study of developmental biology. As a model system Xenopus is well suited for making the most of this potential. However, Xenopus laevis has undergone a genome wide duplication meaning that most genes are represented by two paralogues. This causes a number of problems. Most importantly the presence of duplicated genes mean that a X. laevis microarray will have less or even half the coverage of a similar sized microarray from the closely related but diploid frog Xenopus tropicalis. However, to date, X. laevis is the most commonly used amphibian system for experimental embryology. Therefore, we have tested if a microarray based on sequences from X. tropicalis will work across species using RNA from X. laevis. We produced a pilot oligonucleotide microarray based on sequences from X. tropicalis. The microarray was used to identify genes whose expression levels changed during early X. tropicalis development. The same assay was then carried out using RNA from X. laevis. The cross species experiments gave similar results to those using X. tropicalis RNA. This was true at the whole microarray level and for individual genes, with most genes giving similar results using RNA from X. laevis and X. tropicalis. Furthermore, the overlap in genes identified between a X. laevis and a X. tropicalis set of experiments was only 12% less than the overlap between two sets of X. tropicalis experiments. Therefore researchers can work with X. laevis and still make use of the advantages offered by X. tropicalis microarrays.     

Regulated expression of the X. tropicalis connexin43 promoter

The spatio-temporal expression pattern of the connexin43 gene during Xenopus development has been described (Van der Heyden et al. 2001). To further investigate the regulation and function of connexin43 (Cx43) in amphibians, we have isolated the gene from Xenopus tropicalis (Xt) and determined its structure. The X. tropicalis Cx43 gene displays the typical two exon-one intron connexin configuration, where the first exon is non-coding. The predicted amino acid sequence of the XtCx43 protein is highly homologous to that of X. laevis, chicken and mammals. Expression of XtCx43 cDNA in N2A cells results in gap-junction plaque formation. Promoter activity of a 3.5 kb upstream region of the X. tropicalis Cx43 gene, including exon 1, mimics endogenous timing of expression after injection of reporter constructs in X. laevis embryos.     

Isolation and characterization of sarcomeric actin genes expressed in Xenopus laevis embryos

A Xenopus laevis complementary DNA (cDNA) library prepared from messenger RNAs extracted from embryos has been screened for actin-coding sequences. Two cDNA clones corresponding to an alpha cardiac and an alpha skeletal muscle actin mRNA have been identified and characterized. From a genomic library, we have furthermore isolated the genes that correspond to the characterized cDNAs. In addition we have identified an actin processed gene which seems to be derived from a second type of skeletal muscle actin gene. Southern blot analysis of X. laevis DNA reveals that each of the three genes is present in at least two copies. In Xenopus tropicalis, a similar Southern blot analysis demonstrates that the three alpha actin genes exist as single copy. This result correlates with the genome duplication that has been proposed to have occurred recently in a X. laevis ancestor. A sequence comparison of the X. laevis cardiac and skeletal muscle actin cDNAs shows that the encoded peptides are highly conserved. Nevertheless, the numerous nucleotide changes at silent mutation sites suggest that the genes originated before the amphibia/reptile-bird divergence, more than 350 million years ago. Comparison of the promoters of the cardiac and skeletal actin genes, which are co-expressed in embryos, reveals a few common structural sequence elements.     

Vitellogenin genes and their products in closely and distantly related species of Xenopus

Plasma vitellogenins from two closely related species of Xenopus, X. laevis and X. borealis, and a more ancient species, X. tropicalis, exhibited the same size on gel electrophoresis and were immunologically related. Partial peptide maps of 125I-labelled plasma vitellogenins, however, revealed marked differences in th structure and organisation of vitellogenin in the three Xenopus species. Northern blot hybridisation of liver RNA from oestrogen-treated males and females, probed with cloned vitellogenin cDNA, revealed the presence of mRNA of the same size in the three species of Xenopus, which was absent in untreated male liver. Cell-free translation of total liver RNA showed the presence of functional mRNA coding for vitellogenin subunit of the same size (Mr congruent to 210,000). Restriction endonuclease digestion patterns of genomic DNA from the three Xenopus species, using cloned X. laevis vitellogenin cDNA as the hybridisation probe, revealed significant differences in the organisation of these genes, which occur at a higher multiplicity in X. laevis and X. borealis than in X. tropicalis. Thus, despite a high degree of conservation of size, overall sequence and immunological identity of vitellogenin genes and their products in the three species of Xenopus, there is a substantial structural rearrangement during evolution of Xenopus within this multigene family.