Mitochondrial leak metabolism induces the Spemann-Mangold Organizer via Hif-1α in Xenopus.

An instructive role for metabolism in embryonic patterning is emerging, although a role for mitochondria is poorly defined. We demonstrate that mitochondrial oxidative metabolism establishes the embryonic patterning center, the Spemann-Mangold Organizer, via hypoxia-inducible factor 1α (Hif-1α) in Xenopus. Hypoxia or decoupling ATP production from oxygen consumption expands the Organizer by activating Hif-1α. In addition, oxygen consumption is 20% higher in the Organizer than in the ventral mesoderm, indicating an elevation in mitochondrial respiration. To reconcile increased mitochondrial respiration with activation of Hif-1α, we discovered that the "free" c-subunit ring of the F1Fo ATP synthase creates an inner mitochondrial membrane leak, which decouples ATP production from respiration at the Organizer, driving Hif-1α activation there. Overexpression of either the c-subunit or Hif-1α is sufficient to induce Organizer cell fates even when ß-catenin is inhibited. We propose that mitochondrial leak metabolism could be a general mechanism for activating Hif-1α and Wnt signaling.

SMC5 Plays Independent Roles in Congenital Heart Disease and Neurodevelopmental Disability.

Up to 50% of patients with severe congenital heart disease (CHD) develop life-altering neurodevelopmental disability (NDD). It has been presumed that NDD arises in CHD cases because of hypoxia before, during, or after cardiac surgery. Recent studies detected an enrichment in de novo mutations in CHD and NDD, as well as significant overlap between CHD and NDD candidate genes. However, there is limited evidence demonstrating that genes causing CHD can produce NDD independent of hypoxia. A patient with hypoplastic left heart syndrome and gross motor delay presented with a de novo mutation in SMC5. Modeling mutation of smc5 in Xenopus tropicalis embryos resulted in reduced heart size, decreased brain length, and disrupted pax6 patterning. To evaluate the cardiac development, we induced the conditional knockout (cKO) of Smc5 in mouse cardiomyocytes, which led to the depletion of mature cardiomyocytes and abnormal contractility. To test a role for Smc5 specifically in the brain, we induced cKO in the mouse central nervous system, which resulted in decreased brain volume, and diminished connectivity between areas related to motor function but did not affect vascular or brain ventricular volume. We propose that genetic factors, rather than hypoxia alone, can contribute when NDD and CHD cases occur concurrently.

Membrane potential drives the exit from pluripotency and cell fate commitment via calcium and mTOR.

Transitioning from pluripotency to differentiated cell fates is fundamental to both embryonic development and adult tissue homeostasis. Improving our understanding of this transition would facilitate our ability to manipulate pluripotent cells into tissues for therapeutic use. Here, we show that membrane voltage (Vm) regulates the exit from pluripotency and the onset of germ layer differentiation in the embryo, a process that affects both gastrulation and left-right patterning. By examining candidate genes of congenital heart disease and heterotaxy, we identify KCNH6, a member of the ether-a-go-go class of potassium channels that hyperpolarizes the Vm and thus limits the activation of voltage gated calcium channels, lowering intracellular calcium. In pluripotent embryonic cells, depletion of kcnh6 leads to membrane depolarization, elevation of intracellular calcium levels, and the maintenance of a pluripotent state at the expense of differentiation into ectodermal and myogenic lineages. Using high-resolution temporal transcriptome analysis, we identify the gene regulatory networks downstream of membrane depolarization and calcium signaling and discover that inhibition of the mTOR pathway transitions the pluripotent cell to a differentiated fate. By manipulating Vm using a suite of tools, we establish a bioelectric pathway that regulates pluripotency in vertebrates, including human embryonic stem cells.

Kap-ß2/Transportin mediates ß-catenin nuclear transport in Wnt signaling.

Wnt signaling is essential for many aspects of embryonic development including the formation of the primary embryonic axis. In addition, excessive Wnt signaling drives multiple diseases including cancer, highlighting its importance for disease pathogenesis. ß-catenin is a key effector in this pathway that translocates into the nucleus and activates Wnt responsive genes. However, due to our lack of understanding of ß-catenin nuclear transport, therapeutic modulation of Wnt signaling has been challenging. Here, we took an unconventional approach to address this long-standing question by exploiting a heterologous model system, the budding yeast Saccharomyces cerevisiae, which contains a conserved nuclear transport machinery. In contrast to prior work, we demonstrate that ß-catenin accumulates in the nucleus in a Ran-dependent manner, suggesting the use of a nuclear transport receptor (NTR). Indeed, a systematic and conditional inhibition of NTRs revealed that only Kap104, the ortholog of Kap-ß2/Transportin-1 (TNPO1), was required for ß-catenin nuclear import. We further demonstrate direct binding between TNPO1 and ß-catenin that is mediated by a conserved PY-NLS. Finally, using Xenopus secondary axis and TCF/LEF (T Cell factor/lymphoid enhancer factor family) reporter assays, we demonstrate that our results in yeast can be directly translated to vertebrates. By elucidating the nuclear localization signal in ß-catenin and its cognate NTR, our study suggests new therapeutic targets for a host of human diseases caused by excessive Wnt signaling. Indeed, we demonstrate that a small chimeric peptide designed to target TNPO1 can reduce Wnt signaling as a first step toward therapeutics.

Obtaining Xenopus tropicalis Embryos by In Vitro Fertilization.

Xenopus is a powerful model system for cell and developmental biology in part because frogs produce thousands of eggs and embryos year-round. In vitro fertilization (IVF) is ideal for obtaining developmentally synchronized embryos for microinjection or when natural mating has failed to produce a fertilization. In IVF, females are induced to ovulate, and then eggs are collected by manual expression. After testes are collected from a euthanized male frog, the eggs are fertilized in vitro. The embryos are then treated with cysteine to remove the sticky protective jelly coat. Dejellied embryos are much easier to manipulate during microinjection or when sorting in a Petri dish. The jelly coat is also very difficult to penetrate with an injection needle. After microinjection, embryos are maintained in Petri dishes until desired stages are reached. Although in vitro fertilization in X. laevis and X. tropicalis is similar, critical differences in solutions, handling of testis, response of fertilized eggs directly after introduction of sperm, and developmental timing are required for successful fertilization in X. tropicalis.

Obtaining Xenopus tropicalis Embryos by Natural Mating.

Xenopus is a powerful model system for cell and developmental biology in part because frogs produce thousands of eggs and embryos year-round. Natural matings are a simple and common method to obtain embryos for injection or other experimental use or to raise to adulthood. This method does not require sacrificing a male as in vitro fertilization (IVF) does. Male and female frogs are injected with an ovulation hormone, placed together in a mating bucket, and left for 4-6 h or overnight to mate. Embryos are then collected, treated with cysteine to remove the sticky jelly coat, and used for injections and/or raised to later stages or adulthood. For embryos raised past free-swimming stages, the cysteine step can optionally be skipped, and tadpoles can be allowed to hatch naturally from the jelly coat. Although there are many similarities between natural mating protocols for Xenopus laevis and Xenopus tropicalis, there are key differences such as hormone dosage, timing of ovulation, and embryo incubation temperature. Here we provide a specific protocol for inducing natural matings in X. tropicalis.

DLG5 variants are associated with multiple congenital anomalies including ciliopathy phenotypes.

BACKGROUND: Cilia are dynamic cellular extensions that generate and sense signals to orchestrate proper development and tissue homeostasis. They rely on the underlying polarisation of cells to participate in signalling. Cilia dysfunction is a well-known cause of several diseases that affect multiple organ systems including the kidneys, brain, heart, respiratory tract, skeleton and retina. METHODS: Among individuals from four unrelated families, we identified variants in discs large 5 (DLG5) that manifested in a variety of pathologies. In our proband, we also examined patient tissues. We depleted dlg5 in Xenopus tropicalis frog embryos to generate a loss-of-function model. Finally, we tested the pathogenicity of DLG5 patient variants through rescue experiments in the frog model. RESULTS: Patients with variants of DLG5 were found to have a variety of phenotypes including cystic kidneys, nephrotic syndrome, hydrocephalus, limb abnormalities, congenital heart disease and craniofacial malformations. We also observed a loss of cilia in cystic kidney tissue of our proband. Knockdown of dlg5 in Xenopus embryos recapitulated many of these phenotypes and resulted in a loss of cilia in multiple tissues. Unlike introduction of wildtype DLG5 in frog embryos depleted of dlg5, introduction of DLG5 patient variants was largely ineffective in restoring proper ciliation and tissue morphology in the kidney and brain suggesting that the variants were indeed detrimental to function. CONCLUSION: These findings in both patient tissues and Xenopus shed light on how mutations in DLG5 may lead to tissue-specific manifestations of disease. DLG5 is essential for cilia and many of the patient phenotypes are in the ciliopathy spectrum.

Novel compound heterozygous variants in NHLRC2 in a patient with FINCA syndrome.

Two variants in the ubiquitously expressed NHLRC2 gene have been reported to cause a lethal fibrotic cerebropulmonary disease termed fibrosis, neurodegeneration, and cerebral angiomatosis (FINCA) syndrome in three Finnish children. Our objective was to determine the genetic basis of disease in a new patient with clinical features of FINCA syndrome using whole-exome sequencing (WES) and confirmation by Sanger sequencing. The patient has one known and one novel variant in NHLRC2 (c.442T>G, p.D148Y and c.428C>A, p.H143P, respectively). p.H143P is extremely rare and is not present in the gnomAD database of >140,000 allele sequences from healthy humans. Both variants affect the highly conserved N-terminal thioredoxin (Trx)-like domain of NHLRC2 and are predicted to be damaging. We conclude that a compound heterozygous combination of a known and a novel variant in NHLRC2 causes FINCA syndrome in a 2-year-old Ukrainian patient, underscoring the importance of NHLRC2 as a central regulator of fibrosis.

De novo damaging variants associated with congenital heart diseases contribute to the connectome.

Congenital heart disease (CHD) survivors are at risk for neurodevelopmental disability (NDD), and recent studies identify genes associated with both disorders, suggesting that NDD in CHD survivors may be of genetic origin. Genes contributing to neurogenesis, dendritic development and synaptogenesis organize neural elements into networks known as the connectome. We hypothesized that NDD in CHD may be attributable to genes altering both neural connectivity and cardiac patterning. To assess the contribution of de novo variants (DNVs) in connectome genes, we annotated 229 published NDD genes for connectome status and analyzed data from 3,684 CHD subjects and 1,789 controls for connectome gene mutations. CHD cases had more protein truncating and deleterious missense DNVs among connectome genes compared to controls (OR = 5.08, 95%CI:2.81-9.20, Fisher's exact test P = 6.30E-11). When removing three known syndromic CHD genes, the findings remained significant (OR = 3.69, 95%CI:2.02-6.73, Fisher's exact test P = 1.06E-06). In CHD subjects, the top 12 NDD genes with damaging DNVs that met statistical significance after Bonferroni correction (PTPN11, CHD7, CHD4, KMT2A, NOTCH1, ADNP, SMAD2, KDM5B, NSD2, FOXP1, MED13L, DYRK1A; one-tailed binomial test P ≤ 4.08E-05) contributed to the connectome. These data suggest that NDD in CHD patients may be attributable to genes that alter both cardiac patterning and the connectome.

Novel truncating mutations in CTNND1 cause a dominant craniofacial and cardiac syndrome.

CTNND1 encodes the p120-catenin (p120) protein, which has a wide range of functions, including the maintenance of cell-cell junctions, regulation of the epithelial-mesenchymal transition and transcriptional signalling. Due to advances in next-generation sequencing, CTNND1 has been implicated in human diseases including cleft palate and blepharocheilodontic (BCD) syndrome albeit only recently. In this study, we identify eight novel protein-truncating variants, six de novo, in 13 participants from nine families presenting with craniofacial dysmorphisms including cleft palate and hypodontia, as well as congenital cardiac anomalies, limb dysmorphologies and neurodevelopmental disorders. Using conditional deletions in mice as well as CRISPR/Cas9 approaches to target CTNND1 in Xenopus, we identified a subset of phenotypes that can be linked to p120-catenin in epithelial integrity and turnover, and additional phenotypes that suggest mesenchymal roles of CTNND1. We propose that CTNND1 variants have a wider developmental role than previously described and that variations in this gene underlie not only cleft palate and BCD but may be expanded to a broader velocardiofacial-like syndrome.

Histone H2B monoubiquitination regulates heart development via epigenetic control of cilia motility.

Genomic analyses of patients with congenital heart disease (CHD) have identified significant contribution from mutations affecting cilia genes and chromatin remodeling genes; however, the mechanism(s) connecting chromatin remodeling to CHD is unknown. Histone H2B monoubiquitination (H2Bub1) is catalyzed by the RNF20 complex consisting of RNF20, RNF40, and UBE2B. Here, we show significant enrichment of loss-of-function mutations affecting H2Bub1 in CHD patients (enrichment 6.01, P = 1.67 × 10-03), some of whom had abnormal laterality associated with ciliary dysfunction. In Xenopus, knockdown of rnf20 and rnf40 results in abnormal heart looping, defective development of left-right (LR) asymmetry, and impaired cilia motility. Rnf20, Rnf40, and Ube2b affect LR patterning and cilia synergistically. Examination of global H2Bub1 level in Xenopus embryos shows that H2Bub1 is developmentally regulated and requires Rnf20. To examine gene-specific H2Bub1, we performed ChIP-seq of mouse ciliated and nonciliated tissues and showed tissue-specific H2Bub1 marks significantly enriched at cilia genes including the transcription factor Rfx3 Rnf20 knockdown results in decreased levels of rfx3 mRNA in Xenopus, and exogenous rfx3 can rescue the Rnf20 depletion phenotype. These data suggest that Rnf20 functions at the Rfx3 locus regulating cilia motility and cardiac situs and identify H2Bub1 as an upstream transcriptional regulator controlling tissue-specific expression of cilia genes. Our findings mechanistically link the two functional gene ontologies that have been implicated in human CHD: chromatin remodeling and cilia function.

Visualizing flow in an intact CSF network using optical coherence tomography: implications for human congenital hydrocephalus.

Cerebrospinal fluid (CSF) flow in the brain ventricles is critical for brain development. Altered CSF flow dynamics have been implicated in congenital hydrocephalus (CH) characterized by the potentially lethal expansion of cerebral ventricles if not treated. CH is the most common neurosurgical indication in children effecting 1 per 1000 infants. Current treatment modalities are limited to antiquated brain surgery techniques, mostly because of our poor understanding of the CH pathophysiology. We lack model systems where the interplay between ependymal cilia, embryonic CSF flow dynamics and brain development can be analyzed in depth. This is in part due to the poor accessibility of the vertebrate ventricular system to in vivo investigation. Here, we show that the genetically tractable frog Xenopus tropicalis, paired with optical coherence tomography imaging, provides new insights into CSF flow dynamics and role of ciliary dysfunction in hydrocephalus pathogenesis. We can visualize CSF flow within the multi-chambered ventricular system and detect multiple distinct polarized CSF flow fields. Using CRISPR/Cas9 gene editing, we modeled human L1CAM and CRB2 mediated aqueductal stenosis. We propose that our high-throughput platform can prove invaluable for testing candidate human CH genes to understand CH pathophysiology.

Patients with common variable immunodeficiency with autoimmune cytopenias exhibit hyperplastic yet inefficient germinal center responses.

BACKGROUND: The lack of pathogen-protective, isotype-switched antibodies in patients with common variable immunodeficiency (CVID) suggests germinal center (GC) hypoplasia, yet a subset of patients with CVID is paradoxically affected by autoantibody-mediated autoimmune cytopenias (AICs) and lymphadenopathy. OBJECTIVE: We sought to compare the physical characteristics and immunologic output of GC responses in patients with CVID with AIC (CVID+AIC) and without AIC (CVID-AIC). METHODS: We analyzed GC size and shape in excisional lymph node biopsy specimens from 14 patients with CVID+AIC and 4 patients with CVID-AIC. Using paired peripheral blood samples, we determined how AICs specifically affected B-and T-cell compartments and antibody responses in patients with CVID. RESULTS: We found that patients with CVID+AIC displayed irregularly shaped hyperplastic GCs, whereas GCs were scarce and small in patients with CVID-AIC. GC hyperplasia was also evidenced by an increase in numbers of circulating follicular helper T cells, which correlated with decreased regulatory T-cell frequencies and function. In addition, patients with CVID+AIC had serum endotoxemia associated with a dearth of isotype-switched memory B cells that displayed significantly lower somatic hypermutation frequencies than their counterparts with CVID-AIC. Moreover, IgG+ B cells from patients with CVID+AIC expressed VH4-34-encoded antibodies with unmutated Ala-Val-Tyr and Asn-His-Ser motifs, which recognize both erythrocyte I/i self-antigens and commensal bacteria. CONCLUSIONS: Patients with CVID+AIC do not contain mucosal microbiota and exhibit hyperplastic yet inefficient GC responses that favor the production of untolerized IgG+ B-cell clones that recognize both commensal bacteria and hematopoietic I/i self-antigens.

Cytotoxic T-Lymphocyte-Associated Protein 4 Haploinsufficiency-Associated Inflammation Can Occur Independently of T-Cell Hyperproliferation.

Located contiguously on the long arm of the second chromosome are gene paralogs encoding the immunoglobulin-family co-activation receptors CD28 and cytotoxic T-lymphocyte-associated protein 4 (CTLA4). CD28 and CTLA4 share the same B7 ligands yet each provides opposing proliferative signals to T cells. Herein, we describe for the first time two unrelated subjects with coexisting CD28 and CTLA4 haploinsufficiency due to heterozygous microdeletions of chromosome 2q. Although their clinical phenotype, multi-organ inflammatory disease, is superficially similar to that of CTLA4 haploinsufficient autoimmune lymphoproliferative syndrome type V (ALPS5) patients, we demonstrate our subjects' underlying immunopathology to be distinct. Unlike ALPS5 T cells which hyperproliferate to T-cell receptor-mediated activation and infiltrate organs, T cells from our subjects are hypoproliferative and do not. Instead of T cell infiltrates, biopsies of affected subject tissues demonstrated infiltrates of lineage negative lymphoid cells. This histologic feature correlated with significant increases in circulating type 3 innate lymphoid cells (ILC3s) and ILC3 cytokines, interleukin 22, and interleukin-17A. CTLA4-Ig monotherapy, which we trialed in one subject, was remarkably effective in controlling inflammatory diseases, normalizing ILC3 frequencies, and reducing ILC3 cytokine concentrations.

The Psychiatric Cell Map Initiative: A Convergent Systems Biological Approach to Illuminating Key Molecular Pathways in Neuropsychiatric Disorders.

Although gene discovery in neuropsychiatric disorders, including autism spectrum disorder, intellectual disability, epilepsy, schizophrenia, and Tourette disorder, has accelerated, resulting in a large number of molecular clues, it has proven difficult to generate specific hypotheses without the corresponding datasets at the protein complex and functional pathway level. Here, we describe one path forward-an initiative aimed at mapping the physical and genetic interaction networks of these conditions and then using these maps to connect the genomic data to neurobiology and, ultimately, the clinic. These efforts will include a team of geneticists, structural biologists, neurobiologists, systems biologists, and clinicians, leveraging a wide array of experimental approaches and creating a collaborative infrastructure necessary for long-term investigation. This initiative will ultimately intersect with parallel studies that focus on other diseases, as there is a significant overlap with genes implicated in cancer, infectious disease, and congenital heart defects.

RAPGEF5 Regulates Nuclear Translocation of ß-Catenin.

Canonical Wnt signaling coordinates many critical aspects of embryonic development, while dysregulated Wnt signaling contributes to common diseases, including congenital malformations and cancer. The nuclear localization of ß-catenin is the defining step in pathway activation. However, despite intensive investigation, the mechanisms regulating ß-catenin nuclear transport remain undefined. In a patient with congenital heart disease and heterotaxy, a disorder of left-right patterning, we previously identified the guanine nucleotide exchange factor, RAPGEF5. Here, we demonstrate that RAPGEF5 regulates left-right patterning via Wnt signaling. In particular, RAPGEF5 regulates the nuclear translocation of ß-catenin independently of both ß-catenin cytoplasmic stabilization and the importin ß1/Ran-mediated transport system. We propose a model whereby RAPGEF5 activates the nuclear GTPases, Rap1a/b, to facilitate the nuclear transport of ß-catenin, defining a parallel nuclear transport pathway to Ran. Our results suggest new targets for modulating Wnt signaling in disease states.

A p190BRhoGAP mutation and prolonged RhoB activation in fatal systemic capillary leak syndrome.

We describe a fatal case of pediatric systemic capillary leak (Clarkson's disease) associated with a point mutation in p190BRhoGAP. Dermal microvascular endothelial cells (ECs) isolated from this patient form monolayers with similar levels and distribution of junctional proteins and transendothelial electrical resistance compared with normal human dermal microvascular ECs. However, patient-derived ECs demonstrate a greater increase in permeability and impaired recovery of barrier function in response to tumor necrosis factor (TNF) compared with normal donor EC cultures. TNF transiently activates RhoB in ECs coincident with developing leak, and inactivation of RhoB correlates with barrier recovery. The mutation in p190BRhoGAP impairs RhoB inactivation, and the mutant phenotype of patient-derived ECs is replicated by siRNA knockdown of p190BRhoGAP in normal ECs. These data suggest a previously unknown function for p190BRhoGAP in control of capillary EC barrier function that may also be important in acquired systemic capillary leak associated with critical illness in humans.

Xenopus as a model organism for birth defects-Congenital heart disease and heterotaxy.

Congenital heart disease is the leading cause of birth defects, affecting 9 out of 1000 newborns each year. A particularly severe form of congenital heart disease is heterotaxy, a disorder of left-right development. Despite aggressive surgical management, patients with heterotaxy have poor survival rates and severe morbidity due to their complex congenital heart disease. Recent genetic analysis of affected patients has found novel candidate genes for heterotaxy although their underlying mechanisms remain unknown. In this review, we discuss the importance and challenges of birth defects research including high locus heterogeneity and few second alleles that make defining disease causality difficult. A powerful strategy moving forward is to analyze these candidate genes in a high-throughput human disease model. Xenopus is ideal for these studies. We present multiple examples demonstrating the power of Xenopus in discovering new biology from the analysis of candidate heterotaxy genes such as GALNT11, NEK2 and BCOR. These genes have diverse roles in embryos and have led to a greater understanding of complex signaling pathways and basic developmental biology. It is our hope that the mechanistic analysis of these candidate genes in Xenopus enabled by next generation sequencing of patients will provide clinicians with a greater understanding of patient pathophysiology allowing more precise and personalized medicine, to help patients more effectively in the future.

The heterotaxy gene GALNT11 glycosylates Notch to orchestrate cilia type and laterality.

Heterotaxy is a disorder of left-right body patterning, or laterality, that is associated with major congenital heart disease. The aetiology and mechanisms underlying most cases of human heterotaxy are poorly understood. In vertebrates, laterality is initiated at the embryonic left-right organizer, where motile cilia generate leftward flow that is detected by immotile sensory cilia, which transduce flow into downstream asymmetric signals. The mechanism that specifies these two cilia types remains unknown. Here we show that the N-acetylgalactosamine-type O-glycosylation enzyme GALNT11 is crucial to such determination. We previously identified GALNT11 as a candidate disease gene in a patient with heterotaxy, and now demonstrate, in Xenopus tropicalis, that galnt11 activates Notch signalling. GALNT11 O-glycosylates human NOTCH1 peptides in vitro, thereby supporting a mechanism of Notch activation either by increasing ADAM17-mediated ectodomain shedding of the Notch receptor or by modification of specific EGF repeats. We further developed a quantitative live imaging technique for Xenopus left-right organizer cilia and show that Galnt11-mediated Notch1 signalling modulates the spatial distribution and ratio of motile and immotile cilia at the left-right organizer. galnt11 or notch1 depletion increases the ratio of motile cilia at the expense of immotile cilia and produces a laterality defect reminiscent of loss of the ciliary sensor Pkd2. By contrast, Notch overexpression decreases this ratio, mimicking the ciliopathy primary ciliary dyskinesia. Together our data demonstrate that Galnt11 modifies Notch, establishing an essential balance between motile and immotile cilia at the left-right organizer to determine laterality, and reveal a novel mechanism for human heterotaxy.

Accelerated gene evolution and subfunctionalization in the pseudotetraploid frog Xenopus laevis.

BACKGROUND: Ancient whole genome duplications have been implicated in the vertebrate and teleost radiations, and in the emergence of diverse angiosperm lineages, but the evolutionary response to such a perturbation is still poorly understood. The African clawed frog Xenopus laevis experienced a relatively recent tetraploidization ~40 million years ago. Analysis of the considerable amount of EST sequence available for this species together with the genome sequence of the related diploid Xenopus tropicalis provides a unique opportunity to study the genomic response to whole genome duplication. RESULTS: We identified 2218 gene triplets in which a single gene in X. tropicalis corresponds to precisely two co-orthologous genes in X. laevis--the largest such collection published from any duplication event in animals. Analysis of these triplets reveals accelerated evolution or relaxation of constraint in the peptides of the X. laevis pairs compared with the orthologous sequences in X. tropicalis and other vertebrates. In contrast, single-copy X. laevis genes do not show this acceleration. Duplicated genes can differ substantially in expression levels and patterns. We find no significant difference in gene content in the duplicated set, versus the single-copy set based on molecular and biological function ontologies. CONCLUSION: These results support a scenario in which duplicate genes are retained through a process of subfunctionalization and/or relaxation of constraint on both copies of an ancestral gene.