An EST resource for tilapia based on 17 normalized libraries and assembly of 116,899 sequence tags

Mapping QTL for Sex and Growth Traits in Salt-Tolerant Tilapia (Oreochromis spp. X O. mossambicus)

In aquaculture, growth and sex are economically important traits. To accelerate genetic improvement in increasing growth in salt-tolerant tilapia, we conducted QTL mapping for growth traits and sex with an F2 family, including 522 offspring and two parents. We used 144 polymorphic microsatellites evenly covering the genome of tilapia to genotype the family. QTL analyses were carried out using interval mapping for all individuals, males and females in the family, respectively. Using all individuals, three suggestive QTL for body weight, body length and body thickness respectively were detected in LG20, LG22 and LG12 and explained 2.4% to 3.1% of phenotypic variance (PV). When considering only males, five QTL for body weight were detected on five LGs, and explained 4.1 to 6.3% of PV. Using only females from the F2 family, three QTL for body weight were detected on LG1, LG6 and LG8, and explained 7.9-14.3% of PV. The QTL for body weight in males and females were located in different LGs, suggesting that in salt-tolerant tilapia, different set of genes 'switches' control the growth in males and females. QTL for sex were mapped on LG1 and LG22, indicating multigene sex determination in the salt-tolerant tilapia. This study provides new insights on the locations and effects of QTL for growth traits and sex, and sets the foundation for fine mapping for future marker-assisted selection for growth and sex in salt-tolerant tilapia aquaculture.

Signatures of selection in tilapia revealed by whole genome resequencing

Natural selection and selective breeding for genetic improvement have left detectable signatures within the genome of a species. Identification of selection signatures is important in evolutionary biology and for detecting genes that facilitate to accelerate genetic improvement. However, selection signatures, including artificial selection and natural selection, have only been identified at the whole genome level in several genetically improved fish species. Tilapia is one of the most important genetically improved fish species in the world. Using next-generation sequencing, we sequenced the genomes of 47 tilapia individuals. We identified a total of 1.43 million high-quality SNPs and found that the LD block sizes ranged from 10–100 kb in tilapia. We detected over a hundred putative selective sweep regions in each line of tilapia. Most selection signatures were located in non-coding regions of the tilapia genome. The Wnt signaling, gonadotropin-releasing hormone receptor and integrin signaling pathways were under positive selection in all improved tilapia lines. Our study provides a genome-wide map of genetic variation and selection footprints in tilapia, which could be important for genetic studies and accelerating genetic improvement of tilapia.

Living without DAT: Loss and compensation of the dopamine transporter gene in sauropsids (birds and reptiles)

The dopamine transporter (DAT) is a major regulator of synaptic dopamine (DA) availability. It plays key roles in motor control and motor learning, memory formation, and reward-seeking behavior, is a major target of cocaine and methamphetamines, and has been assumed to be conserved among vertebrates. We have found, however, that birds, crocodiles, and lizards lack the DAT gene. We also found that the unprecedented loss of this important gene is compensated for by the expression of the noradrenaline transporter (NAT) gene, and not the serotonin transporter genes, in dopaminergic cells, which explains the peculiar pharmacology of the DA reuptake activity previously noted in bird striatum. This unexpected pattern contrasts with that of ancestral vertebrates (e.g. fish) and mammals, where the NAT gene is selectively expressed in noradrenergic cells. DA circuits in birds/reptiles and mammals thus operate with an analogous reuptake mechanism exerted by different genes, bringing new insights into gene expression regulation in dopaminergic cells and the evolution of a key molecular player in reward and addiction pathways.

De novo sequencing and characterization of the Bradysia odoriphaga (Diptera: Sciaridae) larval transcriptome

The most serious pestilent threat to the Chinese chive, Allium tuberosum Rottle ex Spreng (Liliaceae) is the Bradysia odoriphaga Yang and Zhang. There is limited genetic research focused on B. odoriphaga, partially due to the lack of genomic resources. The advent of high-throughput sequencing technologies has enabled generation of genomic resources in a short time frame and at minimal costs. In this study, we performed, for the first time, de novo transcriptome sequencing of the B. odoriphaga. Here, 16,829 unigenes were assembled from the total reads, 12,024 of these unigenes were annotated in the NCBI NR protein database, and 9784 were annotated in the Swiss-Prot database. Of these annotated unigenes, 7903 and 5060 unigenes have been assigned to gene ontology categories and clusters of orthologous groups, respectively. Furthermore, 8647 unigenes were mapped to 257 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. We found that 408 unigenes were related to insecticide resistance and metabolism. In addition, 23,122 simple sequence repeats (SSRs) were identified in 11,009 unigenes, and 100 PCR primers of SSR loci were used to validate the assembly quality and polymorphisms. These results provide a good platform for further investigations into the insecticide resistance of B. odoriphaga. Finally, the SSRs identified in B. odoriphaga may be a useful genomic resource.

Casting a wider fish net on animal models in neuropsychiatric research

Neuropsychiatric disorders, such as schizophrenia, are associated with abnormal brain development. In this review, we discuss how studying dimensional components of these disorders, or endophenotypes, in a wider range of animal models will deepen our understanding of how interactions between biological and environmental factors alter the trajectory of neurodevelopment leading to aberrant behavior. In particular, we discuss some of the advantages of incorporating studies of brain and behavior using a range of teleost fish species into current neuropsychiatric research. From the perspective of comparative neurobiology, teleosts share a fundamental pattern of neurodevelopment and functional brain organization with other vertebrates, including humans. These shared features provide a basis for experimentally probing the mechanisms of disease-associated brain abnormalities. Moreover, incorporating information about how behaviors have been shaped by evolution will allow us to better understand the relevance of behavioral variation to determine their physiological underpinnings. We believe that exploiting the conservation in brain development across vertebrate species, and the rich diversity of fish behavior in lab and natural populations will lead to significant new insights and a holistic understanding of the neurobiological systems implicated in neuropsychiatric disorders.

Nutritional status and growth hormone regulate insulin-like growth factor binding protein (igfbp) transcripts in Mozambique tilapia

Growth in teleosts is controlled in large part by the activities of the growth hormone (Gh)/insulin-like growth factor (Igf) system. In this study, we initially identified igf-binding protein (bp)1b, -2b, -4, -5a and -6b transcripts in a tilapia EST library. In Mozambique tilapia (Oreochromis mossambicus), tissue expression profiling of igfbps revealed that igfbp1b and -2b had the highest levels of expression in liver while igfbp4, -5a and -6b were expressed at comparable levels in most other tissues. We compared changes in hepatic igfbp1b, -2b and -5a expression during catabolic conditions (28days of fasting) along with key components of the Gh/Igf system, including plasma Gh and Igf1 and hepatic gh receptor (ghr2), igf1 and igf2 expression. In parallel with elevated plasma Gh and decreased Igf1 levels, we found that hepatic igfbp1b increased substantially in fasted animals. We then tested whether systemic Gh could direct the expression of igfbps in liver. A single intraperitoneal injection of ovine Gh into hypophysectomized tilapia specifically stimulated liver igfbp2b expression along with plasma Igf1 and hepatic ghr2 levels. Our collective data suggest that hepatic endocrine signaling during fasting may involve post-translational regulation of plasma Igf1 via a shift towards the expression of igfbp1b. Thus, Igfbp1b may operate as a molecular switch to restrict Igf1 signaling in tilapia; furthermore, we provide new details regarding isoform-specific regulation of igfbp expression by Gh.

Transcriptome of the Deep-Sea Black Scabbardfish, Aphanopus carbo (Perciformes: Trichiuridae): Tissue-Specific Expression Patterns and Candidate Genes Associated to Depth Adaptation

Deep-sea fishes provide a unique opportunity to study the physiology and evolutionary adaptation to extreme environments. We carried out a high throughput sequencing analysis on a 454 GS-FLX titanium plate using unnormalized cDNA libraries from six tissues of A. carbo. Assemblage and annotations were performed by Newbler and InterPro/Pfam analyses, respectively. The assembly of 544,491 high quality reads provided 8,319 contigs, 55.6% of which retrieved blast hits against the NCBI nonredundant database or were annotated with ESTscan. Comparison of functional genes at both the protein sequences and protein stability levels, associated with adaptations to depth, revealed similarities between A. carbo and other bathypelagic fishes. A selection of putative genes was standardized to evaluate the correlation between number of contigs and their normalized expression, as determined by qPCR amplification. The screening of the libraries contributed to the identification of new EST simple-sequence repeats (SSRs) and to the design of primer pairs suitable for population genetic studies as well as for tagging and mapping of genes. The characterization of the deep-sea fish A. carbo first transcriptome is expected to provide abundant resources for genetic, evolutionary, and ecological studies of this species and the basis for further investigation of depth-related adaptation processes in fishes.

Identification of male-specific amh duplication, sexually differentially expressed genes and microRNAs at early embryonic development of Nile tilapia (Oreochromis niloticus)

Background

The probable influence of genes and the environment on sex determination in Nile tilapia suggests that it should be regarded as a complex trait. Detection of sex determination genes in tilapia has both scientific and commercial importance. The main objective was to detect genes and microRNAs that were differentially expressed by gender in early embryonic development.

Results

Artificial fertilization of Oreochromis niloticus XX females with either sex-reversed ΔXX males or genetically-modified YY ‘supermales’ resulted in all-female and all-male embryos, respectively. RNA of pools of all-female and all-male embryos at 2, 5 and 9 dpf were used as template for a custom Agilent eArray hybridization and next generation sequencing. Fifty-nine genes differentially expressed between genders were identified by a false discovery rate of p < 0.05. The most overexpressed genes were amh and tspan8 in males, and cr/20β-hsd, gpa33, rtn4ipl and zp3 in females (p < 1 × 10⁻⁹). Validation of gene expression using qPCR in embryos and gonads indicated copy number variation in tspan8, gpa33, cr/20β-hsd and amh. Sequencing of amh identified a male-specific duplication of this gene, denoted amhy, differing from the sequence of amh by a 233 bp deletion on exonVII, hence lacking the capability to encode the protein motif that binds to the transforming growth factor beta receptor (TGF-β domain). amh and amhy segregated in the mapping family in full concordance with SD-linked marker on LG23 signifying the QTL for SD. We discovered 831 microRNAs in tilapia embryos of which nine had sexually dimorphic expression patterns by a false discovery rate of p < 0.05. An up-regulated microRNA in males, pma-mir-4585, was characterized with all six predicted target genes including cr/20β-hsd, down-regulated in males.

Conclusions

This study reports the first discovery of sexually differentially expressed genes and microRNAs at a very early stage of tilapia embryonic development, i.e. from 2 dpf. Genes with sexually differential expression patterns are enriched for copy number variation. A novel male-specific duplication of amh, denoted amhy, lacking the TGF-β domain was identified and mapped to the QTL region on LG23 for SD, thus indicating its potential role in SD.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-774) contains supplementary material, which is available to authorized users.

Allelic variant in the anti-Müllerian hormone gene leads to autosomal and temperature-dependent sex reversal in a selected Nile tilapia line

Owing to the demand for sustainable sex-control protocols in aquaculture, research in tilapia sex determination is gaining momentum. The mutual influence of environmental and genetic factors hampers disentangling the complex sex determination mechanism in Nile tilapia (Oreochromis niloticus). Previous linkage analyses have demonstrated quantitative trait loci for the phenotypic sex on linkage groups 1, 3, and 23. Quantitative trait loci for temperature-dependent sex reversal similarly reside on linkage group 23. The anti-Müllerian hormone gene (amh), located in this genomic region, is important for sexual fate in higher vertebrates, and shows sexually dimorphic expression in Nile tilapia. Therefore this study aimed at detecting allelic variants and marker-sex associations in the amh gene. Sequencing identified six allelic variants. A significant effect on the phenotypic sex for SNP ss831884014 (p<0.0017) was found by stepwise logistic regression. The remaining variants were not significantly associated. Functional annotation of SNP ss831884014 revealed a non-synonymous amino acid substitution in the amh protein. Consequently, a fluorescence resonance energy transfer (FRET) based genotyping assay was developed and validated with a representative sample of fish. A logistic linear model confirmed a highly significant effect of the treatment and genotype on the phenotypic sex, but not for the interaction term (treatment: p<0.0001; genotype: p<0.0025). An additive genetic model proved a linear allele substitution effect of 12% in individuals from controls and groups treated at high temperature, respectively. Moreover, the effect of the genotype on the male proportion was significantly higher in groups treated at high temperature, giving 31% more males on average of the three genotypes. In addition, the groups treated at high temperature showed a positive dominance deviation (+11.4% males). In summary, marker-assisted selection for amh variant ss831884014 seems to be highly beneficial to increase the male proportion in Nile tilapia, especially when applying temperature-induced sex reversal.

Evolution of genomic structural variation and genomic architecture in the adaptive radiations of African cichlid fishes

African cichlid fishes are an ideal system for studying explosive rates of speciation and the origin of diversity in adaptive radiation. Within the last few million years, more than 2000 species have evolved in the Great Lakes of East Africa, the largest adaptive radiation in vertebrates. These young species show spectacular diversity in their coloration, morphology and behavior. However, little is known about the genomic basis of this astonishing diversity. Recently, five African cichlid genomes were sequenced, including that of the Nile Tilapia (Oreochromis niloticus), a basal and only relatively moderately diversified lineage, and the genomes of four representative endemic species of the adaptive radiations, Neolamprologus brichardi, Astatotilapia burtoni, Metriaclima zebra, and Pundamila nyererei. Using the Tilapia genome as a reference genome, we generated a high-resolution genomic variation map, consisting of single nucleotide polymorphisms (SNPs), short insertions and deletions (indels), inversions and deletions. In total, around 18.8, 17.7, 17.0, and 17.0 million SNPs, 2.3, 2.2, 1.4, and 1.9 million indels, 262, 306, 162, and 154 inversions, and 3509, 2705, 2710, and 2634 deletions were inferred to have evolved in N. brichardi, A. burtoni, P. nyererei, and M. zebra, respectively. Many of these variations affected the annotated gene regions in the genome. Different patterns of genetic variation were detected during the adaptive radiation of African cichlid fishes. For SNPs, the highest rate of evolution was detected in the common ancestor of N. brichardi, A. burtoni, P. nyererei, and M. zebra. However, for the evolution of inversions and deletions, we found that the rates at the terminal taxa are substantially higher than the rates at the ancestral lineages. The high-resolution map provides an ideal opportunity to understand the genomic bases of the adaptive radiation of African cichlid fishes.

Whole genome comparative analysis of channel catfish (Ictalurus punctatus) with four model fish species

BACKGROUND

Comparative mapping is a powerful tool to study evolution of genomes. It allows transfer of genome information from the well-studied model species to non-model species. Catfish is an economically important aquaculture species in United States. A large amount of genome resources have been developed from catfish including genetic linkage maps, physical maps, BAC end sequences (BES), integrated linkage and physical maps using BES-derived markers, physical map contig-specific sequences, and draft genome sequences. Application of such genome resources should allow comparative analysis at the genome scale with several other model fish species.

RESULTS

In this study, we conducted whole genome comparative analysis between channel catfish and four model fish species with fully sequenced genomes, zebrafish, medaka, stickleback and Tetraodon. A total of 517 Mb draft genome sequences of catfish were anchored to its genetic linkage map, which accounted for 62% of the total draft genome sequences. Based on the location of homologous genes, homologous chromosomes were determined among catfish and the four model fish species. A large number of conserved syntenic blocks were identified. Analysis of the syntenic relationships between catfish and the four model fishes supported that the catfish genome is most similar to the genome of zebrafish.

CONCLUSION

The organization of the catfish genome is similar to that of the four teleost species, zebrafish, medaka, stickleback, and Tetraodon such that homologous chromosomes can be identified. Within each chromosome, extended syntenic blocks were evident, but the conserved syntenies at the chromosome level involve extensive inter-chromosomal and intra-chromosomal rearrangements. This whole genome comparative map should facilitate the whole genome assembly and annotation in catfish, and will be useful for genomic studies of various other fish species.

Generation and analysis of a large-scale expressed sequence Tag database from a full-length enriched cDNA library of developing leaves of Gossypium hirsutum L

Background

Cotton (Gossypium hirsutum L.) is one of the world’s most economically-important crops. However, its entire genome has not been sequenced, and limited resources are available in GenBank for understanding the molecular mechanisms underlying leaf development and senescence.

Methodology/Principal Findings

In this study, 9,874 high-quality ESTs were generated from a normalized, full-length cDNA library derived from pooled RNA isolated from throughout leaf development during the plant blooming stage. After clustering and assembly of these ESTs, 5,191 unique sequences, representative 1,652 contigs and 3,539 singletons, were obtained. The average unique sequence length was 682 bp. Annotation of these unique sequences revealed that 84.4% showed significant homology to sequences in the NCBI non-redundant protein database, and 57.3% had significant hits to known proteins in the Swiss-Prot database. Comparative analysis indicated that our library added 2,400 ESTs and 991 unique sequences to those known for cotton. The unigenes were functionally characterized by gene ontology annotation. We identified 1,339 and 200 unigenes as potential leaf senescence-related genes and transcription factors, respectively. Moreover, nine genes related to leaf senescence and eleven MYB transcription factors were randomly selected for quantitative real-time PCR (qRT-PCR), which revealed that these genes were regulated differentially during senescence. The qRT-PCR for three GhYLSs revealed that these genes express express preferentially in senescent leaves.

Conclusions/Significance

These EST resources will provide valuable sequence information for gene expression profiling analyses and functional genomics studies to elucidate their roles, as well as for studying the mechanisms of leaf development and senescence in cotton and discovering candidate genes related to important agronomic traits of cotton. These data will also facilitate future whole-genome sequence assembly and annotation in G. hirsutum and comparative genomics among Gossypium species.

Cloning and molecular characterization of complement component 1 inhibitor (C1INH) and complement component 8β (C8β) in Nile tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus), one of the most important groups of food fishes in the world, has frequently suffered from serious challenge from pathogens in recent years. Immune responses of Nile tilapia should be understood to protect the aquaculture industry of this fish. The complement system has an important function in recognizing bacteria, opsonizing these pathogens by phagocytes, or killing them by direct lysis. In this study, two Nile tilapia complement component genes, complement component 1 inhibitor (C1INH) and complement component 8β subunit (C8β), were cloned and their expression characteristics were analyzed. C1INH cDNA was found containing a 1791 bp open reading frame (ORF) encoding a putative protein with 597 amino acids, a 101 bp 5'-untranslated region (UTR) and a 236 bp 3'-UTR. The predicted protein structure for this gene consisted of two Ig-like domains and glycosyl hydrolase family-9 active site signature 2. The C8β cDNA consisted of a 1761 bp ORF encoding 587 amino acids, a 15 bp 5'-UTR and a 170 bp 3'-UTR. The predicted protein of C8β contained three motifs, thrombospondin type-1 repeat, membrane attack complex/perforin domain, and LDL-receptor class A. Expression analysis revealed that these two complement genes were highly expressed in the liver, however, were weakly expressed in the gill, heart, brain, kidney, intestine, spleen and dorsal muscle tissues. The present study provided insights into the complement system and immune functions of Nile tilapia.

Comparative genomic analysis of catfish linkage group 8 reveals two homologous chromosomes in zebrafish and other teleosts with extensive inter-chromosomal rearrangements

Background

Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study.

Results

With existing genome resources, 287 unique genes were identified in LG8. Comparative genome analysis indicated that most of these 287 genes on catfish LG8 are located on two homologous chromosomes of zebrafish, medaka, stickleback, and three chromosomes of green-spotted pufferfish. Large numbers of conserved syntenies were identified. Detailed analysis of the conserved syntenies in relation to chromosome level similarities revealed extensive inter-chromosomal and intra-chromosomal rearrangements during evolution. Of the 287 genes, 35 genes were found to be duplicated in the catfish genome, with the vast majority of the duplications being interchromosomal.

Conclusions

Comparative genome analysis is a powerful tool even in the absence of a well-assembled whole genome sequence. In spite of sequence stacking due to low resolution of the linkage and physical maps, conserved syntenies can be identified although the exact gene order and orientation are unknown at present. Through chromosome-level comparative analysis, homologous chromosomes among teleosts can be identified. Syntenic analysis should facilitate annotation of the catfish genome, which in turn, should facilitate functional inference of genes based on their orthology.

A novel beta-defensin antimicrobial peptide in Atlantic cod with stimulatory effect on phagocytic activity

A novel defensin antimicrobial peptide gene was identified in Atlantic cod, Gadus morhua. This three exon/two intron defensin gene codes for a peptide precursor consisting of two domains: a signal peptide of 26 amino acids and a mature peptide of 40 residues. The mature cod defensin has six conserved cysteine residues that form 1–5, 2–4 and 3–6 disulphide bridges. This pattern is typical of beta-defensins and this gene was therefore named cod beta-defensin (defb). The tertiary structure of Defb exhibits an α/β fold with one α helix and β₁β₂β₃ sheets_. RT-PCR analysis indicated that defb transcripts were present mainly in the swim bladder and peritoneum wall but could also be detected at moderate to low levels in skin, head- and excretory kidneys. In situ hybridisation revealed that defb was specifically expressed by cells located in the swim bladder submucosa and the oocytes. During embryonic development, defb gene transcripts were detectable from the golden eye stage onwards and their expression was restricted to the swim bladder and retina. Defb was differentially expressed in several tissues following antigenic challenge with Vibrio anguillarum, being up-regulated up to 25-fold in head kidney. Recombinant Defb displayed antibacterial activity, with a minimal inhibitory concentration of 0.4–0.8 µM and 25–50 µM against the Gram-(+) bacteria Planococcus citreus and Micrococcus luteus, respectively. In addition, Defb stimulated phagocytic activity of cod head kidney leucocytes invitro. These findings imply that beta-defensins may play an important role in the innate immune response of Atlantic cod.

Deep sequencing for de novo construction of a marine fish (Sparus aurata) transcriptome database with a large coverage of protein-coding transcripts

Background

The gilthead sea bream (Sparus aurata) is the main fish species cultured in the Mediterranean area and constitutes an interesting model of research. Nevertheless, transcriptomic and genomic data are still scarce for this highly valuable species. A transcriptome database was constructed by de novo assembly of gilthead sea bream sequences derived from public repositories of mRNA and collections of expressed sequence tags together with new high-quality reads from five cDNA 454 normalized libraries of skeletal muscle (1), intestine (1), head kidney (2) and blood (1).

Results

Sequencing of the new 454 normalized libraries produced 2,945,914 high-quality reads and the de novo global assembly yielded 125,263 unique sequences with an average length of 727 nt. Blast analysis directed to protein and nucleotide databases annotated 63,880 sequences encoding for 21,384 gene descriptions, that were curated for redundancies and frameshifting at the homopolymer regions of open reading frames, and hosted at http://www.nutrigroup-iats.org/seabreamdb. Among the annotated gene descriptions, 16,177 were mapped in the Ingenuity Pathway Analysis (IPA) database, and 10,899 were eligible for functional analysis with a representation in 341 out of 372 IPA canonical pathways. The high representation of randomly selected stickleback transcripts by Blast search in the nucleotide gilthead sea bream database evidenced its high coverage of protein-coding transcripts.

Conclusions

The newly assembled gilthead sea bream transcriptome represents a progress in genomic resources for this species, as it probably contains more than 75% of actively transcribed genes, constituting a valuable tool to assist studies on functional genomics and future genome projects.

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

BACKGROUND

The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates.

RESULTS

The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates.

CONCLUSIONS

Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.

Transcriptomic characterization of the larval stage in gilthead seabream (Sparus aurata) by 454 pyrosequencing

Gilthead seabream (Sparus aurata) is a teleost belonging to the family Sparidae with a high economical relevance in the Mediterranean countries. Although genomic tools have been developed in this species in order to investigate its physiology at the molecular level and consequently its culture, genomic information on post-embryonic development is still scarce. In this study, we have investigated the transcriptome of a marine teleost during the larval stage (from hatching to 60 days after hatching) by the use of 454 pyrosequencing technology. We obtained a total of 68,289 assembled contigs, representing putative transcripts, belonging to 54,606 different clusters. Comparison against all S. aurata expressed sequenced tags (ESTs) from the NCBI database revealed that up to 34,722 contigs, belonging to about 61% of gene clusters, are sequences previously not described. Contigs were annotated through an iterative Blast pipeline by comparison against databases such as NCBI RefSeq from Danio rerio, SwissProt or NCBI teleost ESTs. Our results indicate that we have enriched the number of annotated sequences for this species by more than 50% compared with previously existing databases for the gilthead seabream. Gene Ontology analysis of these novel sequences revealed that there is a statistically significant number of transcripts with key roles in larval development, differentiation, morphology, and growth. Finally, all information has been made available online through user-friendly interfaces such as GBrowse and a Blast server with a graphical frontend.

A novel primer set for multilocus phylogenetic inference in East African cichlid fishes

The cichlid fishes in the East African Great Lakes are a prime model system for the study of adaptive radiation. Therefore, the availability of an elaborate phylogenetic framework is an important prerequisite. Previous phylogenetic hypotheses on East African cichlids are mainly based on mitochondrial and/or fragment-based markers, and, to date, no taxon-rich phylogeny exists that is based on multilocus DNA sequence data. Here, we present the design of an extensive new primer set (24 nuclear makers) for East African cichlids that will be used for multilocus phylogenetic analyses in the future. The primers are designed to work for both Sanger sequencing and next-generation sequencing with the 454 technology. As a proof of principle, we validate these primers in a phylogenetically representative set of 16 cichlid species from Lake Tanganyika and main river systems in the area and provide a basic evaluation of the markers with respect to marker length and diversity indices.

Transcriptome analysis of the Octopus vulgaris central nervous system

Background

Cephalopoda are a class of Mollusca species found in all the world's oceans. They are an important model organism in neurobiology. Unfortunately, the lack of neuronal molecular sequences, such as ESTs, transcriptomic or genomic information, has limited the development of molecular neurobiology research in this unique model organism.

Results

With high-throughput Illumina Solexa sequencing technology, we have generated 59,859 high quality sequences from 12,918,391 paired-end reads. Using BLASTx/BLASTn, 12,227 contigs have blast hits in the Swissprot, NR protein database and NT nucleotide database with E-value cutoff 1e⁻⁵. The comparison between the Octopus vulgaris central nervous system (CNS) library and the Aplysia californica/Lymnaea stagnalis CNS ESTs library yielded 5.93%/13.45% of O. vulgaris sequences with significant matches (1e⁻⁵) using BLASTn/tBLASTx. Meanwhile the hit percentage of the recently published Schistocerca gregaria, Tilapia or Hirudo medicinalis CNS library to the O. vulgaris CNS library is 21.03%–46.19%. We constructed the Phylogenetic tree using two genes related to CNS function, Synaptotagmin-7 and Synaptophysin. Lastly, we demonstrated that O. vulgaris may have a vertebrate-like Blood-Brain Barrier based on bioinformatic analysis.

Conclusion

This study provides a mass of molecular information that will contribute to further molecular biology research on O. vulgaris. In our presentation of the first CNS transcriptome analysis of O. vulgaris, we hope to accelerate the study of functional molecular neurobiology and comparative evolutionary biology.

Transcriptome characterization via 454 pyrosequencing of the annelid Pristina leidyi, an emerging model for studying the evolution of regeneration

Background

The naid annelids contain a number of species that vary in their ability to regenerate lost body parts, making them excellent candidates for evolution of regeneration studies. However, scant sequence data exists to facilitate such studies. We constructed a cDNA library from the naid Pristina leidyi, a species that is highly regenerative and also reproduces asexually by fission, using material from a range of regeneration and fission stages for our library. We then sequenced the transcriptome of P. leidyi using 454 technology.

Results

454 sequencing produced 1,550,174 reads with an average read length of 376 nucleotides. Assembly of 454 sequence reads resulted in 64,522 isogroups and 46,679 singletons for a total of 111,201 unigenes in this transcriptome. We estimate that over 95% of the transcripts in our library are present in our transcriptome. 17.7% of isogroups had significant BLAST hits to the UniProt database and these include putative homologs of a number of genes relevant to regeneration research. Although many sequences are incomplete, the mean sequence length of transcripts (isotigs) is 707 nucleotides. Thus, many sequences are large enough to be immediately useful for downstream applications such as gene expression analyses. Using in situ hybridization, we show that two Wnt/β-catenin pathway genes (homologs of frizzled and β-catenin) present in our transcriptome are expressed in the regeneration blastema of P. leidyi, demonstrating the usefulness of this resource for regeneration research.

Conclusions

454 sequencing is a rapid and efficient approach for identifying large numbers of genes in an organism that lacks a sequenced genome. This transcriptome dataset will be a valuable resource for molecular analyses of regeneration in P. leidyi and will serve as a starting point for comparisons to non-regenerating naids. It also contributes significantly to the still limited genomic resources available for annelids and lophotrochozoans more generally.

A high-resolution map of the Nile tilapia genome: a resource for studying cichlids and other percomorphs

Background

The Nile tilapia (Oreochromis niloticus) is the second most farmed fish species worldwide. It is also an important model for studies of fish physiology, particularly because of its broad tolerance to an array of environments. It is a good model to study evolutionary mechanisms in vertebrates, because of its close relationship to haplochromine cichlids, which have undergone rapid speciation in East Africa. The existing genomic resources for Nile tilapia include a genetic map, BAC end sequences and ESTs, but comparative genome analysis and maps of quantitative trait loci (QTL) are still limited.

Results

We have constructed a high-resolution radiation hybrid (RH) panel for the Nile tilapia and genotyped 1358 markers consisting of 850 genes, 82 markers corresponding to BAC end sequences, 154 microsatellites and 272 single nucleotide polymorphisms (SNPs). From these, 1296 markers could be associated in 81 RH groups, while 62 were not linked. The total size of the RH map is 34,084 cR₃₅₀₀ and 937,310 kb. It covers 88% of the entire genome with an estimated inter-marker distance of 742 Kb. Mapping of microsatellites enabled integration to the genetic map. We have merged LG8 and LG24 into a single linkage group, and confirmed that LG16-LG21 are also merged. The orientation and association of RH groups to each chromosome and LG was confirmed by chromosomal in situ hybridizations (FISH) of 55 BACs. Fifty RH groups were localized on the 22 chromosomes while 31 remained small orphan groups. Synteny relationships were determined between Nile tilapia, stickleback, medaka and pufferfish.

Conclusion

The RH map and associated FISH map provide a valuable gene-ordered resource for gene mapping and QTL studies. All genetic linkage groups with their corresponding RH groups now have a corresponding chromosome which can be identified in the karyotype. Placement of conserved segments indicated that multiple inter-chromosomal rearrangements have occurred between Nile tilapia and the other model fishes. These maps represent a valuable resource for organizing the forthcoming genome sequence of Nile tilapia, and provide a foundation for evolutionary studies of East African cichlid fishes.

Gender and sexual behavior modulate the composition of serum lipocalins in Nile tilapia (Oreochromis niloticus)

In tilapia species, plasma lipoproteins with high electrophoretic mobility function in intra- and intergender communication. Blood samples taken at onset and peak of daily sexual activity from dominant and subordinate Oreochromis niloticus males and females were fractionated by native gel electrophoresis and the fast-migrating proteins were subjected to mass spectrometry. Mining the sequence data of the Cichlid Genome Consortium, we identified 11 proteins from the lipocalin super-family and characterized their genes' structures. Phylogenetic and structural analyses subdivided these genes into two classes: (I) 3-coding-exon apolipoproteins and (II) more complex 6-coding-exon sulfide-bond-containing lipocalins. Five apolipoproteins and PTGDSL1, TBTBP, and MSP proteins were modulated by gender and sexual behavior. PTGDSL1 protein was only observed in the plasma serum of dominant males. However, the cysteine residue in the position that is crucial for synthetase activity in mammalian prostaglandin D synthetases was not conserved in PTGDSL1 or PTGDSL2 proteins. In line with previous reports suggesting their involvement in male functions as pheromone transporters, TBTBP and MSP proteins were not detected in females at the onset of daily activity. Their increasing amount in males was concordant with the increase in apolipoproteins AFP4L, APOA4a, APOA4b, APO14kD and APOC2, which were detected exclusively in dominant males, indicating a possible role in mobilization of the energy required to maintain their social hierarchy.

SNP marker detection and genotyping in tilapia

We have generated a unique resource consisting of nearly 175 000 short contig sequences and 3569 SNP markers from the widely cultured GIFT (Genetically Improved Farmed Tilapia) strain of Nile tilapia (Oreochromis niloticus). In total, 384 SNPs were selected to monitor the wider applicability of the SNPs by genotyping tilapia individuals from different strains and different geographical locations. In all strains and species tested (O. niloticus, O. aureus and O. mossambicus), the genotyping assay was working for a similar number of SNPs (288-305 SNPs). The actual number of polymorphic SNPs was, as expected, highest for individuals from the GIFT population (255 SNPs). In the individuals from an Egyptian strain and in individuals caught in the wild in the basin of the river Volta, 197 and 163 SNPs were polymorphic, respectively. A pairwise calculation of Nei's genetic distance allowed the discrimination of the individual strains and species based on the genotypes determined with the SNP set. We expect that this set will be widely applicable for use in tilapia aquaculture, e.g. for pedigree reconstruction. In addition, this set is currently used for assaying the genetic diversity of native Nile tilapia in areas where tilapia is, or will be, introduced in aquaculture projects. This allows the tracing of escapees from aquaculture and the monitoring of effects of introgression and hybridization.

Linkage and Physical Mapping of Sex Region on LG23 of Nile Tilapia (Oreochromis niloticus)

Evidence supports that sex determination (SD) in tilapia is controlled by major genetic factors that may interact with minor genetic as well as environmental factors, thus implying that SD should be analyzed as a quantitative trait. Quantitative trait loci (QTL) for SD in Oreochromis niloticus were previously detected on linkage groups (LG) 1 and 23. Twenty-one short single repeats (SSR) of >12 TGs and one single nucleotide polymorphism were identified using the unpublished tilapia genome sequence on LG23. All markers showed two segregating alleles in a mapping family that was obtained by a cross between O. niloticus male (XY) and sex-reversed female (ΔXY) yielding 29 females (XX) and 61 males (XY and YY). Interval mapping analysis mapped the QTL peak between SSR markers ARO172 and ARO177 with a maximum F value of 78.7 (P < 7.6 × 10(-14)). Twelve adjacent markers found in this region were homozygous in females and either homozygous for the alternative allele or heterozygous in males. This segment was defined as the sex region (SR). The SR encompasses 1.5 Mbp on a single tilapia scaffold (no. 101) harboring 51 annotated genes. Among 10 candidate genes for SD that were tested for gene expression, anti-Müllerian hormone (Amh), which is located in the center of the SR, showed the highest overexpression in male vs. female embryos at 3 to 7 days postfertilization.

Generation of ESTs for flowering gene discovery and SSR marker development in upland cotton

Background

Upland cotton, Gossypium hirsutum L., is one of the world's most important economic crops. In the absence of the entire genomic sequence, a large number of expressed sequence tag (EST) resources of upland cotton have been generated and used in several studies. However, information about the flower development of this species is rare.

Methodology/Principal Findings

To clarify the molecular mechanism of flower development in upland cotton, 22,915 high-quality ESTs were generated and assembled into 14,373 unique sequences consisting of 4,563 contigs and 9,810 singletons from a normalized and full-length cDNA library constructed from pooled RNA isolated from shoot apexes, squares, and flowers. Comparative analysis indicated that 5,352 unique sequences had no high-degree matches to the cotton public database. Functional annotation showed that several upland cotton homologs with flowering-related genes were identified in our library. The majority of these genes were specifically expressed in flowering-related tissues. Three GhSEP (G. hirsutum L. SEPALLATA) genes determining floral organ development were cloned, and quantitative real-time PCR (qRT-PCR) revealed that these genes were expressed preferentially in squares or flowers. Furthermore, 670 new putative microsatellites with flanking sequences sufficient for primer design were identified from the 645 unigenes. Twenty-five EST–simple sequence repeats were randomly selected for validation and transferability testing in 17 Gossypium species. Of these, 23 were identified as true-to-type simple sequence repeat loci and were highly transferable among Gossypium species.

Conclusions/Significance

A high-quality, normalized, full-length cDNA library with a total of 14,373 unique ESTs was generated to provide sequence information for gene discovery and marker development related to upland cotton flower development. These EST resources form a valuable foundation for gene expression profiling analysis, functional analysis of newly discovered genes, genetic linkage, and quantitative trait loci analysis.

Positive Darwinian selection drives the evolution of the morphology-related gene, EPCAM, in particularly species-rich lineages of African cichlid fishes

The study of genetic evolution within the context of adaptive radiations offers insights to genes and selection pressures that result in rapid morphological change. Cichlid fishes are very species-rich and variable in coloration, behavior, and morphology, and so provide a classical model system for studying the genetics of adaptive radiation. In this study, we researched the evolution of the epithelial cell adhesion molecule (EPCAM), a candidate gene for the adaptive evolution of morphology broadly, and skin development specifically, in fishes. We compared EPCAM gene sequences from a rapidly speciating African cichlid lineage (the haplochromines), a species-poor African lineage (Nile tilapia Oreochromis niloticus), and a very young adaptive radiation in the Neotropics (sympatric crater lake Midas cichlids, Amphilophus sp.). Our results, based on a hierarchy of evolutionary analyses of nucleotide substitution, demonstrate that there are different selection pressures on the EPCAM gene among the cichlid lineages. Several waves of positive natural selection were identified not only on the terminal branches, but also on ancestral branches. Interestingly, significant positive or directional selection was found in the haplochromine cichlids only but not the comparatively species-poor tilapia lineage. We hypothesize that the strong signal of selection in the ancestral African cichlid lineage coincided with the transition from riverine to lacustrine habitat. The two neotropical species for which we collected new sequence data were invariant in the EPCAM locus. Our results suggest that functional changes promoted by positive Darwinian selection are widespread in the EPCAM gene during African cichlid evolution.

Comparative transcriptomics of Eastern African cichlid fishes shows signs of positive selection and a large contribution of untranslated regions to genetic diversity

The hundreds of endemic species of cichlid fishes in the East African Great Lakes Tanganyika, Malawi, and Victoria are a prime model system in evolutionary biology. With five genomes currently being sequenced, eastern African cichlids also represent a forthcoming genomic model for evolutionary studies of genotype-to-phenotype processes in adaptive radiations. Here we report the functional annotation and comparative analyses of transcriptome data sets for two eastern African cichlid species, Astatotilapia burtoni and Ophthalmotilapia ventralis, representatives of the modern haplochromines and ectodines, respectively. Nearly 647,000 expressed sequence tags were assembled in more than 46,000 contigs for each species using the 454 sequencing technology, largely expanding the current sequence data set publicly available for these cichlids. Total predicted coverage of their proteome diversity is approximately 50% for both species. Comparative qualitative and quantitative analyses show very similar transcriptome data for the two species in terms of both functional annotation and relative abundance of gene ontology terms expressed. Average genetic distance between species is 1.75% when all transcript types are considered including nonannotated sequences, 1.33% for annotated sequences only including untranslated regions, and decreases to nearly half, 0.95%, for coding sequences only, suggesting a large contribution of noncoding regions to their genetic diversity. Comparative analyses across the two species, tilapia and the outgroup medaka based on an overlapping data set of 1,216 genes (∼526 kb) demonstrate cichlid-specific signature of disruptive selection and provide a set of candidate genes that are putatively under positive selection. Overall, these data sets offer the genetic platform for future comparative analyses in light of the upcoming genomes for this taxonomic group.

Evolution of microRNAs and the diversification of species

MicroRNAs (miRNAs) are ancient, short noncoding RNA molecules that regulate the transcriptome through post-transcriptional mechanisms. miRNA riboregulation is involved in a diverse range of biological processes, and misregulation is implicated in disease. It is generally thought that miRNAs function to canalize cellular outputs, for instance as “fail-safe” repressors of gene misexpression. Genomic surveys in humans have revealed reduced genetic polymorphism and the signature of negative selection for both miRNAs themselves and the target sequences to which they are predicted to bind. We investigated the evolution of miRNAs and their binding sites across cichlid fishes from Lake Malawi (East Africa), where hundreds of diverse species have evolved in the last million years. Using low-coverage genome sequence data, we identified 100 cichlid miRNA genes with mature regions that are highly conserved in other animal species. We computationally predicted target sites on the 3′-untranslated regions (3′-UTRs) of cichlid genes to which miRNAs may bind and found that these sites possessed elevated single nucleotide polymorphism (SNP) densities. Furthermore, polymorphic sites in predicted miRNA targets showed higher minor allele frequencies on average and greater genetic differentiation between Malawi lineages when compared with a neutral expectation and nontarget 3′-UTR SNPs. Our data suggest that divergent selection on miRNA riboregulation may have contributed to the diversification of cichlid species and may similarly play a role in rapid phenotypic evolution of other natural systems.

Comparative physical maps derived from BAC end sequences of tilapia (Oreochromis niloticus)

Background

The Nile tilapia is the second most important fish in aquaculture. It is an excellent laboratory model, and is closely related to the African lake cichlids famous for their rapid rates of speciation. A suite of genomic resources has been developed for this species, including genetic maps and ESTs. Here we analyze BAC end-sequences to develop comparative physical maps, and estimate the number of genome rearrangements, between tilapia and other model fish species.

Results

We obtained sequence from one or both ends of 106,259 tilapia BACs. BLAST analysis against the genome assemblies of stickleback, medaka and pufferfish allowed identification of homologies for approximately 25,000 BACs for each species. We calculate that rearrangement breakpoints between tilapia and these species occur about every 3 Mb across the genome. Analysis of 35,000 clones previously assembled into contigs by restriction fingerprints allowed identification of longer-range syntenies.

Conclusions

Our data suggest that chromosomal evolution in recent teleosts is dominated by alternate loss of gene duplicates, and by intra-chromosomal rearrangements (~one per million years). These physical maps are a useful resource for comparative positional cloning of traits in cichlid fishes. The paired BAC end sequences from these clones will be an important resource for scaffolding forthcoming shotgun sequence assemblies of the tilapia genome.