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.

The dishevelled associated activator of morphogenesis protein 2 (Daam2) regulates neural tube closure.

BACKGROUND: The Wnt signaling pathway is highly conserved in metazoans and regulates a large array of cellular processes including motility, polarity and fate determination, and stem cell homeostasis. Modulation of the actin cytoskeleton via the non-canonical Wnt pathway regulate cell polarity and cell migration that are required for proper vertebrate gastrulation and subsequent neurulation. However, the mechanism(s) of how the non-canonical pathway mediates actin cytoskeleton modulation is not fully understood. RESULTS: Herein, we characterize the role of the Formin-homology protein; dishevelled associated activator of morphogenesis 2 (Daam2) protein in the Wnt signaling pathway. Co-immunoprecipitation assays confirm the binding of Daam2 to dishevelled2 (Dvl2) as well as the domains within these proteins required for interaction; additionally, the interaction between Daam2 and Dvl2 was Wnt-regulated. Sub-cellular localization studies reveal Daam2 is cytoplasmic and regulates the cellular actin cytoskeleton by modulating actin filament formation. During Xenopus development, a knockdown or loss of Daam2 specifically produces neural tube closure defects indicative of a role in non-canonical signaling. Additionally, our studies did not identify any role for Daam2 in canonical Wnt signaling in mammalian culture cells or the Xenopus embryo. CONCLUSIONS: Our studies together identify Daam2 as a component of the non-canonical Wnt pathway and Daam2 is a regulator of neural tube morphogenesis during vertebrate development.

Mink1 regulates spemann organizer cell fate in the xenopus gastrula via Hmga2.

Congenital Heart Disease (CHD) is the most common birth defect and leading cause of infant mortality, yet molecular mechanisms explaining CHD remain mostly unknown. Sequencing studies are identifying CHD candidate genes at a brisk rate including MINK1, a serine/threonine kinase. However, a plausible molecular mechanism connecting CHD and MINK1 is unknown. Here, we reveal that mink1 is required for proper heart development due to its role in left-right patterning. Mink1 regulates canonical Wnt signaling to define the cell fates of the Spemann Organizer and the Left-Right Organizer, a ciliated structure that breaks bilateral symmetry in the vertebrate embryo. To identify Mink1 targets, we applied an unbiased proteomics approach and identified the high mobility group architectural transcription factor, Hmga2. We report that Hmga2 is necessary and sufficient for regulating Spemann's Organizer. Indeed, we demonstrate that Hmga2 can induce Spemann Organizer cell fates even when ß-catenin, a critical effector of the Wnt signaling pathway, is depleted. In summary, we discover a transcription factor, Hmga2, downstream of Mink1 that is critical for the regulation of Spemann's Organizer, as well as the LRO, defining a plausible mechanism for CHD.

AP2A2 mutation and defective endocytosis in a Malian family with hereditary spastic paraplegia.

Hereditary spastic paraplegia (HSP) comprises a large group of neurogenetic disorders characterized by progressive lower extremity spasticity. Neurological evaluation and genetic testing were completed in a Malian family with early-onset HSP. Three children with unaffected consanguineous parents presented with symptoms consistent with childhood-onset complicated HSP. Neurological evaluation found lower limb weakness, spasticity, dysarthria, seizures, and intellectual disability. Brain MRI showed corpus callosum thinning with cortical and spinal cord atrophy, and an EEG detected slow background in the index patient. Whole exome sequencing identified a homozygous missense variant in the adaptor protein (AP) complex 2 alpha-2 subunit (AP2A2) gene. Western blot analysis showed reduced levels of AP2A2 in patient-iPSC derived neuronal cells. Endocytosis of transferrin receptor (TfR) was decreased in patient-derived neurons. In addition, we observed increased axon initial segment length in patient-derived neurons. Xenopus tropicalis tadpoles with ap2a2 knockout showed cerebral edema and progressive seizures. Immunoprecipitation of the mutant human AP-2-appendage alpha-C construct showed defective binding to accessory proteins. We report AP2A2 as a novel genetic entity associated with HSP and provide functional data in patient-derived neuron cells and a frog model. These findings expand our understanding of the mechanism of HSP and improve the genetic diagnosis of this condition.

Biallelic CRELD1 variants cause a multisystem syndrome, including neurodevelopmental phenotypes, cardiac dysrhythmias, and frequent infections.

PURPOSE: We sought to delineate a multisystem disorder caused by recessive cysteine-rich with epidermal growth factor-like domains 1 (CRELD1) gene variants. METHODS: The impact of CRELD1 variants was characterized through an international collaboration utilizing next-generation DNA sequencing, gene knockdown, and protein overexpression in Xenopus tropicalis, and in vitro analysis of patient immune cells. RESULTS: Biallelic variants in CRELD1 were found in 18 participants from 14 families. Affected individuals displayed an array of phenotypes involving developmental delay, early-onset epilepsy, and hypotonia, with about half demonstrating cardiac arrhythmias and some experiencing recurrent infections. Most harbored a frameshift in trans with a missense allele, with 1 recurrent variant, p.(Cys192Tyr), identified in 10 families. X tropicalis tadpoles with creld1 knockdown displayed developmental defects along with increased susceptibility to induced seizures compared with controls. Additionally, human CRELD1 harboring missense variants from affected individuals had reduced protein function, indicated by a diminished ability to induce craniofacial defects when overexpressed in X tropicalis. Finally, baseline analyses of peripheral blood mononuclear cells showed similar proportions of immune cell subtypes in patients compared with healthy donors. CONCLUSION: This patient cohort, combined with experimental data, provide evidence of a multisystem clinical syndrome mediated by recessive variants in CRELD1.

Expanding EMC foldopathies: Topogenesis deficits alter the neural crest.

The endoplasmic reticulum (ER) membrane protein complex (EMC) is essential for the insertion of a wide variety of transmembrane proteins into the plasma membrane across cell types. Each EMC is composed of Emc1-7, Emc10, and either Emc8 or Emc9. Recent human genetics studies have implicated variants in EMC genes as the basis for a group of human congenital diseases. The patient phenotypes are varied but appear to affect a subset of tissues more prominently than others. Namely, craniofacial development seems to be commonly affected. We previously developed an array of assays in Xenopus tropicalis to assess the effects of emc1 depletion on the neural crest, craniofacial cartilage, and neuromuscular function. We sought to extend this approach to additional EMC components identified in patients with congenital malformations. Through this approach, we determine that EMC9 and EMC10 are important for neural crest development and the development of craniofacial structures. The phenotypes observed in patients and our Xenopus model phenotypes similar to EMC1 loss of function likely due to a similar mechanism of dysfunction in transmembrane protein topogenesis.

DAAM2 Variants Cause Nephrotic Syndrome via Actin Dysregulation.

The discovery of >60 monogenic causes of nephrotic syndrome (NS) has revealed a central role for the actin regulators RhoA/Rac1/Cdc42 and their effectors, including the formin INF2. By whole-exome sequencing (WES), we here discovered bi-allelic variants in the formin DAAM2 in four unrelated families with steroid-resistant NS. We show that DAAM2 localizes to the cytoplasm in podocytes and in kidney sections. Further, the variants impair DAAM2-dependent actin remodeling processes: wild-type DAAM2 cDNA, but not cDNA representing missense variants found in individuals with NS, rescued reduced podocyte migration rate (PMR) and restored reduced filopodia formation in shRNA-induced DAAM2-knockdown podocytes. Filopodia restoration was also induced by the formin-activating molecule IMM-01. DAAM2 also co-localizes and co-immunoprecipitates with INF2, which is intriguing since variants in both formins cause NS. Using in vitro bulk and TIRF microscopy assays, we find that DAAM2 variants alter actin assembly activities of the formin. In a Xenopus daam2-CRISPR knockout model, we demonstrate actin dysregulation in vivo and glomerular maldevelopment that is rescued by WT-DAAM2 mRNA. We conclude that DAAM2 variants are a likely cause of monogenic human SRNS due to actin dysregulation in podocytes. Further, we provide evidence that DAAM2-associated SRNS may be amenable to treatment using actin regulating compounds.

Paired Box 9 (PAX9), the RNA polymerase II transcription factor, regulates human ribosome biogenesis and craniofacial development.

Dysregulation of ribosome production can lead to a number of developmental disorders called ribosomopathies. Despite the ubiquitous requirement for these cellular machines used in protein synthesis, ribosomopathies manifest in a tissue-specific manner, with many affecting the development of the face. Here we reveal yet another connection between craniofacial development and making ribosomes through the protein Paired Box 9 (PAX9). PAX9 functions as an RNA Polymerase II transcription factor to regulate the expression of proteins required for craniofacial and tooth development in humans. We now expand this function of PAX9 by demonstrating that PAX9 acts outside of the cell nucleolus to regulate the levels of proteins critical for building the small subunit of the ribosome. This function of PAX9 is conserved to the organism Xenopus tropicalis, an established model for human ribosomopathies. Depletion of pax9 leads to craniofacial defects due to abnormalities in neural crest development, a result consistent with that found for depletion of other ribosome biogenesis factors. This work highlights an unexpected layer of how the making of ribosomes is regulated in human cells and during embryonic development.

A nonsense variant in Rap Guanine Nucleotide Exchange Factor 5 (RAPGEF5) is associated with equine familial isolated hypoparathyroidism in Thoroughbred foals.

Idiopathic hypocalcemia in Thoroughbred (TB) foals causes tetany and seizures and is invariably fatal. Based upon the similarity of this disease with human familial hypoparathyroidism and occurrence only in the TB breed, we conducted a genetic investigation on two affected TB foals. Familial hypoparathyroidism was identified, and pedigree analysis suggested an autosomal recessive (AR) mode of inheritance. We performed whole-genome sequencing of the two foals, their unaffected dams and four unaffected, unrelated TB horses. Both homozygosity mapping and an association analysis were used to prioritize potential genetic variants. Of the 2,808 variants that significantly associated with the phenotype using an AR mode of inheritance (P<0.02) and located within a region of homozygosity, 1,507 (54%) were located in a 9.7 Mb region on chr4 (44.9-54.6 Mb). Within this region, a nonsense variant (RAPGEF5 c.2624C>A,p.Ser875*) was significantly associated with the hypoparathyroid phenotype (Pallelic = 0.008). Affected foals were homozygous for the variant, with two additional affected foals subsequently confirmed in 2019. Necropsies of all affected foals failed to identify any histologically normal parathyroid glands. Because the nonsense mutation in RAPGEF5 was near the C-terminal end of the protein, the impact on protein function was unclear. Therefore, we tested the variant in our Xenopus overexpression model and demonstrated RAPGEF5 loss-of-function. This RAPGEF5 variant represents the first genetic variant for hypoparathyroidism identified in any domestic animal species.

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.

Nucleoporin NUP205 plays a critical role in cilia and congenital disease.

Ciliopathies affect a variety of tissues during development including the heart, kidneys, respiratory tract, and retina. Though an increasing number of monogenic causes of ciliopathies have been described, many remain unexplained. Recently, recessive variants in NUP93 and NUP205 encoding two proteins of the inner ring of the nuclear pore complex were implicated as causes of steroid resistant nephrotic syndrome. In addition, we previously found that the inner ring nucleoporins NUP93 and NUP188 function in proper left-right patterning in developing embryos via a role at the cilium. Here, we describe the role of an additional inner ring nucleoporin NUP205 in cilia biology and establishment of normal organ situs. Using knockdown in Xenopus, we show that Nup205 depletion results in loss of cilia and abnormal cardiac morphology. Furthermore, by transmission electron microscopy, we observe a loss of cilia and mispositioning of intracellular ciliary structures such as basal bodies and rootlets upon depleting inner ring nucleoporins. We describe a model wherein NUP93 interacting with either NUP188 or NUP205 is necessary for cilia. We thus provide evidence that dysregulation of inner ring nucleoporin genes that have been identified in patients may contribute to pathogenesis through cilia dysfunction.

A retrospective cohort analysis of the Yale pediatric genomics discovery program.

The Pediatric Genomics Discovery Program (PGDP) at Yale uses next-generation sequencing (NGS) and translational research to evaluate complex patients with a wide range of phenotypes suspected to have rare genetic diseases. We conducted a retrospective cohort analysis of 356 PGDP probands evaluated between June 2015 and July 2020, querying our database for participant demographics, clinical characteristics, NGS results, and diagnostic and research findings. The three most common phenotypes among the entire studied cohort (n = 356) were immune system abnormalities (n = 105, 29%), syndromic or multisystem disease (n = 103, 29%), and cardiovascular system abnormalities (n = 62, 17%). Of 216 patients with final classifications, 77 (36%) received new diagnoses and 139 (64%) were undiagnosed; the remaining 140 patients were still actively being investigated. Monogenetic diagnoses were found in 67 (89%); the largest group had variants in known disease genes but with new contributions such as novel variants (n = 31, 40%) or expanded phenotypes (n = 14, 18%). Finally, five PGDP diagnoses (8%) were suggestive of novel gene-to-phenotype relationships. A broad range of patients can benefit from single subject studies combining NGS and functional molecular analyses. All pediatric providers should consider further genetics evaluations for patients lacking precise molecular diagnoses.

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.

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.

WDR5 regulates left-right patterning via chromatin-dependent and -independent functions.

Congenital heart disease (CHD) is a major cause of infant mortality and morbidity, yet the genetic causes and mechanisms remain opaque. In a patient with CHD and heterotaxy, a disorder of left-right (LR) patterning, a de novo mutation was identified in the chromatin modifier gene WDR5 WDR5 acts as a scaffolding protein in the H3K4 methyltransferase complex, but a role in LR patterning is unknown. Here, we show that Wdr5 depletion leads to LR patterning defects in Xenopus via its role in ciliogenesis. Unexpectedly, we find a dual role for WDR5 in LR patterning. First, WDR5 is expressed in the nuclei of monociliated cells of the LR organizer (LRO) and regulates foxj1 expression. LR defects in wdr5 morphants can be partially rescued with the addition of foxj1 Second, WDR5 localizes to the bases of cilia. Using a mutant form of WDR5, we demonstrate that WDR5 also has an H3K4-independent role in LR patterning. Guided by the patient phenotype, we identify multiple roles for WDR5 in LR patterning, providing plausible mechanisms for its role in ciliopathies like heterotaxy and CHD.

RPSA, a candidate gene for isolated congenital asplenia, is required for pre-rRNA processing and spleen formation in Xenopus.

A growing number of tissue-specific inherited disorders are associated with impaired ribosome production, despite the universal requirement for ribosome function. Recently, mutations in RPSA, a protein component of the small ribosomal subunit, were discovered to underlie approximately half of all isolated congenital asplenia cases. However, the mechanisms by which mutations in this ribosome biogenesis factor lead specifically to spleen agenesis remain unknown, in part due to the lack of a suitable animal model for study. Here we reveal that RPSA is required for normal spleen development in the frog, Xenopus tropicalis Depletion of Rpsa in early embryonic development disrupts pre-rRNA processing and ribosome biogenesis, and impairs expression of the key spleen patterning genes nkx2-5, bapx1 and pod1 in the spleen anlage. Importantly, we also show that whereas injection of human RPSA mRNA can rescue both pre-rRNA processing and spleen patterning, injection of human mRNA bearing a common disease-associated mutation cannot. Together, we present the first animal model of RPSA-mediated asplenia and reveal a crucial requirement for RPSA in pre-rRNA processing and molecular patterning during early Xenopus development.

Functional testing for variant prioritization in a family with long QT syndrome.

Next-generation sequencing platforms are being increasingly applied in clinical genetic settings for evaluation of families with suspected heritable disease. These platforms potentially improve the diagnostic yield beyond that of disease-specific targeted gene panels, but also increase the number of rare or novel genetic variants that may confound precise diagnostics. Here, we describe a functional testing approach used to interpret the results of whole exome sequencing (WES) in a family presenting with syncope and sudden death. One individual had a prolonged QT interval on electrocardiogram (ECG) and carried a diagnosis of long QT syndrome (LQTS), but a second individual did not meet criteria for LQTS. Filtering WES results for uncommon variants with arrhythmia association identified four for further analyses. In silico analyses indicated that two of these variants, KCNH2 p.(Cys555Arg) and KCNQ1 p.(Arg293Cys), were likely to be causal in this family's LQTS. We subsequently performed functional characterization of these variants in a heterologous expression system. The expression of KCNQ1-Arg293Cys did not show a deleterious phenotype but KCNH2-Cys555Arg demonstrated a loss-of-function phenotype that was partially dominant. Our stepwise approach identified a precise genetic etiology in this family, which resulted in the establishment of a LQTS diagnosis in the second individual as well as an additional asymptomatic family member, enabling personalized clinical management. Given its ability to aid in the diagnosis, the application of functional characterization should be considered as a value adjunct to in silico analyses of WES.

Expansion of NEUROD2 phenotypes to include developmental delay without seizures.

De novo heterozygous variants in the brain-specific transcription factor Neuronal Differentiation Factor 2 (NEUROD2) have been recently associated with early-onset epileptic encephalopathy and developmental delay. Here, we report an adolescent with developmental delay without seizures who was found to have a novel de novo heterozygous NEUROD2 missense variant, p.(Leu163Pro). Functional testing using an in vivo assay of neuronal differentiation in Xenopus laevis tadpoles demonstrated that the patient variant of NEUROD2 displays minimal protein activity, strongly suggesting a loss of function effect. In contrast, a second rare NEUROD2 variant, p.(Ala235Thr), identified in an adolescent with developmental delay but lacking parental studies for inheritance, showed normal in vivo NEUROD2 activity. We thus provide clinical, genetic, and functional evidence that NEUROD2 variants can lead to developmental delay without accompanying early-onset seizures, and demonstrate how functional testing can complement genetic data when determining variant pathogenicity.

Alkylglycerol monooxygenase, a heterotaxy candidate gene, regulates left-right patterning via Wnt signaling.

Congenital heart disease (CHD) is a major cause of morbidity in the pediatric population yet its genetic and molecular causes remain poorly defined. Previously, we identified AGMO as a candidate heterotaxy disease gene, a disorder of left-right (LR) patterning that can have a profound effect on cardiac function. AGMO is the only known alkylglycerol monooxygenase, an orphan tetrahydrobiopterin dependent enzyme that cleaves the ether linkage in alkylglycerols. However, whether AGMO plays a role in LR patterning was unexplored. Here we reveal that Agmo is required for correct development of the embryonic LR axis in Xenopus embryos recapitulating the patient's heterotaxy phenotype. Mechanistically, we demonstrate that Agmo is a regulator of canonical Wnt signaling, required during gastrulation for normal formation of the left - right organizer. Mutational analysis demonstrates that this function is dependent on Agmo's alkylglycerol monooxygenase activity. Together, our findings identify Agmo as a regulator of canonical Wnt signaling, demonstrate a role for Agmo in embryonic axis formation, and provide insight into the poorly understood developmental requirements for ether lipid cleavage.

A chromosome-scale genome assembly and dense genetic map for Xenopus tropicalis.

The Western clawed frog Xenopus tropicalis is a diploid model system for both frog genetics and developmental biology, complementary to the paleotetraploid X. laevis. Here we report a chromosome-scale assembly of the X. tropicalis genome, improving the previously published draft genome assembly through the use of new assembly algorithms, additional sequence data, and the addition of a dense genetic map. The improved genome enables the mapping of specific traits (e.g., the sex locus or Mendelian mutants) and the characterization of chromosome-scale synteny with other tetrapods. We also report an improved annotation of the genome that integrates deep transcriptome sequence from diverse tissues and stages. The exon-intron structures of these genes are highly conserved relative to both X. laevis and human, as are chromosomal linkages ("synteny") and local gene order. A network analysis of developmental gene expression will aid future studies.