Evolutionary Conservation of Regulatory Elements in Vertebrate Hox Gene Clusters

Comparisons of DNA sequences among evolutionarily distantly related genomes permit identification of conserved functional regions in noncoding DNA. Hox genes are highly conserved in vertebrates, occur in clusters, and are uninterrupted by other genes. We aligned (PipMaker) the nucleotide sequences of the HoxA clusters of tilapia, pufferfish, striped bass, zebrafish, horn shark, human, and mouse, which are separated by approximately 500 million years of evolution. In support of our approach, several identified putative regulatory elements known to regulate the expression of Hox genes were recovered. The majority of the newly identified putative regulatory elements contain short fragments that are almost completely conserved and are identical to known binding sites for regulatory proteins (Transfac database). The regulatory intergenic regions located between the genes that are expressed most anteriorly in the embryo are longer and apparently more evolutionarily conserved than those at the other end of Hox clusters. Different presumed regulatory sequences are retained in either the Aα or Aβ duplicated Hox clusters in the fish lineages. This suggests that the conserved elements are involved in different gene regulatory networks and supports the duplication-deletion-complementation model of functional divergence of duplicated genes.     

Differential expression of neuroligin genes in the nervous system of zebrafish

The establishment and maturation of appropriate synaptic connections is crucial in the development of neuronal circuits. Cellular adhesion is believed to play a central role in this process. Neuroligins are neuronal cell adhesion molecules that are hypothesized to act in the initial formation and maturation of synaptic connections. In order to establish the zebrafish as a model to investigate the in vivo role of Neuroligin proteins in nervous system development, we identified the zebrafish orthologs of neuroligin family members and characterized their expression. Zebrafish possess seven neuroligin genes. Synteny analysis and sequence comparisons show that NLGN2, NLGN3, and NLGN4X are duplicated in zebrafish, but NLGN1 has a single zebrafish ortholog. All seven zebrafish neuroligins are expressed in complex patterns in the developing nervous system and in the adult brain. The spatial and temporal expression patterns of these genes suggest that they occupy a role in nervous system development and maintenance. Developmental Dynamics 239:703–714, 2010. © 2010 Wiley‐Liss, Inc.     

Functional Divergence of Two Zebrafish Midkine Growth Factors Following Fish-Specific Gene Duplication

In mammals, the unique midkine (mdk) gene encodes a secreted heparin-binding growth factor with neurotrophic activity. Here, we show the presence of two functional mdk genes named mdka and mdkb in zebrafish and rainbow trout. Both midkine proteins are clearly different from the related pleiotrophin, which was also identified in zebrafish and other fishes. Zebrafish mdka and mdkb genes map to linkage groups LG7 and LG25, respectively, both presenting synteny to human chromosome 11, in which the unique human ortholog mdk is located. At least four other genes unique in mammals are also present as duplicates on LG7 and LG25. Phylogenetic and divergence analyses suggested that LG7/LG25 paralogs including mdka and mdkb have been formed at approximately the same time, early during the evolution of the fish lineage. Hence, zebrafish and rainbow trout mdka and mdkb might have been generated by an ancient block duplication, and might be remnants of the proposed fish-specific whole-genome duplication. In contrast to the ubiquitous expression of their mammalian counterpart, zebrafish mdka and mdkb are expressed in spatially restricted, mostly nonoverlapping patterns during embryonic development and strongly in distinct domains in the adult brain. Ectopic ubiquitous expression of both mdk genes in early zebrafish embryos caused completely distinct effects on neural crest and floorplate development. These data indicate that mdka and mdkb underwent functional divergence after duplication. This provides an outstanding model to analyze the molecular mechanisms that lead to differences in pathways regulating the formation of homologous embryonic structures in different vertebrates.     

Repertoire and evolution of TNF superfamily in Crassostrea gigas: Implications for expansion and diversification of this superfamily in Mollusca

Tumor necrosis factor superfamily (TNFSF) members represent a group of cytokines participating in diverse immunological, pathological and developmental pathways. However, compared with deuterostomia and cnidaia, the composition and evolution of TNF homologous in protostomia are still not well understood. In the present study, a total of 81 TNF superfamily (TNFSF) genes from 15 mollusk species, including 23 TNFSF genes from Crassostrea gigas, were surveyed by genome-wide bioinformatics analysis. The phylogenetic analysis showed that 14 out of 23 C. gigas TNFSF genes in five clades exhibited orthologous relationships with Pinctada fucata TNFSF genes. Moreover, there were 15 C. gigas TNFSF genes located in oyster-specific clusters, which were contributed by small-scaled tandem and/or segmental duplication events in oyster. By comparing the sequences of duplicated TNFSF pairs, exon loss and variant in exon/intron length were revealed as the major modes of divergence in gene structure. Most of the duplicated C. gigas TNFSF pairs were evolved under purifying selection with consistent tissue expression patterns, implying functional constraint shaped diversification. This study demonstrated the expansion and early divergence of TNF superfamily in C. gigas, which provides potential insight into revealing the evolution and function of this superfamily in mollusk.     

Organization and expression analysis of the zebrafish hepcidin gene, an antimicrobial peptide gene conserved among vertebrates

Hepcidin is an antimicrobial peptide and iron-regulatory molecule that is conserved among vertebrates. Mutations or over-expression of the human hepcidin gene have been found in patients with hemochromatosis and refractory anemia. To further understand the function and regulation of hepcidin, animal models are needed. We sequenced cDNA, genes and upstream regions of zebrafish hepcidin and analyzed gene expression by kinetic PCR. Zebrafish hepcidin genes consist of two introns and three exons that encode a prepropeptide (91 amino acids). The amino acid sequences and gene organization were remarkably conserved between zebrafish and other species. Elevated gene expression was observed in abdominal organs, skin, and heart in fish that developed signs of infection following bacterial injection. Zebrafish may be a suitable model organism for further study of hepcidin gene regulation.     

15000 unique zebrafish EST clusters and their future use in microarray for profiling gene expression patterns during embryogenesis

A total of 15590 unique zebrafish EST clusters from two cDNA libraries have been identified. Most significantly, only 22% (3437) of the 15590 unique clusters matched 2805 (of 15200) clusters in the Danio rerio UniGene database, indicating that our EST set is complementary to the existing ESTs in the public database and will be invaluable in assisting the annotation of genes based on the upcoming zebrafish genome sequence. Blast search showed that 7824 of our unique clusters matched 6710 known or predicted proteins in the nonredundant database. A cDNA microarray representing approximately 3100 unique zebrafish cDNA clusters has been generated and used to profile the gene expression patterns across six different embryonic stages (cleavage, blastula, gastrula, segmentation, pharyngula, and hatching). Analysis of expression data using K-means clustering revealed that genes coding for muscle-specific proteins displayed similar expression patterns, confirming that the coordinate gene expression is important for myogenesis. Our results demonstrate that the combination of microarray technology with the zebrafish model system can provide useful information on how genes are coordinated in a genetic network to control zebrafish embryogenesis and can help to identify novel genes that are important for organogenesis.     

Discovery of an unusual alternative splicing pathway of the immunoglobulin heavy chain in a teleost fish, Danio rerio

In present study, we identified a novel membrane immunoglobulin M isotype from zebrafish (Danio rerio), which was designated as mIgM-2, adding a new member to the Immunoglobulin family in teleost fish. The full length of cloned mIgM-2 cDNA was 611 bp, encoding 150 amino acids. The putative mIgM-2 protein sequence consists of one constant region and a trans-membrane region. Phylogenetic analysis showed that mIgM-2 grouped with the known zebrafish IgM sequences. The mIgM-2 mRNA was widely expressed in immune-related tissues including intestine, kidney and skin. In vivo stimulation with LPS significantly up-regulates the expression of mIgM-2. Our results will add new insight into the immunoglobulin class diversity of teleost fish, and to better understand the evolutionary history of adaptive immunity from fish to mammals as a whole.     

A global assembly of cotton ESTs

Approximately 185,000 Gossypium EST sequences comprising >94,800,000 nucleotides were amassed from 30 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including drought stress and pathogen challenges. These libraries were derived from allopolyploid cotton (Gossypium hirsutum; A_T and D_T genomes) as well as its two diploid progenitors, Gossypium arboreum (A genome) and Gossypium raimondii (D genome). ESTs were assembled using the Program for Assembling and Viewing ESTs (PAVE), resulting in 22,030 contigs and 29,077 singletons (51,107 unigenes). Further comparisons among the singletons and contigs led to recognition of 33,665 exemplar sequences that represent a nonredundant set of putative Gossypium genes containing partial or full-length coding regions and usually one or two UTRs. The assembly, along with their UniProt BLASTX hits, GO annotation, and Pfam analysis results, are freely accessible as a public resource for cotton genomics. Because ESTs from diploid and allotetraploid Gossypium were combined in a single assembly, we were in many cases able to bioinformatically distinguish duplicated genes in allotetraploid cotton and assign them to either the A or D genome. The assembly and associated information provide a framework for future investigation of cotton functional and evolutionary genomics.     

Gene conversion and the evolution of protocadherin gene cluster diversity

The synaptic cell adhesion molecules encoded by the protocadherin gene cluster are hypothesized to provide a molecular code involved in the generation of synaptic complexity in the developing brain. Variation in copy number and sequence content of protocadherin cluster genes among vertebrate species could reflect adaptive differences in protocadherin function. We have completed an analysis of zebrafish protocadherin cluster genes. Zebrafish have two unlinked protocadherin clusters, DrPcdh1 and DrPcdh2. Like mammalian protocadherin clusters, DrPcdh1 has both α and γ variable and constant region exons. A consensus protocadherin promoter motif sequence identified in mammals is also conserved in zebrafish. Few orthologous relationships, however, are apparent between zebrafish and mammalian protocadherin proteins. Here we show that protocadherin cluster genes in human, mouse, rat, and zebrafish are subject to striking gene conversion events. These events are restricted to regions of the coding sequence, particularly the coding sequences of ectodomain 6 and the cytoplasmic domain. Diversity among paralogs is restricted to particular ectodomains that are excluded from conversion events. Conversion events are also strongly correlated with an increase in third-position GC content. We propose that the combination of lineage-specific duplication, restricted gene conversion, and adaptive variation in diversified ectodomains drives vertebrate protocadherin cluster evolution.     

Transcriptome profiling of zebrafish infected with Streptococcus suis

Streptococcus suis is an important pathogen in swine, and it also represents an emerging zoonotic agent. Zebrafish as a model for the evaluation of virulence of S. suis has been demonstrated before. Here, an Affymetrix Zebrafish GeneChip was used to identify alterations in gene expression of zebrafish injected with S. suis serotype 2 strain HA9801. The results showed that 189 genes were differentially expressed, of which 125 genes were upregulated and 64 genes were downregulated. Gene Ontology category and KEGG pathway were analyzed for differentially expressed genes. Upregulated genes were involved in response to bacterium, immune response, inflammatory response, complement activation, defense response. Three genes (encoding serum amyloid protein A, matrix metalloproteinase 9 and apoptosis-related cysteine protease) and genes involved in the regulation of IL-6 biosynthetic process, which have previously been implicated in the response to S. suis infection in other organisms, were also upregulated. Downregulated genes played roles in glycolysis, carbohydrate metabolic process, amino acids metabolism, behavior and muscle. The reliability of the data obtained from the microarray was verified by performing quantitative real-time PCR on 12 representative genes. The data may provide further validation of this model, which will contribute to understanding of S. suis pathogenic mechanisms.     

Close Split of Sorghum and Maize Genome Progenitors

It is generally believed that maize (Zea mays L. ssp. mays) arose as a tetraploid; however, the two progenitor genomes cannot be unequivocally traced within the genome of modern maize. We have taken a new approach to investigate the origin of the maize genome. We isolated and sequenced large genomic fragments from the regions surrounding five duplicated loci from the maize genome and their orthologous loci in sorghum, and then we compared these sequences with the orthologous regions in the rice genome. Within the studied segments, we identified 11 genes that were conserved in location, order, and orientation. We performed phylogenetic and distance analyses and examined the patterns of estimated times of divergence for sorghum and maize gene orthologs and also the time of divergence for maize orthologs. Our results support a tetraploid origin of maize. This analysis also indicates contemporaneous divergence of the ancestral sorghum genome and the two maize progenitor genomes about 11.9 million years ago (Mya). On the basis of a putative conversion event detected for one of the genes, tetraploidization must have occurred before 4.8 Mya, and therefore, preceded the major maize genome expansion by gene amplification and retrotransposition.     

Discovery of zebrafish (Danio rerio) interleukin-23 alpha (IL-23α) chain, a subunit important for the formation of IL-23, a cytokine involved in the development of Th17 cells and inflammation

This paper reports the cloning and sequencing of interleukin (IL)-23 p19 subunit for the first time within a non-mammalian species, the zebrafish (Danio rerio), which was discovered using a synteny approach. In addition, amino acid sequences were for IL-23 p19 subunits were also predicted from the stickleback, Fugu and Tetraodon genomes and included in this investigation. The zebrafish IL-23 p19 cDNA consisted of a 66 bp 5′ UTR, a 249 bp 3′ UTR and a single open reading frame of 567 bp giving a predicted 188 aa IL-23 p19 molecule. Multiple alignment of zebrafish IL-23 p19, with other known IL-23 p19 and IL-12 p35 amino acid sequences revealed areas of amino acid conservation, such as the presence of four predicted α-helixes, cysteines important for disulphide bond formation and the conservation of a tryptophan known to interact with the receptor. Amino acid homologies and phylogenetic analysis confirmed the relationship of the fish IL-23 p19 subunits with their mammalian homologues. All the teleost fish IL-23 p19 subunits were found to have 4 exons and 3 introns similar to that of human and mouse IL-23 p19 and a limited degree of synteny was found between the organisms for the regions containing the IL-23 p19 genes with only PAB-dependent poly(A)-specific ribonuclease subunit 2 (PAN2) and IL-23 p19 found in the same order on human chromosome 12 and all the fish genomes looked at. Lastly using real-time PCR, constitutive expression of IL-12 p40 and IL-23 p19 was observed in the kidney, liver, gut and muscle with IL-12 p40 expression higher than IL-23 p19. As soon as an hour after stimulation with LPS, there was an increase of IL-23 p19 in zebrafish leukocytes and an increase of IL-1β, IL-12 p40 and IL-23 p19 expression was found after infection of zebrafish for 1 or 6 days with Mycobaterium marinum strain E11.     

Duplicated major histocompatibility complex class II genes in the tongue sole (Cynoglossus semilaevis)

The major histocompatibility complex (MHC) molecule plays an important role in the vertebrate immune system. However, we have a limited understanding of the MHC genomic structure in teleosts. Using gene cloning and family analysis, we isolate the MHC class II genes in the tongue sole (Cynoglossus semilaevis) and find that both class II A and class II B genes are duplicated (named Cyse‐DAA and Cyse‐DBA, Cyse‐DAB and Cyse‐DBB, respectively). The class II A genes consist of four exons with a highly conserved genomic structure, but each gene has unique and defining exon 2 and intron 2 sequences. The class II B genes have a conserved six‐exon genomic structure, with intron 3 splitting the β2 encoding region into two exons. Each class II B gene has unique variations in exon 2 and intron 1 sequences. The two class II A genes have similar expression patterns among tissues, with high levels in spleen and gill. Both class II B genes have similar patterns, with high expression in spleen, gill and intestine. The alleles of MHC class II have wide distribution and reliable inheritance in the families analysed. This indicates that the duplicated MHC genes are all classical class II genes. The class II gene duplication with divergent exon and intron sequences, but similar expression patterns in tongue sole provides new insights into MHC evolution.     

Expression of strawberry notch family genes during zebrafish embryogenesis

Our previous study suggested a possible role for Sbno1, a mouse homologue of strawberry notch gene during brain development. In this report, we cloned the zebrafish homologues of sbno, and examined their expression pattern during embryogenesis by whole‐mount in situ hybridization. Zebrafish have three sbno genes: one Sbno1 homologue and two Sbno2 homologues, sbno2a and sbno2b. We observed that the expression of sbno1 and sbno2a was initially ubiquitous and gradually became predominant in the central nervous system as development progressed. The expression of sbno2b was observed in non‐neural tissues in contrast to the other two genes. sbno1 and sbno2a exhibited higher expression in distinct regions within the nervous system of pharyngula‐stage embryos, suggesting possible differing roles for sbno1 and sbno2a during later stages of embryogenesis. Together, the observed gene expression patterns suggest an important role of sbno‐family genes during development of the vertebrate central nervous system. Developmental Dynamics 239:1789–1796, 2010. © 2010 Wiley‐Liss, Inc.     

Expression of sox11 gene duplicates in zebrafish suggests the reciprocal loss of ancestral gene expression patterns in development.

To investigate the role of sox genes in vertebrate development, we have isolated sox11 from zebrafish (Danio rerio). Two distinct classes of sox11-related cDNAs were identified, sox11a and sox11b. The predicted protein sequences shared 75% identity. In a gene phylogeny, both sox11a and sox11b cluster with human, mouse, chick, and Xenopus Sox11, indicating that zebrafish, like Xenopus, has two orthologues of tetrapod Sox11. The work reported here investigates the evolutionary origin of these two gene duplicates and the consequences of their duplication for development. The sox11a and sox11b genes map to linkage groups 17 and 20, respectively, together with other loci whose orthologues are syntenic with human SOX11, suggesting that during the fish lineage, a large chromosome region sharing conserved syntenies with mammals has become duplicated. Studies in mouse and chick have shown that Sox11 is expressed in the central nervous system during development. Expression patterns of zebrafish sox11a and sox11b confirm that they are expressed in the developing nervous system, including the forebrain, midbrain, hindbrain, eyes, and ears from an early stage. Other sites of expression include the fin buds and somites. The two sox genes, sox11a and sox11b, are expressed in both overlapping and distinct sites. Their expression patterns suggest that sox11a and sox11b may share the developmental domains of the single Sox11 gene present in mouse and chick. For example, zebrafish sox11a is expressed in the anterior somites, and zebrafish sox11b is expressed in the posterior somites, but the single Sox11 gene of mouse is expressed in all the somites. Thus, the zebrafish duplicate genes appear to have reciprocally lost expression domains present in the sox11 gene of the last common ancestor of tetrapods and zebrafish. This splitting of the roles of Sox11 between two paralogues suggests that regulatory elements governing the expression of the sox11 gene in the common ancestor of zebrafish and tetrapods may have been reciprocally mutated in the zebrafish gene duplicates. This is consistent with duplicate gene evolution via a duplication-degeneration-complementation process.     

Effect of peanut addition to the cafeteria diet on adiposity and inflammation in zebrafish (Danio rerio)

Zebrafish (Danio rerio) emerged as a model of diet-induced obesity because of its genetic homology to humans. Peanut is a rich source of monounsaturated fatty acid and its consumption is associated with decreased inflammatory markers and obesity control. This study evaluated the effects of peanut addition to the cafeteria diet (CAF) by analysis of fatty acids into the head, adiposity and expression of TNF, IL6 and FASN genes using zebrafish as experimental model. The zebrafish were maintained in tanks for 60 days and treated with standard (ST) and CAF diets, respectively. The CAF diet increased the oleic acid content in zebrafish heads, however the body weight, body mass index, adipose tissue and expression of inflammatory and lipid metabolism genes did not differ between the groups. This study suggests that the addition of peanut in the CAF diet can control weight gain, the inflammatory markers and lipid metabolism in zebrafish model.