Comparative genomic analysis of teleost fish <Emphasis Type="Italic">bmal</Emphasis> genes

Bmal1 (Brain and muscle ARNT like 1) gene is a key circadian clock gene. Tetrapods also have the second Bmal gene, Bmal2. Fruit fly has only one bmal1/cycle gene. Interrogation of the five teleost fish genome sequences coupled with phylogenetic and splice site analyses found that zebrafish have two bmal1 genes, bmal1a and bmal1b, and bmal2a; Japanese pufferfish (fugu), green spotted pufferfish (tetraodon) and Japanese medaka fish each have two bmal2 genes, bmal2a and bmal2b, and bmal1a; and three-spine stickleback have bmal1a and bmal2b. Syntenic analysis further indicated that zebrafish bmal1a/bmal1b, and fugu, tetraodon and medaka bmal2a/bmal2b are ancient duplicates. Although the dN/dS ratios of these four fish bmal duplicates are all <1, implicating they have been under purifying selection, the Tajima relative rate test showed that fugu, tetraodon and medaka bmal2a/bmal2b have asymmetric evolutionary rates, suggesting that one of these duplicates have been subject to positive selection or relaxed functional constraint. These results support the notion that teleost fish bmal genes were derived from the fish-specific genome duplication (FSGD), divergent resolution following the duplication led to retaining different ancient bmal duplicates in different fishes, which could have shaped the evolution of the complex teleost fish timekeeping mechanisms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Twist controls skeletal development and dorsoventral patterning by regulating runx2 in zebrafish

Background

Twist1a and twist1b are the principal components of twists that negatively regulate a number of cellular signaling events. Expression of runx2 and downstream targets is essential for skeletal development and ventral organizer formation and specification in early vertebrate embryos, but what controls ventral activity of maternal runx2 and how twists function in zebrafish embryogenesis still remain unclear.

Methodology/Principal Findings

By studying the loss of twist induced by injection of morpholino-oligonucleotide in zebrafish, we found that twist1a and twist1b, but not twist2 or twist3, were required for proper skeletal development and dorsoventral patterning in early embryos. Overexpression of twist1a or twist1b following mRNA injection resulted in deteriorated skeletal development and formation of typical dorsalized embryos, whereas knockdown of twist1a and twist1b led to the formation of abnormal embryos with enhanced skeletal formation and typical ventralized patterning. Overexpression of twist1a or twist1b decreased the expression of runx2b, whereas twist1a and twist1b knockdown increased runx2b expression. We have further demonstrated that phenotypes induced by twist1a and twist1b knockdown were rescued by runx2b knockdown.

Conclusions/Significance

Together, these results suggest that twist1a and twist1b control skeletal development and dorsoventral patterning by regulating runx2b in zebrafish and provide potential targets for the treatment of diseases or syndromes associated with decreased skeletal development.     

Characterization of <Emphasis Type="Italic">twist</Emphasis> and <Emphasis Type="Italic">snail</Emphasis> gene expression during mesoderm and nervous system development in the polychaete annelid <Emphasis Type="Italic">Capitella</Emphasis> sp. I

To investigate the evolutionary history of mesoderm in the bilaterian lineage, we are studying mesoderm development in the polychaete annelid, Capitella sp. I, a representative lophotrochozoan. In this study, we focus on the Twist and Snail families as candidate mesodermal patterning genes and report the isolation and in situ expression patterns of two twist homologs (CapI-twt1 and CapI-twt2) and two snail homologs (CapI-sna1 and CapI-sna2) in Capitella sp. I. CapI-twt1 is expressed in a subset of mesoderm derivatives during larval development, while CapI-twt2 shows more general mesoderm expression at the same stages. Neither twist gene is detected before the completion of gastrulation. The two snail genes have very distinct expression patterns. At cleavage and early gastrula stages, CapI-sna1 is broadly expressed in precursors of all three germ layers and becomes restricted to cells around the closing blastopore during late gastrulation; CapI-sna2 expression is not detected at these stages. After gastrulation, both snail genes are expressed in the developing central nervous system (CNS) at stages when neural precursor cells are internalized, and CapI-sna1 is also expressed laterally within the segmental mesoderm. Based on the expression patterns in this study, we suggest a putative function for Capitella sp. I twist genes in mesoderm differentiation and for snail genes in regulating CNS development and general cell migration during gastrulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Natural selection in the <Emphasis Type="Italic">TLR-</Emphasis>related genes in the course of primate evolution

The innate immune system constitutes the front line of host defense against pathogens. Toll-like receptors (TLRs) recognize molecules derived from pathogens and play crucial roles in the innate immune system. Here, we provide evidence that the TLR-related genes have come under natural selection pressure in the course of primate evolution. We compared the nucleotide sequences of 16 TLR-related genes, including TLRs (TLR1–10), MYD88, TILAP, TICAM1, TICAM2, MD2, and CD14, among seven primate species. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving regions in the TLR-related genes. The genomic segments encoding the intracellular Toll/interleukin 1 receptor domains, which exhibited lower rates of non-synonymous substitution, have undergone purifying selection. In contrast, TLR4, which carried a high proportion of non-synonymous substitutions in the part of extracellular domain spanning 200 amino acids, was found to have been the suggestive target of positive Darwinian selection in primate evolution. However, sequence analyses from 25 primate species, including eight hominoids, six Old World monkeys, eight New World monkeys, and three prosimians, showed no evidence that the pressure of positive Darwinian selection has shaped the pattern of sequence variations in TLR4 among New World monkeys and prosimians. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Molecular cloning and gene expression of the prox1a and prox1b genes in the medaka,Oryzias latipes

Prox1 is a prospero-related homeobox gene. Prox1 is expressed in various internal organs and is related to those differentiations. Small fishes such as the zebrafish and the medaka are useful model animals in the clarification of the mechanism of development. The zebrafish prox1 is also identified, and it contributes to clarifying the function of prox1. However, it is necessary to note that many genes are duplicated in teleost fishes. In this study, we identified the orthologs of the mammalian prox1 gene in the medaka. The gene was also duplicated in the medaka, and we named it prox1a and prox1b. In silico analysis from the perspective of synteny indicated that medaka prox1a was similar to the prox1 gene of other vertebrates. Medaka prox1a was expressed in all internal organs that we have examined by RT-PCR. In contrast, medaka prox1b expression was limited to the brain, heart, liver, kidney, thymus, gill, testis, and ovary. This suggests that the two prox1 genes do not have a complementary relationship. In addition, we examined their expression patterns during embryonic development using whole-mount in situ hybridization. The expression pattern of prox1a showed a pattern similar to that of zebrafish prox1. In contrast, medaka prox1b was expressed asymmetrically in part of the central nervous system, especially strongly in the right side of the habenula.     

Use of mitochondrial and nuclear genes to infer the origin of two endemic pigeons from the Canary Islands

DNA nucleotide sequences from two mitochondrial genes (cytochrome b and NADH dehydrogenase subunit 2) and the nuclear intron 7 of β-fibrinogen were obtained to infer the phylogenetic origin of the two endemic Canarian pigeons: Bolle’s Pigeon (Columba bollii) and Laurel Pigeon (C. junoniae). Phylogenetic analyses of mitochondrial and nuclear genes based on maximum parsimony, maximum likelihood and Bayesian inference all converged into a congruent topology: C. bollii clusters together with the Wood Pigeon (C. palumbus) which is common in Europe and Asia, while C. junoniae was found near the base of the clade that includes other species of the genus Columba from the Old World. Laurel Pigeon probably represents an old lineage that might have colonized the Canary Islands a long time ago (20 My) while Bolle’s Pigeon might have arrived on the archipelago much later during the Upper Miocene (5 My). Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative analysis of her genes during fish somitogenesis suggests a mouse/chick-like mode of oscillation in medaka

Somitogenesis is the key developmental step, which divides the vertebrate body axis into segmentally repeated structures. It requires an intricate process of pre-patterning, which is driven by an oscillator mechanism consisting of the Delta–Notch pathway and various hairy- and Enhancer of split-related (her) genes. The subset of her genes, which are necessary to set up the segmentation clock, reveal a complex scenario of interactions. To understand which her genes are essential core players in this process, we compared the expression patterns of somitogenesis-relevant her genes in zebrafish and medaka (Oryzias latipes). Most of the respective medaka genes (Ol-her) are duplicated like what has been shown for zebrafish (Dr-her) and pufferfish genes (Fr-her). However, zebrafish genes show some additional copies and significant differences in expression patterns. For the paralogues Dr-her1 and Dr-her11, only one copy exists in the medaka (Ol-her1/11), which combines the expression patterns found for both zebrafish genes. In contrast to Dr-her5, the medaka orthologue appears to play a role in somitogenesis because it is expressed in the presomitic mesoderm (PSM). PSM expression also suggests a role for both Ol-her13 genes, homologues of mouse Hes6 (mHes6), in this process, which would be consistent with a conserved mHes6 homologue gear in the segmentation clock exclusively in lower vertebrates. Members of the mHes5 homologue group seem to be involved in somite formation in all vertebrates (e.g. Dr- and Ol-her12), although different paralogues are additionally recruited in zebrafish (e.g. Dr-her15) and medaka (e.g. Ol-her4). We found that the linkage between duplicates is strongly conserved between pufferfish and medaka and less well conserved in zebrafish. Nevertheless, linkage and orientation of several her duplicates are identical in all three species. Therefore, small-scale duplications must have happened before whole genome duplication occurred in a fish ancestor. Expression of multiple stripes in the intermediate PSM, characteristic for the zebrafish orthologues, is absent in all somitogenesis-related her genes of the medaka. In fact, the expression mode of Ol-her1/11 and Ol-her5 indicates dynamism similar to the hairy clock genes in chicken and mouse. This suggests that Danio rerio shows a rather derived clock mode when compared to other fish species and amniotes or that, alternatively, the clock mode evolved independently in zebrafish, medaka and mouse or chicken.An erratum to this article can be found at     

Comparative genomic analysis of the major histocompatibility complex class I region in the teleost genus <Emphasis Type="Italic">Oryzias</Emphasis>

The major histocompatibility complex (MHC) class I region of teleosts harbors a tight cluster of the class IA genes and several other genes directly involved in class I antigen presentation. Moreover, the dichotomous haplotypic lineages (termed d- and N- lineages) of the proteasome subunit beta genes, PSMB8 and PSMB10, are present in this region of the medaka, Oryzias latipes. To understand the evolution of the Oryzias MHC class I region at the nucleotide sequence level, we analyzed bacterial artificial chromosome clones covering the MHC class I region containing the d- lineage of Oryzias luzonensis and the d- and N- lineages of Oryzias dancena. Comparison among these three elucidated sequences and the published sequences of the d- and N- lineages of O. latipes indicated that the order and orientation of the encoded genes were completely conserved among these five genomic regions, except for the class IA genes, which showed species-specific variation in copy number. The PSMB8 and PSMB10 genes showed trans-species dimorphism. The remaining regions flanking the PSMB10, PSMB8, and class IA genes showed high degrees of sequence conservation at interspecies as well as intraspecies levels. Thus, the three independent evolutionary patterns under apparently distinctive selective pressures are recognized in the Oryzias MHC class I region. Electronic Supplementary Material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Duplicate <Emphasis Type="Italic">dmbx1</Emphasis>genes regulate progenitor cell cycle and differentiation during zebrafish midbrain and retinal development

Background  

The Dmbx1 gene is important for the development of the midbrain and hindbrain, and mouse gene targeting experiments reveal that this gene is required for mediating postnatal and adult feeding behaviours. A single Dmbx1 gene exists in terrestrial vertebrate genomes, while teleost genomes have at least two paralogs. We compared the loss of function of the zebrafish dmbx1a and dmbx1b genes in order to gain insight into the molecular mechanism by which dmbx1 regulates neurogenesis, and to begin to understand why these duplicate genes have been retained in the zebrafish genome.     

<Emphasis Type="Italic">kitb</Emphasis>, a second zebrafish ortholog of mouse <Emphasis Type="Italic">Kit</Emphasis>

The large numbers of duplicated pairs of genes in zebrafish compared to their mammalian counterparts has lead to the notion that expression of zebrafish co-orthologous pairs in some cases can together describe the expression of their mammalian counterpart. Here, we explore this notion by identification and analysis of a second zebrafish ortholog of the mammalian Kit receptor tyrosine kinase (kitb). We show that in embryos, kitb is expressed in a non-overlapping pattern to that of kita, in the anterior ventral mesoderm, Rohon-beardRohon–Beard neurons, the otic vesicle, and trigeminal ganglia. The expression pattern of kita and kitb in zebrafish together approximates that of Kit in mouse, with the exception that neither zebrafish kit gene is expressed in primordial germ cells, a site of kit expression in the mouse embryo. In addition, zebrafish kita is expressed in a site of zebrafish primitive hematopoiesis but not required for blood development, and we fail to detect kitb expression in sites of zebrafish hematopoiesis. Thus, the expression and function of zebrafish kit genes cannot be described as a simple partition of the expression and function of mouse Kit. We discuss the possibility that these unaccounted for expression domains and functions are derived from more ancestral gene duplications and partitioning instead of the relatively recent teleost teleost-specific duplication. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Unexpected functional redundancy between Twist and Slug (Snail2) and their feedback regulation of NF-κB via Nodal and Cerberus

A NF-κB-Twist-Snail network controls axis and mesoderm formation in Drosophila. Using translation-blocking morpholinos and hormone-regulated proteins, we demonstrate the presence of an analogous network in the early Xenopus embryo. Loss of twist (twist1) function leads to a reduction of mesoderm and neural crest markers, an increase in apoptosis, and a decrease in snail1 (snail) and snail2 (slug) mRNA levels. Injection of snail2 mRNA rescues twist's loss of function phenotypes and visa versa. In the early embryo NF-κB/RelA regulates twist, snail2, and snail1 mRNA levels; similarly Nodal/Smad2 regulate twist, snail2, snail1, and relA RNA levels. Both Twist and Snail2 negatively regulate levels of cerberus RNA, which encodes a Nodal, bone morphogenic protein (BMP), and Wnt inhibitor. Cerberus's anti-Nodal activity inhibits NF-κB activity and decreases relA RNA levels. These results reveal both conserved and unexpected regulatory interactions at the core of a vertebrate's mesodermal specification network.     

Evolution of teleostean hatching enzyme genes and their paralogous genes

We isolated genes for hatching enzymes and their paralogs having two cysteine residues at their N-terminal regions in addition to four cysteines conserved in all the astacin family proteases. Genes for such six-cysteine-containing astacin proteases (C6AST) were searched out in the medaka genome database. Five genes for MC6AST1 to 5 were found in addition to embryo-specific hatching enzyme genes. RT-PCR and whole-mount in situ hybridization evidenced that MC6AST1 was expressed in embryos and epidermis of almost all adult tissues examined, while MC6AST2 and 3 were in mesenterium, intestine, and testis. MC6AST4 and 5 were specifically expressed in jaw. In addition, we cloned C6AST cDNA homologs from zebrafish, ayu, and fugu. The MC6AST1 to 5 genes were classified into three groups in the phylogenetic positions, and the expression patterns and hatching enzymes were clearly discriminated from other C6ASTs. Analysis of the exon–intron structures clarified that genes for hatching enzymes MHCE and MAHCE were intron-less, while other MC6AST genes were basically the same as the gene for another hatching enzyme MLCE. In the basal Teleost, the C6AST genes having the ancestral exon–intron structure (nine exon/eight intron structure) first appeared by duplication and chromosomal translocation. Thereafter, maintaining such ancestral exon–intron structure, the LCE gene was newly diversified in Euteleostei, and the MC6AST1 to 5 gene orthologs were duplicated and diversified independently in respective fish lineages. The HCE gene lost all introns in Euteleostei, whereas in the lineage to zebrafish, it was translocated from chromosome to chromosome and lost some of its introns.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.The nucleotide sequence data reported in the present paper will appear in the DDBJ/EMBL/GenBank nucleotide sequence databases with accession numbers from AB256940 to AB256952.     

Comparative evolution of the mitochondrial cytochrome <Emphasis Type="Italic">b</Emphasis> gene and nuclear <Emphasis Type="Italic">S7</Emphasis> ribosomal protein gene intron 1 in sinipercid fishes and their relatives

The extant sinipercids are a group of freshwater percoid fishes endemic to East Asia. A recent mitochondrial cytochrome b phylogeny of sinipercids has challenged some aspects for their traditional taxonomy and molecular phylogeny, especially for the monophyly of Sinipercidae. In this study, we analyzed mitochondrial cytochrome b and nuclear encoded S7 ribosomal protein gene intron 1 for 10 sinipercid species and 11 related species to compare the phylogenetic signal and nucleotide substitution properties of these two gene sequences. The length of S7 intron 1 ranged from 461 to 719 bp, but alignment was not difficult, and the indels, the proportion of which in the total nucleotides ranged from 3.76 to 45.83%, were phylogenetically informative. Our results indicate that: (1) the relative rate presented by cyt b is five times that of S7 intron 1; (2) the proportion of phylogenetic information is higher in S7 than in cyt b; (3) S7 intron 1 has more base composition bias, but more uniform nucleotide substitution properties; (4) the overall ratio between transitions and transversions in S7 intron 1 is lower than in cyt b. Maximum parsimony and Bayesian analyses of aligned S7 intron 1 and the combined S7 and cyt b dataset resulted in phylogenies that contained the previously identified genera Siniperca and Coreoperca, whereas the monophyly of Coreoperca cannot be corroborated by separate cyt b analysis. The monophyly of Sinipercidae is not supported in separate and combined dataset analyses, although the alternative hypothesis cannot be significantly rejected based on approximately unbiased tests and Shimodaira–Hasegawa tests. Overall, maximum parsimony analyses result in trees with a lack of phylogenetic resolution in deep nodes, and the signal from S7 intron 1 conflicts the cyt b signal in the combined dataset analyses. The reasons for the poor performance of cyt b to S7 intron 1 in the phylogeny are discussed.     

<Emphasis Type="Italic">Dwarf 88</Emphasis>, a novel putative esterase gene affecting architecture of rice plant

Rice architecture is an important agronomic trait that affects grain yield. We characterized a tillering dwarf mutant d88 derived from Oryza sativa ssp. japonica cultivar Lansheng treated with EMS. The mutant had excessive shorter tillers and smaller panicles and seeds compared to the wild-type. A reduction in number and size of parenchyma cells around stem marrow cavity as well as a delay in the elongation of parenchyma cells caused slender tillers and dwarfism in the d88 mutant. The D88 gene was isolated via map-based cloning and identified to encode a putative esterase. The gene was expressed in most rice organs, with especially high levels in the vascular tissues. The mutant carried a nucleotide substitution in the first exon of the gene that led to the substitution of arginine for glycine, which presumably disrupted the functionally conserved N-myristoylation domain of the protein. The function of the gene was confirmed by complementation test and antisense analysis. D88, thus, represents a new category of genes that regulates cell growth and organ development and consequently plant architecture. The potential relationship between the tiller formation associated genes and D88 is discussed and future identification of the substrate for D88 may lead to the characterization of new pathways regulating plant development. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Cytochrome-<Emphasis Type="Italic">b</Emphasis> variation in <Emphasis Type="Italic">Apis mellifera</Emphasis> samples and its association with COI–COII patterns

Five Mbo I (Mbo-A, Mbo-M, Mbo-C¹, Mbo-C² and Mbo-C³) and Hinf I (Hinf-1 to Hinf-5) patterns were observed in Apis mellifera samples after restriction of a 485 bp fragment of the mitochondrial cytochrome-b (cyt-b) gene. Associating the cyt-b Restriction fragment length polymorphism (RFLP) pattern of each sample to its respective previously established COI–COII (Dra I sites) pattern, five restriction patterns (Mbo-C¹, Mbo-C², Mbo-C³, Hinf-1 and Hinf-4) were observed in samples of maternal origin associated to the evolutionary branch C. No deletions or insertions were observed and the nucleotide substitution rate was estimated at 5.4%. Higher nucleotide diversity was observed among the branch C-haplotypes when compared with A and M lineages. Further studies are needed to confirm if the cyt-b + COI–COII haplotypes help to assign certain phylogeographic patterns to the branch C and to clarify phylogenetic relationships among A. mellifera subspecies.