Identification, characterization, and mapping of expressed sequence tags from an embryonic zebrafish heart cDNA library

Dynamics of gene expression patterns during early development of the European seabass (Dicentrarchus labrax)

Larval and embryonic stages are the most critical period in the life cycle of marine fish. Key developmental events occur early in development and are influenced by external parameters like stress, temperature, salinity, and photoperiodism. Any failure may cause malformations, developmental delays, poor growth, and massive mortalities. Advanced understanding of molecular processes underlying marine larval development may lead to superior larval rearing conditions. Today, the new sequencing and bioinformatic methods allow transcriptome screens comprising messenger (mRNA) and microRNA (miRNA) with the scope of detecting differential expression for any species of interest. In the present study, we applied Illumina technology to investigate the transcriptome of early developmental stages of the European seabass ( Dicentrarchus labrax). The European seabass, in its natural environment, is a euryhaline species and has shown high adaptation processes in early life phases. During its embryonic and larval phases the European seabass lives in a marine environment and as a juvenile it migrates to coastal zones, estuaries, and lagoons. Investigating the dynamics of gene expression in its early development may shed light on factors promoting phenotypic plasticity and may also contribute to the improvement and advancement of rearing methods of the European seabass, a species of high economic importance in European and Mediterranean aquaculture. We present the identification, characterization, and expression of mRNA and miRNA, comprising paralogous genes and differentially spliced transcripts from early developmental stages of the European seabass. We further investigated the detection of possible interactions of miRNA with mRNA.

Cardiac transcriptome and dilated cardiomyopathy genes in zebrafish

BACKGROUND

Genetic studies of cardiomyopathy and heart failure have limited throughput in mammalian models. Adult zebrafish have been recently pursued as a vertebrate model with higher throughput, but genetic conservation must be tested.

METHODS AND RESULTS

We conducted transcriptome analysis of zebrafish heart and searched for fish homologues of 51 known human dilated cardiomyopathy-associated genes. We also identified genes with high cardiac expression and genes with differential expression between embryonic and adult stages. Among tested genes, 30 had a single zebrafish orthologue, 14 had 2 homologues, and 5 had ≥3 homologues. By analyzing the expression data on the basis of cardiac abundance and enrichment hypotheses, we identified a single zebrafish gene for 14 of 19 multiple-homologue genes and 2 zebrafish homologues of high priority for ACTC1. Of note, our data suggested vmhc and vmhcl as functional zebrafish orthologues for human genes MYH6 and MYH7, respectively, which are established molecular markers for cardiac remodeling.

CONCLUSIONS

Most known genes for human dilated cardiomyopathy have a corresponding zebrafish orthologue, which supports the use of zebrafish as a conserved vertebrate model. Definition of the cardiac transcriptome and fetal gene program will facilitate systems biology studies of dilated cardiomyopathy in zebrafish.

Gene expression throughout a vertebrate's embryogenesis

Background

Describing the patterns of gene expression during embryonic development has broadened our understanding of the processes and patterns that define morphogenesis. Yet gene expression patterns have not been described throughout vertebrate embryogenesis. This study presents statistical analyses of gene expression during all 40 developmental stages in the teleost Fundulus heteroclitus using four biological replicates per stage.

Results

Patterns of gene expression for 7,000 genes appear to be important as they recapitulate developmental timing. Among the 45% of genes with significant expression differences between pairs of temporally adjacent stages, significant differences in gene expression vary from as few as five to more than 660. Five adjacent stages have disproportionately more significant changes in gene expression (> 200 genes) relative to other stages: four to eight and eight to sixteen cell stages, onset of circulation, pre and post-hatch, and during complete yolk absorption. The fewest differences among adjacent stages occur during gastrulation. Yet, at stage 16, (pre-mid-gastrulation) the largest number of genes has peak expression. This stage has an over representation of genes in oxidative respiration and protein expression (ribosomes, translational genes and proteases). Unexpectedly, among all ribosomal genes, both strong positive and negative correlations occur. Similar correlated patterns of expression occur among all significant genes.

Conclusions

These data provide statistical support for the temporal dynamics of developmental gene expression during all stages of vertebrate development.

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.

Comparative analysis of the testis and ovary transcriptomes in zebrafish by combining experimental and computational tools

Studies on the zebrafish model have contributed to our understanding of several important developmental processes, especially those that can be easily studied in the embryo. However, our knowledge on late events such as gonad differentiation in the zebrafish is still limited. Here we provide an analysis on the gene sets expressed in the adult zebrafish testis and ovary in an attempt to identify genes with potential role in (zebra)fish gonad development and function. We produced 10,533 expressed sequence tags (ESTs) from zebrafish testis or ovary and downloaded an additional 23,642 gonad-derived sequences from the zebrafish EST database. We clustered these sequences together with over 13,000 kidney-derived zebrafish ESTs to study partial transcriptomes for these three organs. We searched for genes with gonad-specific expression by screening macroarrays containing at least 2600 unique cDNA inserts with testis-, ovary- and kidney-derived cDNA probes. Clones hybridizing to only one of the two gonad probes were selected, and subsequently screened with computational tools to identify 72 genes with potentially testis-specific and 97 genes with potentially ovary-specific expression, respectively. PCR-amplification confirmed gonad-specificity for 21 of the 45 clones tested (all without known function). Our study, which involves over 47,000 EST sequences and specialized cDNA arrays, is the first analysis of adult organ transcriptomes of zebrafish at such a scale. The study of genes expressed in adult zebrafish testis and ovary will provide useful information on regulation of gene expression in teleost gonads and might also contribute to our understanding of the development and differentiation of reproductive organs in vertebrates.

Identification of a molecular marker for type A spermatogonia by microarray analysis using gonadal cells from pvasa-GFP transgenic rainbow trout (Oncorhynchus mykiss)

The spermatogonia of fish can be classified as being either undifferentiated type A spermatogonia or differentiated type B spermatogonia. Although type A spermatogonia, which contain spermatogonial stem cells, have been demonstrated to be a suitable material for germ cell transplantation, no molecular markers for distinguishing between type A and type B spermatogonia in fish have been developed to date. We therefore sought to develop a molecular marker for type A spermatogonia in rainbow trout. Using GFP-dependent flow cytometry (FCM), enriched fractions of type A and type B spermatogonia, testicular somatic cells, and primordial germ cells were prepared from rainbow trout possessing the green fluorescent protein (GFP) gene driven by trout vasa regulatory regions (pvasa-GFP rainbow trout). The gene-expression profiles of each cell fraction were then compared with a microarray containing cDNAs representing 16,006 genes from several salmonid species. Genes exhibiting high expression for type A spermatogonia relative to above-mentioned other types of gonadal cells were identified and subjected to RT-PCR and quatitative PCR analysis. Since only the rainbow trout notch1 homologue showed significantly high expression in the type A spermatogonia-enriched fraction, we propose that notch1 may be a useful molecular marker for type A spermatogonia. The combination of GFP-dependent FCM and microarray analysis of pvasa-GFP rainbow trout can therefore be applied to the identification of potentially useful molecular markers of germ cells in fish.

Characterisation of purified parvalbumin from five fish species and nucleotide sequencing of this major allergen from Pacific pilchard, Sardinops sagax

IgE-mediated allergic reaction to seafood is a common cause of food allergy including anaphylactic reactions. Parvalbumin, the major fish allergen, has been shown to display IgE cross-reactivity among fish species consumed predominantly in Europe and the Far East. However, cross-reactivity studies of parvalbumin from fish species widely consumed in the Southern hemisphere are limited as is data relating to immunological and molecular characterisation. In this study, antigenic cross-reactivity and the presence of oligomers and isomers of parvalbumin from five highly consumed fish species in Southern Africa were assessed by immunoblotting using purified parvalbumin and crude fish extracts. Pilchard (Sardinops sagax) parvalbumin was found to display the strongest IgE reactivity among 10 fish-allergic consumers. The cDNA sequence of the beta-form of pilchard parvalbumin was determined and designated Sar sa 1.0101 (accession number FM177701 EMBL/GenBank/DDBJ databases). Oligomeric forms of parvalbumin were observed in all fish species using a monoclonal anti-parvalbumin antibody and subject's sera. Isoforms varied between approximately 10-13 kDa. A highly cross-reactive allergenic isoform of parvalbumin was identified and sequenced, providing a successful primary step towards the generation of a recombinant form that could be used for diagnostic and potential therapeutic use in allergic individuals.

Transcriptomic and proteomic analysis of global ischemia and cardioprotection in the rabbit heart

Cardioplegia is used to partially alleviate the effects of surgically induced global ischemia injury; however, the molecular mechanisms involved in this cardioprotection remain to be elucidated. To improve the understanding of the molecular processes modulating the effects of global ischemia and the cardioprotection afforded by cardioplegia, we constructed rabbit heart cDNA libraries and isolated, sequenced, and identified a compendium of nonredundant cDNAs for use in transcriptomic and proteomic analyses. New Zealand White rabbits were used to compare the effects of global ischemia and cardioplegia compared with control (nonischemic) hearts. The effects of RNA and protein synthesis on the cardioprotection afforded by cardioplegia were investigated separately by preperfusion with either alpha-amanitin or cycloheximide. Our results demonstrate that cardioplegia partially ameliorates the effects of global ischemia and that the cardioprotection is modulated by RNA- and protein-dependent mechanisms. Transcriptomic and proteomic enrichment analyses indicated that global ischemia downregulated genes/proteins associated with mitochondrial function and energy production, cofactor catabolism, and the generation of precursor metabolites of energy. In contrast, cardioplegia significantly increased differentially expressed genes/proteins associated with the mitochondrion and mitochondrial function and significantly upregulated the biological processes of muscle contraction, involuntary muscle contraction, carboxylic acid and fatty acid catabolic processes, fatty acid beta-oxidation, and fatty acid metabolic processes.

Identification and analysis of a prepro-chicken gonadotropin releasing hormone II (preprocGnRH-II) precursor in the Asian seabass, Lates calcarifer, based on an EST-based assessment of its brain transcriptome

Using a novel library of 5637 expressed sequence tags (ESTs) from the brain tissue of the Asian seabass (Lates calcarifer), we first characterized the brain transcriptome for this economically important species. The ESTs generated from the brain of L. calcarifer yielded 2410 unique transcripts (UTs) which comprise of 982 consensi and 1428 singletons. Based on database similarity, 1005 UTs (41.7%) can be assigned putative functions and were grouped into 12 functional categories related to the brain function. Amongst others, we have identified genes that are putatively involved in energy metabolism, ion pumps and channels, synapse related genes, neurotransmitter and its receptors, stress induced genes and hormone related genes. Subsequently we selected a putative preprocGnRH-II precursor for further characterization. The complete cDNA sequence of the gene obtained was found to code for an 85-amino acid polypeptide that significantly matched preprocGnRH-II precursor sequences from other vertebrates, and possesses structural characteristics that are similar to that of other species, consisting of a signal peptide (23 residues), a GnRH decapeptide (10 residues), an amidation/proteolytic-processing signal (glycine-lysine-argine) and a GnRH associated peptide (GAP) (49 residues). Phylogenetic analysis showed that this putative L. calcarifer preprocGnRH-II sequence is a member of the subcohort Euteleostei and divergent from the sequences of the subcohort Otocephalan. These findings provide compelling evidence that the putative L. calcarifer preprocGnRH-II precursor obtained in this study is orthologous to that of other vertebrates. The functional prediction of this preprocGnRH-II precursor sequence through in silico analyses emphasizes the effectiveness of the EST approach in gene identification in L. calcarifer.

A method for the construction of equalized directional cDNA libraries from hydrolyzed total RNA

BACKGROUND

The transcribed sequences of a cell, the transcriptome, represent the trans-acting fraction of the genetic information, yet eukaryotic cDNA libraries are typically made from only the poly-adenylated fraction. The non-coding or translated but non-polyadenylated RNAs are therefore not represented. The goal of this study was to develop a method that would more completely represent the transcriptome in a useful format, avoiding over-representation of some of the abundant, but low-complexity non-translated transcripts.

RESULTS

We developed a combination of self-subtraction and directional cloning procedures for this purpose. Libraries were prepared from partially degraded (hydrolyzed) total RNA from three different species. A restriction endonuclease site was added to the 3' end during first-strand synthesis using a directional random-priming technique. The abundant non-polyadenylated rRNA and tRNA sequences were largely removed by using self-subtraction to equalize the representation of the various RNA species. Sequencing random clones from the libraries showed that 87% of clones were in the forward orientation with respect to known or predicted transcripts. 70% matched identified or predicted translated RNAs in the sequence databases. Abundant mRNAs were less frequent in the self-subtracted libraries compared to a non-subtracted mRNA library. 3% of the sequences were from known or hypothesized ncRNA loci, including five matches to miRNA loci.

CONCLUSION

We describe a simple method for making high-quality, directional, random-primed, cDNA libraries from small amounts of degraded total RNA. This technique is advantageous in situations where a cDNA library with complete but equalized representation of transcribed sequences, whether polyadenylated or not, is desired.

Transcriptome analysis of zebrafish embryogenesis using microarrays

Zebrafish (Danio rerio) is a well-recognized model for the study of vertebrate developmental genetics, yet at the same time little is known about the transcriptional events that underlie zebrafish embryogenesis. Here we have employed microarray analysis to study the temporal activity of developmentally regulated genes during zebrafish embryogenesis. Transcriptome analysis at 12 different embryonic time points covering five different developmental stages (maternal, blastula, gastrula, segmentation, and pharyngula) revealed a highly dynamic transcriptional profile. Hierarchical clustering, stage-specific clustering, and algorithms to detect onset and peak of gene expression revealed clearly demarcated transcript clusters with maximum gene activity at distinct developmental stages as well as co-regulated expression of gene groups involved in dedicated functions such as organogenesis. Our study also revealed a previously unidentified cohort of genes that are transcribed prior to the mid-blastula transition, a time point earlier than when the zygotic genome was traditionally thought to become active. Here we provide, for the first time to our knowledge, a comprehensive list of developmentally regulated zebrafish genes and their expression profiles during embryogenesis, including novel information on the temporal expression of several thousand previously uncharacterized genes. The expression data generated from this study are accessible to all interested scientists from our institute resource database (http://giscompute.gis.a-star.edu.sg/~govind/zebrafish/data_download.html).

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.

ZebrafishIRX1bin the embryonic cardiac ventricle

The synchronous contraction of the vertebrate heart requires a conduction system. While coordinated contraction of the cardiac chambers is observed in zebrafish larvae, no histological evidence yet has been found for the existence of a cardiac conduction system in this tractable teleost. The homeodomain transcription factor gene IRX1 has been shown in the mouse embryo to be a marker of cells that give rise to the distinctive cardiac ventricular conduction system. Here, I demonstrate that zebrafish IRX1b is expressed in a restricted subset of ventricular myocytes within the embryonic zebrafish heart. IRX1b expression occurs as the electrical maturation of the heart is taking place, in a location analogous to the initial expression domain of mouse IRX1. The gene expression pattern of IRX1b is altered in silent heart genetic mutant embryos and in embryos treated with the endothelin receptor antagonist bosentan. Furthermore, injection of a morpholino oligonucleotide targeted to block IRX1b translation slows the heart rate.

Gene expression profile of zebrafish exposed to hypoxia during development

Understanding how vertebrates respond to hypoxia can have important clinical implications. Fish have evolved the ability to survive long exposure to low oxygen levels. However, little is known about the specific changes in gene expression that result from hypoxia. In this study we used a zebrafish cDNA microarray to examine the expression of >4,500 genes in zebrafish embryos exposed to 24 h of hypoxia during development. We tested the hypotheses that hypoxia changes gene expression profile of the zebrafish embryos and that these changes can be reverted by reexposure to a normoxic (20.8% O(2)) environment. Our data were consistent with both of these hypotheses: indicating that zebrafish embryos undergo adaptive changes in gene expression in response to hypoxia. Our study provides a striking genetic portrait of the zebrafish embryos' adaptive responses to hypoxic stress and demonstrates the utility of the microarray technology as a tool for analyzing complex developmental processes in the zebrafish.

Construction of a zebrafish cDNA microarray: gene expression profiling of the zebrafish during development

Vertebrate embryogenesis is a complex process controlled by a transcriptional hierarchy that coordinates the action of thousands of genes. To identify and analyze the expression patterns of these genes, we constructed a zebrafish cDNA microarray containing 4512 unique genes identified from zebrafish embryonic heart, adult hearts, and skeletal muscle cDNA libraries. We examined the patterns of gene expression during development in the zebrafish between five time points relative to 12h post-fertilization (hpf). Differentially expressed genes can be grouped into two categories, early genes that are expressed at 5hpf and genes expressed at 48/72/120hpf. Furthermore, we report the utilization of cDNA microarray technology to investigate the adaptive molecular responses of zebrafish to hypoxia during development. Our study provides the first utilization of cDNA microarray in the zebrafish and reveals dynamic changes in levels of gene expression in relation to development and survival of the zebrafish embryos under hypoxic stress.

Expressed sequence tag analysis of expression profiles of zebrafish testis and ovary

In the present study, two gonad cDNA libraries from zebrafish testes and ovaries were constructed and a total of 1025 expressed sequence tag (EST) clones were generated from the two libraries: 501 from the testis library and 524 from the ovary library. A total of 641 of the EST clones were identified to share significant sequence identity with known sequences in GenBank, representing at least 478 different zebrafish genes. In order to understand the molecular compositions of the two gonad organs, the expression profiles of the identified clones in these two gonad cDNA libraries were analyzed. Both gonad libraries have a higher portion of clones for nuclear proteins and a lower portion for proteins in translational machinery, cytoskeleton and mitochondria than our previously characterized whole-adult cDNA library. Most abundant cDNA clones in the two gonad libraries were identified and over 10% of ovary clones were found to encode egg membrane proteins (zona pellucida or ZP proteins). Furthermore, the testis library showed a more even distribution of cDNA clones with relatively fewer abundant clones that tend to contribute redundant clones in EST projects; thus, the testis library can supply more unique and novel cDNA sequences in a zebrafish EST project. Another aim of this study is to identify cDNA clones that can be used as molecular markers for the analysis of the gonad development in zebrafish. Eleven potential clones were selected to analyze their expression patterns by Northern blot hybridization. Most of them showed a specific or predominant expression in the expected testis or ovary tissue. At last, four of the clones were found, by section in situ hybridization, to be expressed specifically in the germ cells of the testis or ovary and thus they are suitable molecular markers for analyses of spermatogenesis and oogenesis.

Cardiac patterning and morphogenesis in zebrafish

Development of the embryonic vertebrate heart requires the precise coordination of pattern formation and cell movement. Taking advantage of the availability of zebrafish mutations that disrupt cardiogenesis, several groups have identified key regulators of specific aspects of cardiac patterning and morphogenesis. Several genes, including gata5, fgf8, bmp2b, one-eyed pinhead, and hand2, have been shown to be relevant to the patterning events that regulate myocardial differentiation. Studies of mutants with morphogenetic defects have indicated at least six genes that are essential for cardiac fusion and heart tube assembly, including casanova, bonnie and clyde, gata5, one-eyed pinhead, hand2, miles apart, and heart and soul. Furthermore, analysis of the jekyll gene has indicated its important role during the morphogenesis of the atrioventricular valve. Altogether, these data provide a substantial foundation for future investigations of cardiac patterning, cardiac morphogenesis, and the relationship between these processes.

Immune-type receptor genes in zebrafish share genetic and functional properties with genes encoded by the mammalian leukocyte receptor cluster

An extensive, highly diversified multigene family of novel immune-type receptor (nitr) genes has been defined in Danio rerio (zebrafish). The genes are predicted to encode type I transmembrane glycoproteins consisting of extracellular variable (V) and V-like C2 (V/C2) domains, a transmembrane region and a cytoplasmic tail. All of the genes examined encode immunoreceptor tyrosine-based inhibition motifs in the cytoplasmic tail. Radiation hybrid panel mapping and analysis of a deletion mutant line (b240) indicate that a minimum of approximately 40 nitr genes are contiguous in the genome and span approximately 0.6 Mb near the top of zebrafish linkage group 7. One flanking region of the nitr gene complex shares conserved synteny with a region of mouse chromosome 7, which shares conserved synteny with human 19q13.3-q13.4 that encodes the leukocyte receptor cluster. Antibody-induced crosslinking of Nitrs that have been introduced into a human natural killer cell line inhibits the phosphorylation of mitogen-activated protein kinase that is triggered by natural killer-sensitive tumor target cells. Nitrs likely represent intermediates in the evolution of the leukocyte receptor cluster.

Genomics and the pathophysiology of heart failure

Heart failure is not a single disease entity, but a syndrome with various causes, including hypertension, ischemic and congenital heart disease, cardiomyopathy, and myocarditis. Because of the multiple etiologies and secondary adaptations contributing to heart failure, the study of the cellular and molecular mechanisms underlying the development and progression of this syndrome has been rather challenging. Much has been learned about the remodeling processes in heart failure, which involve complex interactions among numerous mediators in signaling and regulatory pathways. The Human Genome Project and related projects have provided a preliminary database for a genome-wide analysis of complex polygenic disorders such as heart failure. With the aid of expressed sequence tag technology and microarray applications, both known and previously uncharacterized genes involved in the induction and regression of cardiac hypertrophy and its progression to heart failure can be analyzed simultaneously. Deciphering the complexity of sequence-structure-function relationships in heart failure is a goal for the future, and will require advances in structural biology, proteomics, and computational technology. In this review, we summarize the cellular and molecular aspects of heart failure, and how recent applications of genomic technologies have been successful in achieving a more complete portrait of gene expression in this pathologic state.