From expression cloning to gene modeling: the development of Xenopus gene sequence resources

Using Xenopus to discover new candidate genes involved in BOR and other congenital hearing loss syndromes

Hearing in infants is essential for brain development, acquisition of verbal language skills, and development of social interactions. Therefore, it is important to diagnose hearing loss soon after birth so that interventions can be provided as early as possible. Most newborns in the United States are screened for hearing deficits and commercially available next-generation sequencing hearing loss panels often can identify the causative gene, which may also identify congenital defects in other organs. One of the most prevalent autosomal dominant congenital hearing loss syndromes is branchio-oto-renal syndrome (BOR), which also presents with defects in craniofacial structures and the kidney. Currently, mutations in three genes, SIX1, SIX5, and EYA1, are known to be causative in about half of the BOR patients that have been tested. To uncover new candidate genes that could be added to congenital hearing loss genetic screens, we have combined the power of Drosophila mutants and protein biochemical assays with the embryological advantages of Xenopus, a key aquatic animal model with a high level of genomic similarity to human, to identify potential Six1 transcriptional targets and interacting proteins that play a role during otic development. We review our transcriptomic, yeast 2-hybrid, and proteomic approaches that have revealed a large number of new candidates. We also discuss how we have begun to identify how Six1 and co-factors interact to direct developmental events necessary for normal otic development.

Genome evolution in the allotetraploid frog Xenopus laevis

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We demonstrate the allotetraploid origin of X. laevis by partitioning its genome into two homeologous subgenomes, marked by distinct families of “fossil” transposable elements. Based on the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged ~34 million years ago (Mya) and combined to form an allotetraploid ~17–18 Mya. 56% of all genes are retained in two homeologous copies. Protein function, gene expression, and the amount of flanking conserved sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Accurate Profiling of Gene Expression and Alternative Polyadenylation with Whole Transcriptome Termini Site Sequencing (WTTS-Seq)

Construction of next-generation sequencing (NGS) libraries involves RNA manipulation, which often creates noisy, biased, and artifactual data that contribute to errors in transcriptome analysis. In this study, a total of 19 whole transcriptome termini site sequencing (WTTS-seq) and seven RNA sequencing (RNA-seq) libraries were prepared from Xenopus tropicalis adult and embryo samples to determine the most effective library preparation method to maximize transcriptomics investigation. We strongly suggest that appropriate primers/adaptors are designed to inhibit amplification detours and that PCR overamplification is minimized to maximize transcriptome coverage. Furthermore, genome annotation must be improved so that missing data can be recovered. In addition, a complete understanding of sequencing platforms is critical to limit the formation of false-positive results. Technically, the WTTS-seq method enriches both poly(A)+ RNA and complementary DNA, adds 5′- and 3′-adaptors in one step, pursues strand sequencing and mapping, and profiles both gene expression and alternative polyadenylation (APA). Although RNA-seq is cost prohibitive, tends to produce false-positive results, and fails to detect APA diversity and dynamics, its combination with WTTS-seq is necessary to validate transcriptome-wide APA.

Towards the bridging of molecular genetics data across Xenopus species

BACKGROUND

The clawed African frog Xenopus laevis has been one of the main vertebrate models for studies in developmental biology. However, for genetic studies, Xenopus tropicalis has been the experimental model of choice because it shorter life cycle and due to a more tractable genome that does not result from genome duplication as in the case of X. laevis. Today, although still organized in a large number of scaffolds, nearly 85% of X. tropicalis and 89% of X. laevis genomes have been sequenced. There is expectation for a comparative physical map that can be used as a Rosetta Stone between X. laevis genetic studies and X. tropicalis genomic research.

RESULTS

In this work, we have mapped using coarse-grained alignment the 18 chromosomes of X. laevis, release 9.1, on the 10 reference scaffolds representing the haploid genome of X. tropicalis, release 9.0. After validating the mapping with theoretical data, and estimating reference averages of genome sequence identity, 37 to 44% between the two species, we have carried out a synteny analysis for 2,112 orthologous genes. We found that 99.6% of genes are in the same organization.

CONCLUSIONS

Taken together, our results make possible to establish the correspondence between 62 and 65.5% of both genomes, percentage of identity, synteny and automatic annotation of transcripts of both species, providing a new and more comprehensive tool for comparative analysis of these two species, by allowing to bridge molecular genetics data among them.

Using Xenopus to discover new genes involved in branchiootorenal spectrum disorders

Congenital hearing loss is an important clinical problem because, without early intervention, affected children do not properly acquire language and consequently have difficulties developing social skills. Although most newborns in the US are screened for hearing deficits, even earlier diagnosis can be made with prenatal genetic screening. Genetic screening that identifies the relevant mutated gene can also warn about potential congenital defects in organs not related to hearing. We will discuss efforts to identify new candidate genes that underlie the Branchiootorenal spectrum disorders in which affected children have hearing deficits and are also at risk for kidney defects. Mutations in two genes, SIX1 and EYA1, have been identified in about half of the patients tested. To uncover new candidate genes, we have used the aquatic animal model, Xenopus laevis, to identify genes that are part of the developmental genetic pathway of Six1 during otic and kidney development. We have already identified a large number of potential Six1 transcriptional targets and candidate co-factor proteins that are expressed at the right time and in the correct tissues to interact with Six1 during development. We discuss the advantages of using this system for gene discovery in a human congenital hearing loss syndrome.

The Xenopus ORFeome: A resource that enables functional genomics

Functional characterisation of proteins and large-scale, systems-level studies are enabled by extensive sets of cloned open reading frames (ORFs) in an easily-accessible format that enables many different applications. Here we report the release of the first stage of the Xenopus ORFeome, which contains 8673 ORFs from the Xenopus Gene Collection (XGC) for Xenopus laevis, cloned into a Gateway® donor vector enabling rapid in-frame transfer of the ORFs to expression vectors. This resource represents an estimated 7871 unique genes, approximately 40% of the non-redundant X. laevis gene complement, and includes 2724 genes where the human ortholog has an association with disease. Transfer into the Gateway system was validated by 5' and 3' end sequencing of the entire collection and protein expression of a set of test clones. In a parallel process, the underlying ORF predictions from the original XGC collection were re-analysed to verify quality and full-length status, identifying those proteins likely to exhibit truncations when translated. These data are integrated into Xenbase, the Xenopus community database, which associates genomic, expression, function and human disease model metadata to each ORF, enabling end-users to search for ORFeome clones with links to commercial distributors of the collection. When coupled with the experimental advantages of Xenopus eggs and embryos, the ORFeome collection represents a valuable resource for functional genomics and disease modelling.

Coexistence of Y, W, and Z sex chromosomes in Xenopus tropicalis

Homomorphic sex chromosomes and rapid turnover of sex-determining genes can complicate establishing the sex chromosome system operating in a given species. This difficulty exists in Xenopus tropicalis, an anuran quickly becoming a relevant model for genetic, genomic, biochemical, and ecotoxicological research. Despite the recent interest attracted by this species, little is known about its sex chromosome system. Direct evidence that females are the heterogametic sex, as in the related species Xenopus laevis, has yet to be presented. Furthermore, X. laevis' sex-determining gene, DM-W, does not exist in X. tropicalis, and the sex chromosomes in the two species are not homologous. Here we identify X. tropicalis' sex chromosome system by integrating data from (i) breeding sex-reversed individuals, (ii) gynogenesis, (iii) triploids, and (iv) crosses among several strains. Our results indicate that at least three different types of sex chromosomes exist: Y, W, and Z, observed in YZ, YW, and ZZ males and in ZW and WW females. Because some combinations of parental sex chromosomes produce unisex offspring and other distorted sex ratios, understanding the sex-determination systems in X. tropicalis is critical for developing this flexible animal model for genetics and ecotoxicology.

Microscale imaging of cilia-driven fluid flow

Cilia-driven fluid flow is important for multiple processes in the body, including respiratory mucus clearance, gamete transport in the oviduct, right-left patterning in the embryonic node, and cerebrospinal fluid circulation. Multiple imaging techniques have been applied toward quantifying ciliary flow. Here, we review common velocimetry methods of quantifying fluid flow. We then discuss four important optical modalities, including light microscopy, epifluorescence, confocal microscopy, and optical coherence tomography, that have been used to investigate cilia-driven flow.

Comparative genomic analysis of slc39a12/ZIP12: insight into a zinc transporter required for vertebrate nervous system development

The zinc transporter ZIP12, which is encoded by the gene slc39a12, has previously been shown to be important for neuronal differentiation in mouse Neuro-2a neuroblastoma cells and primary mouse neurons and necessary for neurulation during Xenopus tropicalis embryogenesis. However, relatively little is known about the biochemical properties, cellular regulation, or the physiological role of this gene. The hypothesis that ZIP12 is a zinc transporter important for nervous system function and development guided a comparative genetics approach to uncover the presence of ZIP12 in various genomes and identify conserved sequences and expression patterns associated with ZIP12. Ortholog detection of slc39a12 was conducted with reciprocal BLAST hits with the amino acid sequence of human ZIP12 in comparison to the human paralog ZIP4 and conserved local synteny between genomes. ZIP12 is present in the genomes of almost all vertebrates examined, from humans and other mammals to most teleost fish. However, ZIP12 appears to be absent from the zebrafish genome. The discrimination of ZIP12 compared to ZIP4 was unsuccessful or inconclusive in other invertebrate chordates and deuterostomes. Splice variation, due to the inclusion or exclusion of a conserved exon, is present in humans, rats, and cows and likely has biological significance. ZIP12 also possesses many putative di-leucine and tyrosine motifs often associated with intracellular trafficking, which may control cellular zinc uptake activity through the localization of ZIP12 within the cell. These findings highlight multiple aspects of ZIP12 at the biochemical, cellular, and physiological levels with likely biological significance. ZIP12 appears to have conserved function as a zinc uptake transporter in vertebrate nervous system development. Consequently, the role of ZIP12 may be an important link to reported congenital malformations in numerous animal models and humans that are caused by zinc deficiency.

Engineering Xenopus embryos for phenotypic drug discovery screening

Many rare human inherited diseases remain untreatable despite the fact that the disease causing genes are known and adequate mouse disease models have been developed. In vivo phenotypic drug screening relies on isolating drug candidates by their ability to produce a desired therapeutic phenotype in whole organisms. Embryos of zebrafish and Xenopus frogs are abundant, small and free-living. They can be easily arrayed in multi-well dishes and treated with small organic molecules. With the development of novel genome modification tools, such a zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR/Cas, it is now possible to efficiently engineer non-mammalian models of inherited human diseases. Here, we will review the rapid progress made in adapting these novel genome editing tools to Xenopus. The advantages of Xenopus embryos as in vivo models to study human inherited diseases will be presented and their utility for drug discovery screening will be discussed. Being a tetrapod, Xenopus complements zebrafish as an indispensable non-mammalian animal model for the study of human disease pathologies and the discovery of novel therapeutics for inherited diseases.

Combining different mRNA capture methods to analyze the transcriptome: analysis of the Xenopus laevis transcriptome

mRNA sequencing (mRNA-seq) is a commonly used technique to survey gene expression from organisms with fully sequenced genomes. Successful mRNA-seq requires purification of mRNA away from the much more abundant ribosomal RNA, which is typically accomplished by oligo-dT selection. However, mRNAs with short poly-A tails are captured poorly by oligo-dT based methods. We demonstrate that combining mRNA capture via oligo-dT with mRNA capture by the 5’ 7-methyl guanosine cap provides a more complete view of the transcriptome and can be used to assay changes in mRNA poly-A tail length on a genome-wide scale. We also show that using mRNA-seq reads from both capture methods as input for de novo assemblers provides a more complete reconstruction of the transcriptome than either method used alone. We apply these methods of mRNA capture and de novo assembly to the transcriptome of Xenopus laevis, a well-studied frog that currently lacks a finished sequenced genome, to discover transcript sequences for thousands of mRNAs that are currently absent from public databases. The methods we describe here will be broadly applicable to many organisms and will provide insight into the transcriptomes of organisms with sequenced and unsequenced genomes.

Efficient high-throughput sequencing of a laser microdissected chromosome arm

BACKGROUND

Genomic sequence assemblies are key tools for a broad range of gene function and evolutionary studies. The diploid amphibian Xenopus tropicalis plays a pivotal role in these fields due to its combination of experimental flexibility, diploid genome, and early-branching tetrapod taxonomic position, having diverged from the amniote lineage ~360 million years ago. A genome assembly and a genetic linkage map have recently been made available. Unfortunately, large gaps in the linkage map attenuate long-range integrity of the genome assembly.

RESULTS

We laser dissected the short arm of X. tropicalis chromosome 7 for next generation sequencing and computational mapping to the reference genome. This arm is of particular interest as it encodes the sex determination locus, but its genetic map contains large gaps which undermine available genome assemblies. Whole genome amplification of 15 laser-microdissected 7p arms followed by next generation sequencing yielded ~35 million reads, over four million of which uniquely mapped to the X. tropicalis genome. Our analysis placed more than 200 previously unmapped scaffolds on the analyzed chromosome arm, providing valuable low-resolution physical map information for de novo genome assembly.

CONCLUSION

We present a new approach for improving and validating genetic maps and sequence assemblies. Whole genome amplification of 15 microdissected chromosome arms provided sufficient high-quality material for localizing previously unmapped scaffolds and genes as well as recognizing mislocalized scaffolds.

An intact brachyury function is necessary to prevent spurious axial development in Xenopus laevis

We have previously shown that the member of the HES family hairy2 induces the ectopic expression of dorsal markers when it is overexpressed in the ventral side of Xenopus embryos. Intriguingly, hairy2 represses the mesoderm transcription factor brachyury (bra) throughout its domain in the marginal zone. Here we show that in early gastrula, bra and hairy2 are expressed in complementary domains. Overexpression of bra repressed hairy2. Interference of bra function with a dominant-negative construct expanded the hairy2 domain and, like hairy2 overexpression, promoted ectopic expression of dorsal axial markers in the ventral side and induced secondary axes without head and notochord. Hairy2 depletion rescued the ectopic dorsal development induced by interference of bra function. We concluded that an intact bra function is necessary to exclude hairy2 expression from the non-organiser field, to impede the ectopic specification of dorsal axial fates and the appearance of incomplete secondary axes. This evidence supports a previously unrecognised role for bra in maintaining the dorsal fates inhibited in the ventral marginal zone, preventing the appearance of trunk duplications.

High-throughput sequencing will metamorphose the analysis of thyroid hormone receptor function during amphibian development

Amphibian metamorphosis is marked by dramatic thyroid hormone (T(3))-induced changes including de novo morphogenesis, tissue remodeling, and organ resorption through programmed cell death. These changes involve cascades of gene regulation initiated by thyroid hormone (TH). TH functions by regulating gene expression through TH receptors (TR). TR are DNA-binding transcription factors that belong to the steroid hormone receptor superfamily. In the absence of ligand, TR can repress gene expression by recruiting a corepressor complex, whereas liganded TR recruits a coactivator complex for gene activation. Earlier studies have led us to propose a dual function model for TR during development. In premetamorphic tadpoles, unliganded TR represses transcription involving corepressors. During metamorphosis, endogenous T(3) allows TR to activate gene expression. To fully understand the diversity of T(3) effects during metamorphosis, whole genome analysis of transcriptome and mechanism of TR action should be carried out. To this end, the new sequencing technologies have dramatically changed how fundamental questions in biology are being addressed and is now making the transition from technology development to being a standard for genomic and functional genomic analysis. This review focuses on the applications of high-throughput technologies to the field of amphibian metamorphosis.

Exon capture and bulk segregant analysis: rapid discovery of causative mutations using high-throughput sequencing

Background

Exome sequencing has transformed human genetic analysis and may do the same for other vertebrate model systems. However, a major challenge is sifting through the large number of sequence variants to identify the causative mutation for a given phenotype. In models like Xenopus tropicalis, an incomplete and occasionally incorrect genome assembly compounds this problem. To facilitate cloning of X. tropicalis mutants identified in forward genetic screens, we sought to combine bulk segregant analysis and exome sequencing into a single step.

Results

Here we report the first use of exon capture sequencing to identify mutations in a non-mammalian, vertebrate model. We demonstrate that bulk segregant analysis coupled with exon capture sequencing is not only able to identify causative mutations but can also generate linkage information, facilitate the assembly of scaffolds, identify misassembles, and discover thousands of SNPs for fine mapping.

Conclusion

Exon capture sequencing and bulk segregant analysis is a rapid, inexpensive method to clone mutants identified in forward genetic screens. With sufficient meioses, this method can be generalized to any model system with a genome assembly, polished or unpolished, and in the latter case, it also provides many critical genomic resources.