Comparative genomic analysis reveals the evolutionary conservation of Pax gene family

Hox gene expression during development of the phoronid Phoronopsis harmeri

Background

Phoronida is a small group of marine worm-like suspension feeders, which together with brachiopods and bryozoans form the clade Lophophorata. Although their development is well studied on the morphological level, data regarding gene expression during this process are scarce and restricted to the analysis of relatively few transcription factors. Here, we present a description of the expression patterns of Hox genes during the embryonic and larval development of the phoronid Phoronopsis harmeri.

Results

We identified sequences of eight Hox genes in the transcriptome of Ph. harmeri and determined their expression pattern during embryonic and larval development using whole mount in situ hybridization. We found that none of the Hox genes is expressed during embryonic development. Instead their expression is initiated in the later developmental stages, when the larval body is already formed. In the investigated initial larval stages the Hox genes are expressed in the non-collinear manner in the posterior body of the larvae: in the telotroch and the structures that represent rudiments of the adult worm. Additionally, we found that certain head-specific transcription factors are expressed in the oral hood, apical organ, preoral coelom, digestive system and developing larval tentacles, anterior to the Hox-expressing territories.

Conclusions

The lack of Hox gene expression during early development of Ph. harmeri indicates that the larval body develops without positional information from the Hox patterning system. Such phenomenon might be a consequence of the evolutionary intercalation of the larval form into an ancestral life cycle of phoronids. The observed Hox gene expression can also be a consequence of the actinotrocha representing a “head larva”, which is composed of the most anterior body region that is devoid of Hox gene expression. Such interpretation is further supported by the expression of head-specific transcription factors. This implies that the Hox patterning system is used for the positional information of the trunk rudiments and is, therefore, delayed to the later larval stages. We propose that a new body form was intercalated to the phoronid life cycle by precocious development of the anterior structures or by delayed development of the trunk rudiment in the ancestral phoronid larva.

The drug target genes show higher evolutionary conservation than non-target genes

Although evidence indicates that drug target genes share some common evolutionary features, there have been few studies analyzing evolutionary features of drug targets from an overall level. Therefore, we conducted an analysis which aimed to investigate the evolutionary characteristics of drug target genes. We compared the evolutionary conservation between human drug target genes and non-target genes by combining both the evolutionary features and network topological properties in human protein-protein interaction network. The evolution rate, conservation score and the percentage of orthologous genes of 21 species were included in our study. Meanwhile, four topological features including the average shortest path length, betweenness centrality, clustering coefficient and degree were considered for comparison analysis. Then we got four results as following: compared with non-drug target genes, 1) drug target genes had lower evolutionary rates; 2) drug target genes had higher conservation scores; 3) drug target genes had higher percentages of orthologous genes and 4) drug target genes had a tighter network structure including higher degrees, betweenness centrality, clustering coefficients and lower average shortest path lengths. These results demonstrate that drug target genes are more evolutionarily conserved than non-drug target genes. We hope that our study will provide valuable information for other researchers who are interested in evolutionary conservation of drug targets.

Phylogenetic analysis and expression patterns of Pax genes in the onychophoran Euperipatoides rowelli reveal a novel bilaterian Pax subfamily

Pax family genes encode a class of transcription factors that regulate various developmental processes. To shed light on the evolutionary history of these genes in Panarthropoda (Onychophora + Tardigrada + Arthropoda), we analyzed the Pax repertoire in the embryonic and adult transcriptomes of the onychophoran Euperipatoides rowelli. Our data revealed homologs of all five major bilaterian Pax subfamilies in this species, including Pax2/5/8, Pax4/6, Pox-neuro, Pax1/9/Pox-meso, and Pax3/7. In addition, we identified a new Pax member, pax-α, which does not fall into any other known Pax subfamily but instead clusters in the heterogenic Pax-α/β clade containing deuterostome, ecdysozoan, and lophotrochozoan gene sequences. These findings suggest that the last common bilaterian ancestor possessed six rather than five Pax genes, which have been retained in the panarthropod lineage. The expression data of Pax orthologs in the onychophoran embryo revealed distinctive patterns, some of which might be related to their ancestral roles in the last common panarthropod ancestor, whereas others might be specific to the onychophoran lineage. The derived roles include, for example, an involvement of pax2/5/8, pox-neuro, and pax3/7 in onychophoran nephridiogenesis, and an additional function of pax2/5/8 in the formation of the ventral and preventral organs. Furthermore, our transcriptomic analyses suggest that at least some Pax genes, including pax6 and pax-α, are expressed in the adult onychophoran head, although the corresponding functions remain to be clarified. The remarkable diversity of the Pax expression patterns highlights the functional and evolutionary plasticity of these genes in panarthropods.

The function of DrPax1b gene in the embryonic development of zebrafish

Vertebrate Pax1 gene is a member of Pax gene family and encodes a transcription factor associated with crucial roles in the development of pharyngeal pouch, scletrotome and limb bud. In zebrafish, the genome contains two Pax1 paralogs, DrPax1a and DrPax1b, which share high sequence similarity with other Pax1 genes. To elucidate the function of zebrafish DrPax1b gene, we first examined the gene expression pattern and found that it was mainly expressed in the endodermal pharyngeal pouch, caudal somites, notochord, and fin bud. Then, we performed knockdown experiments using antisense morpholino oligonucleotides, which lead to the defects in the vertebral column, tail, pharyngeal skeleton, and pectoral fin. Additionally, we also found that the mouse MmPax1 mRNA, but not the amphioxus AmphiPax1/9 mRNA, could rescue the MO-induced defects. Furthermore, sequence alignment revealed that the N-terminal region of vertebrate Pax1 and amphioxus Pax1/9 were highly conserved, whereas their C-terminal regions were relatively divergent. However, the chimeric Am(N)Dr(C)Pax1, Mm(N)Dr(C)Pax1 and Dr(N)Mm(C)Pax1 mRNA could partially rescue the defects, while the Dr(N)Am(C)Pax1 mRNA could not. In conclusion, our data demonstrate a conserved function of DrPax1b in the development of the vertebral column, pectoral fin and pharyngeal skeleton formation in zebrafish and also provide critical insight into the functional evolution of Pax1 gene by changing its C-terminal sequence.

Conserved localization of Pax6 and Pax7 transcripts in the brain of representatives of sarcopterygian vertebrates during development supports homologous brain regionalization

Many of the genes involved in brain patterning during development are highly conserved in vertebrates and similarities in their expression patterns help to recognize homologous cell types or brain regions. Among these genes, Pax6 and Pax7 are expressed in regionally restricted patterns in the brain and are essential for its development. In the present immunohistochemical study we analyzed the distribution of Pax6 and Pax7 cells in the brain of six representative species of tetrapods and lungfishes, the closest living relatives of tetrapods, at several developmental stages. The distribution patterns of these transcription factors were largely comparable across species. In all species only Pax6 was expressed in the telencephalon, including the olfactory bulbs, septum, striatum, and amygdaloid complex. In the diencephalon, Pax6 and Pax7 were distinct in the alar and basal parts, mainly in prosomeres 1 and 3. Pax7 specifically labeled cells in the optic tectum (superior colliculus) and Pax6, but not Pax7, cells were found in the tegmentum. Pax6 was found in most granule cells of the cerebellum and Pax7 labeling was detected in cells of the ventricular zone of the rostral alar plate and in migrated cells in the basal plate, including the griseum centrale and the interpeduncular nucleus. Caudally, Pax6 cells formed a column, whereas the ventricular zone of the alar plate expressed Pax7. Since the observed Pax6 and Pax7 expression patterns are largely conserved they can be used to identify subdivisions in the brain across vertebrates that are not clearly discernible with classical techniques.

Expression patterns of Pax6 and Pax7 in the adult brain of a urodele amphibian, Pleurodeles waltl

Expression patterns of Pax6, Pax7, and, to a lesser extent, Pax3 genes were analyzed by a combination of immunohistochemical techniques in the central nervous system of adult specimens of the urodele amphibian Pleurodeles waltl. Only Pax6 was found in the telencephalon, specifically the olfactory bulbs, striatum, septum, and lateral and central parts of the amygdala. In the diencephalon, Pax6 and Pax7 were distinct in the alar and basal parts, respectively, of prosomere 3. The distribution of Pax6, Pax7, and Pax3 cells correlated with the three pretectal domains. Pax7 specifically labeled cells in the dorsal mesencephalon, mainly in the optic tectum, and Pax6 cells were the only cells found in the tegmentum. Large populations of Pax7 cells occupied the rostral rhombencephalon, along with lower numbers of Pax6 and Pax3 cells. Pax6 was found in most granule cells of the cerebellum. Pax6 cells also formed a column of scattered neurons in the reticular formation and were found in the octavolateral area. The rhombencephalic ventricular zone of the alar plate expressed Pax7. Dorsal Pax7 cells and ventral Pax6 cells were found along the spinal cord. Our results show that the expression of Pax6 and Pax7 is widely maintained in the brains of adult urodeles, in contrast to the situation in other tetrapods. This discrepancy could be due to the generally pedomorphic features of urodele brains. Although the precise role of these transcription factors in adult brains remains to be determined, our findings support the idea that they may also function in adult urodeles.

Contrasting evolutionary dynamics of the developmental regulator PAX9, among bats, with evidence for a novel post-transcriptional regulatory mechanism

Morphological evolution can be the result of natural selection favoring modification of developmental signaling pathways. However, little is known about the genetic basis of such phenotypic diversity. Understanding these mechanisms is difficult for numerous reasons, yet studies in model organisms often provide clues about the major developmental pathways involved. The paired-domain gene, PAX9, is known to be a key regulator of development, particularly of the face and teeth. In this study, using a comparative genetics approach, we investigate PAX9 molecular evolution among mammals, focusing on craniofacially diversified (Phyllostomidae) and conserved (Vespertilionidae) bat families, and extend our comparison to other orders of mammal. Open-reading frame analysis disclosed signatures of selection, in which a small percentage of residues vary, and lineages acquire different combinations of variation through recurrent substitution and lineage specific changes. A few instances of convergence for specific residues were observed between morphologically convergent bat lineages. Bioinformatic analysis for unknown PAX9 regulatory motifs indicated a novel post-transcriptional regulatory mechanism involving a Musashi protein. This regulation was assessed through fluorescent reporter assays and gene knockdowns. Results are compatible with the hypothesis that the number of Musashi binding-elements in PAX9 mRNA proportionally regulates protein translation rate. Although a connection between morphology and binding element frequency was not apparent, results indicate this regulation would vary among craniofacially divergent bat species, but be static among conserved species. Under this model, Musashi's regulatory control of alternative human PAX9 isoforms would also vary. The presence of Musashi-binding elements within PAX9 of all mammals examined, chicken, zebrafish, and the fly homolog of PAX9, indicates this regulatory mechanism is ancient, originating basal to much of the animal phylogeny.

Identification, evolution and expression of a CD36 homolog in the basal chordate amphioxus Branchiostoma japonicum

CD36, as one member of scavenger receptor class B (SRB) family, is a transmembrane glycoprotein and has been associated with diverse normal physiological processes and pathological conditions. However, little is known about it in amphioxus, a model organism for insights into the origin and evolution of vertebrates. In this paper, CD36 homologs in amphioxus were identified. Evolutionary analysis suggested that amphioxus BfCD36F-a/b, which were more similar to vertebrate CD36, might represent the primitive form before the splitting of CD36, SRB1 and SRB2 genes during evolution. Then the BjCD36F-a cDNA was cloned from Branchiostoma japonicum using RACE technology. Real-time PCR and in situ hybridization revealed the expression of BjCD36F-a in all the tissues detected with the highest expression in the hepatic caecum. The BjCD36F-a expression was obviously up-regulated after feeding and down-regulated during fasting, indicating a role of BjCD36F-a in feeding regulation. Besides, the up-regulation expression of BjCD36F-a transcripts was also found after either Lipoteichoic acid (LTA) treatment in the BjCD36F-a-transfected FG cells or Escherichia coli (E. coli) challenge in vivo, implying an immune-related function for BjCD36F-a. Collectively, we identify and characterize a conserved gene that is important in the fundamental process of immune and nutritional regulation. These are the first such data in amphioxus, laying a foundation for further study of their physiological functions.

The function and developmental expression of alternatively spliced isoforms of amphioxus and Xenopus laevis Pax2/5/8 genes: revealing divergence at the invertebrate to vertebrate transition

Pax genes encode highly conserved transcription factors vital for metazoan development. Pax transcripts, particularly those in Group II (Pax2/5/8), are extensively alternatively spliced. This study compares the transcriptional activation capacity and developmental stage-specific expression of major isoforms of Group II Pax proteins in amphioxus (Branchiostoma floridae) and in Xenopus laevis. The comparison reveals considerable divergence of splice forms between the lineages, with the X. laevis Group II Pax genes (Pax2, Pax5, and Pax8) possessing a greater repertoire of regulated and functionally distinct splice forms than the single amphioxus gene (Pax2/5/8). Surprisingly, some apparently conserved splice forms are expressed at quite different levels during development in the two organisms and present different capacities to activate transcription. However, despite this divergence, the combinatorial transcriptional activation capacity of the isoforms present in early X. laevis and amphioxus development are broadly similar. This suggests that the some of the conserved functional roles, implied by the expression of Group II Pax genes in homologous tissues of amphioxus and X. laevis embryos, may depend upon the combination of isoforms expressed in a particular tissue at a particular time in development. Thus, during early development, the evolutionary constraint on the net effect of several isoforms co-expressed in a given tissue may be more strict than that on specific isoforms. This flexibility may facilitate the appearance of new exons and splicing patterns in the vertebrate duplicates, leading to isoforms with subtly distinct functions critical to the subsequent development of vertebrate-specific cell types and structures.

Systematic elucidation and in vivo validation of sequences enriched in hindbrain transcriptional control

Illuminating the primary sequence encryption of enhancers is central to understanding the regulatory architecture of genomes. We have developed a machine learning approach to decipher motif patterns of hindbrain enhancers and identify 40,000 sequences in the human genome that we predict display regulatory control that includes the hindbrain. Consistent with their roles in hindbrain patterning, MEIS1, NKX6-1, as well as HOX and POU family binding motifs contributed strongly to this enhancer model. Predicted hindbrain enhancers are overrepresented at genes expressed in hindbrain and associated with nervous system development, and primarily reside in the areas of open chromatin. In addition, 77 (0.2%) of these predictions are identified as hindbrain enhancers on the VISTA Enhancer Browser, and 26,000 (60%) overlap enhancer marks (H3K4me1 or H3K27ac). To validate these putative hindbrain enhancers, we selected 55 elements distributed throughout our predictions and six low scoring controls for evaluation in a zebrafish transgenic assay. When assayed in mosaic transgenic embryos, 51/55 elements directed expression in the central nervous system. Furthermore, 30/34 (88%) predicted enhancers analyzed in stable zebrafish transgenic lines directed expression in the larval zebrafish hindbrain. Subsequent analysis of sequence fragments selected based upon motif clustering further confirmed the critical role of the motifs contributing to the classifier. Our results demonstrate the existence of a primary sequence code characteristic to hindbrain enhancers. This code can be accurately extracted using machine-learning approaches and applied successfully for de novo identification of hindbrain enhancers. This study represents a critical step toward the dissection of regulatory control in specific neuronal subtypes.

Expression of the pair-rule gene homologs runt, Pax3/7, even-skipped-1 and even-skipped-2 during larval and juvenile development of the polychaete annelid Capitella teleta does not support a role in segmentation

Background

Annelids and arthropods each possess a segmented body. Whether this similarity represents an evolutionary convergence or inheritance from a common segmented ancestor is the subject of ongoing investigation.

Methods

To investigate whether annelids and arthropods share molecular components that control segmentation, we isolated orthologs of the Drosophila melanogaster pair-rule genes, runt, paired (Pax3/7) and eve, from the polychaete annelid Capitella teleta and used whole mount in situ hybridization to characterize their expression patterns.

Results

When segments first appear, expression of the single C. teleta runt ortholog is only detected in the brain. Later, Ct-runt is expressed in the ventral nerve cord, foregut and hindgut. Analysis of Pax genes in the C. teleta genome reveals the presence of a single Pax3/7 ortholog. Ct-Pax3/7 is initially detected in the mid-body prior to segmentation, but is restricted to two longitudinal bands in the ventral ectoderm. Each of the two C. teleta eve orthologs has a unique and complex expression pattern, although there is partial overlap in several tissues. Prior to and during segment formation, Ct-eve1 and Ct-eve2 are both expressed in the bilaterial pair of mesoteloblasts, while Ct-eve1 is expressed in the descendant mesodermal band cells. At later stages, Ct-eve2 is expressed in the central and peripheral nervous system, and in mesoderm along the dorsal midline. In late stage larvae and adults, Ct-eve1 and Ct-eve2 are expressed in the posterior growth zone.

Conclusions

C. teleta eve, Pax3/7 and runt homologs all have distinct expression patterns and share expression domains with homologs from other bilaterians. None of the pair-rule orthologs examined in C. teleta exhibit segmental or pair-rule stripes of expression in the ectoderm or mesoderm, consistent with an independent origin of segmentation between annelids and arthropods.