Comparative analysis of the genomic organization of Pax9 and its conserved physical association with Nkx2-9 in the human, mouse, and pufferfish genomes

Hemichordate genomes and deuterostome origins

Acorn worms, also known as enteropneust (literally, 'gut-breathing') hemichordates, are marine invertebrates that share features with echinoderms and chordates. Together, these three phyla comprise the deuterostomes. Here we report the draft genome sequences of two acorn worms, Saccoglossus kowalevskii and Ptychodera flava. By comparing them with diverse bilaterian genomes, we identify shared traits that were probably inherited from the last common deuterostome ancestor, and then explore evolutionary trajectories leading from this ancestor to hemichordates, echinoderms and chordates. The hemichordate genomes exhibit extensive conserved synteny with amphioxus and other bilaterians, and deeply conserved non-coding sequences that are candidates for conserved gene-regulatory elements. Notably, hemichordates possess a deuterostome-specific genomic cluster of four ordered transcription factor genes, the expression of which is associated with the development of pharyngeal 'gill' slits, the foremost morphological innovation of early deuterostomes, and is probably central to their filter-feeding lifestyle. Comparative analysis reveals numerous deuterostome-specific gene novelties, including genes found in deuterostomes and marine microbes, but not other animals. The putative functions of these genes can be linked to physiological, metabolic and developmental specializations of the filter-feeding ancestor.

Genome-wide association mapping in dogs enables identification of the homeobox gene, NKX2-8, as a genetic component of neural tube defects in humans

Neural tube defects (NTDs) is a general term for central nervous system malformations secondary to a failure of closure or development of the neural tube. The resulting pathologies may involve the brain, spinal cord and/or vertebral column, in addition to associated structures such as soft tissue or skin. The condition is reported among the more common birth defects in humans, leading to significant infant morbidity and mortality. The etiology remains poorly understood but genetic, nutritional, environmental factors, or a combination of these, are known to play a role in the development of NTDs. The variable conditions associated with NTDs occur naturally in dogs, and have been previously reported in the Weimaraner breed. Taking advantage of the strong linkage-disequilibrium within dog breeds we performed genome-wide association analysis and mapped a genomic region for spinal dysraphism, a presumed NTD, using 4 affected and 96 unaffected Weimaraners. The associated region on canine chromosome 8 (pgenome  =3.0 × 10(-5)), after 100,000 permutations, encodes 18 genes, including NKX2-8, a homeobox gene which is expressed in the developing neural tube. Sequencing NKX2-8 in affected Weimaraners revealed a G to AA frameshift mutation within exon 2 of the gene, resulting in a premature stop codon that is predicted to produce a truncated protein. The exons of NKX2-8 were sequenced in human patients with spina bifida and rare variants (rs61755040 and rs10135525) were found to be significantly over-represented (p=0.036). This is the first documentation of a potential role for NKX2-8 in the etiology of NTDs, made possible by investigating the molecular basis of naturally occurring mutations in dogs.

A novel g.-1258G>A mutation in a conserved putative regulatory element of PAX9 is associated with autosomal dominant molar hypodontia

Mutations in the transcription factor PAX9 which plays a critical role in the switching of odontogenic potential from the epithelium to the mesenchyme during tooth development cause autosomal dominant non-syndromic hypodontia primarily affecting molars. Linkage analysis on a family segregating autosomal dominant molar hypodontia with markers flanking and within PAX9 yielded a maximum multipoint LOD score of 3.6. No sequence variants were detected in the coding or 5'- and 3'-untranslated regions (UTRs) of PAX9. However, we identified a novel g.-1258G>A sequence variant in all affected individuals of the family but not in the unaffected family members or in 3088 control chromosomes. This mutation is within a putative 5'-regulatory sequence upstream of PAX9 highly conserved in primates, somewhat conserved in ungulates and carnivores but not conserved in rodents. Bioinformatics analysis of the sequence determined that there was no abolition or creation of a putative binding site for known transcription factors. Based on our previous findings that haploinsufficiency for PAX9 leads to hypodontia, we postulate that the g.-1258G>A variant reduces the expression of PAX9 which underlies the hypodontia phenotype in this family.

Organization of conserved elements near key developmental regulators in vertebrate genomes

Sequence conservation has traditionally been used as a means to target functional regions of complex genomes. In addition to its use in identifying coding regions of genes, the recent availability of whole genome data for a number of vertebrates has permitted high-resolution analyses of the noncoding "dark matter" of the genome. This has resulted in the identification of a large number of highly conserved sequence elements that appear to be preserved in all bony vertebrates. Further positional analysis of these conserved noncoding elements (CNEs) in the genome demonstrates that they cluster around genes involved in developmental regulation. This chapter describes the identification and characterization of these elements, with particular reference to their composition and organization.

Cloning of zebrafish nkx6.2 and a comprehensive analysis of the conserved transcriptional response to Hedgehog/Gli signaling in the zebrafish neural tube

Sonic Hedgehog (Shh) signaling helps pattern the vertebrate neural tube, in part by regulating the dorsal/ventral expression of a number of homeodomain containing transcription factors. These Hh responsive genes have been divided into two classes, with Class II genes being activated by Hh signaling and Class I genes being repressed by Hh signaling. While the transcriptional response to varying Hh levels is well defined in chick and mouse, it is only partially described in zebrafish, despite the fact that zebrafish has emerged as a powerful genetic system for the study of neural patterning. To better characterize the Hh response in the zebrafish neural tube, we cloned the zebrafish Class II Hh target genes nkx2.9 and nkx6.2. We then analyzed the expression of a number of Class I and Class II Hh responsive genes in wild type, Hh mutant, and Hh over-expressing zebrafish embryos. We show that expression of Class I and Class II genes is highly conserved in the vertebrate neural tube. Further, ventral-most Class II gene expression was completely lost in all Hh pathway mutants analyzed, indicating high levels of Hh signaling are blocked in all of these mutants. In contrast, more dorsally expressed genes were variably affected in different Hh pathway mutants, indicating mid-levels of Hh signaling are differentially affected. This comprehensive expression study provides an important tool for the characterization of Hh signaling in zebrafish and provides a sensitive assay for determining the degree to which newly identified zebrafish mutants affect Hh signaling.

Linking disease-associated genes to regulatory networks via promoter organization

Pathway- or disease-associated genes may participate in more than one transcriptional co-regulation network. Such gene groups can be readily obtained by literature analysis or by high-throughput techniques such as microarrays or protein-interaction mapping. We developed a strategy that defines regulatory networks by in silico promoter analysis, finding potentially co-regulated subgroups without a priori knowledge. Pairs of transcription factor binding sites conserved in orthologous genes (vertically) as well as in promoter sequences of co-regulated genes (horizontally) were used as seeds for the development of promoter models representing potential co-regulation. This approach was applied to a Maturity Onset Diabetes of the Young (MODY)-associated gene list, which yielded two models connecting functionally interacting genes within MODY-related insulin/glucose signaling pathways. Additional genes functionally connected to our initial gene list were identified by database searches with these promoter models. Thus, data-driven in silico promoter analysis allowed integrating molecular mechanisms with biological functions of the cell.

A genetic linkage map for Arctic char (Salvelinus alpinus): evidence for higher recombination rates and segregation distortion in hybrid versus pure strain mapping parents

We constructed a genetic linkage map for Arctic char (Salvelinus alpinus) using two backcrosses between genetically divergent strains. Forty-six linkage groups (expected = 39-41) and 19 homeologous affinities (expected = 25) were identified using 184 microsatellites, 129 amplified fragment length polymorphisms (AFLPs), 13 type I gene markers, and one phenotypic marker, SEX. Twenty-six markers remain unlinked. Female map distance (9.92 Morgans) was substantially higher than male map distance (3.90 Morgans) based on the most complete parental information (i.e., the F1 hybrids). Female recombination rates were often significantly higher than those of males across all pairwise comparisons within homologous chromosomal segments (average female to male ratios within families was 1.69:1). The female hybrid parent had significantly higher recombination rates than the pure strain female parent. Segregation distortion was detected in four linkage groups (4, 8, 13, 20) for both families. In family 3, only the largest fish were sampled for genotyping, suggesting that segregation distortion may represent regions possessing influences on growth. In family 2, almost all cases showing segregation distortion involved markers in the female hybrid parent.

Identification of Cis-regulatory Elements in the Mouse Pax9/Nkx2-9 Genomic Region: Implication for Evolutionary Conserved Synteny

We previously reported close physical linkage between Pax9 and Nkx2-9 in the human, mouse, and pufferfish (Fugu rubripes) genomes. In this study, we analyzed cis-regulatory elements of the two genes by comparative sequencing in the three species and by transgenesis in the mouse. We identified two regions including conserved noncoding sequences that possessed specific enhancer activities for expression of Pax9 in the medial nasal process and of Nkx2-9 in the ventral neural tube. Remarkably, the latter contained the consensus Gli-binding motif. Interestingly, the identified Pax9 cis-regulatory sequences were located in an intron of the neighboring gene Slc25a21. Close examination of an extended genomic interval around Pax9 revealed the presence of strong synteny conservation in the human, mouse, and Fugu genomes. We propose such an intersecting organization of cis-regulatory sequences in multigenic regions as a possible mechanism that maintains evolutionary conserved synteny.

Aberrant Pax1 and Pax9 expression in Jarcho-Levin syndrome: report of two Caucasian siblings and literature review

We report two consecutive Caucasian male siblings of nonconsanguineous parents autopsied at 22 and 13 weeks gestational age both with prenatal diagnosis of Jarcho-Levin syndrome (JLS). Segmentation anomalies of the vertebrae and ribs encompass a spectrum of syndromes with or without associated anomalies of other developmental fields, and include spondylothoracic dysostosis (STD), JLS, Casamassima-Morton-Nance (CMN) syndrome, and spondylocostal dysostosis (SCD), among others. In both these new JLS cases the autopsies confirmed that there were severe developmental alterations in the thoracic and vertebral skeleton (including "crab-like" thorax), accompanied in the older fetus by renal defects. Because vertebral development is controlled by a limited number of master genes including Pax1 and Pax9, we analyzed protein expression from these genes in these two cases compared to age-matched controls. Immunochemical analysis showed a significant reduction in levels of Pax1 and Pax9 protein expression in chondrocytes of the vertebral column. Implications for the etiology and pathogenesis of JLS and related disorders are discussed.

Regions of microsynteny in Magnaporthe grisea and Neurospora crassa

A bacterial artificial chromosome (BAC) clone containing 110,467 bp of genomic DNA from Magnaporthe grisea was sequenced, annotated, and compared to the genomes of Neurospora crassa, Candida albicans, and Saccharomyces cerevisiae. Twenty-six open reading frames (ORFs), involved in multiple biochemical pathways, were identified in the BAC sequence. A region of 53 kb, containing 18 of the 26 ORFs, was found to be syntenic to a portion of the N. crassa genome. Subregions of complete colinearity as well as interrupted colinearity were present. No synteny was evident with either C. albicans or S. cerevisiae. The identification of syntenic regions containing highly conserved genes across two genera that have been evolutionarily separated for approximately 200 million years elicits many biological questions as to the function and identity of these genes.