A first generation BAC-based physical map of the rainbow trout genome

Fine scale diversity in the lava: genetic and phenotypic diversity in small populations of Arctic charr Salvelinus alpinus

BACKGROUND

A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure.

RESULTS

These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect.

CONCLUSIONS

Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable.

Genomic architecture and repertoire of the rainbow trout immunoglobulin light chain genes

The genomic loci encoding the four immunoglobulin light chains (IgL1, IgL2, IgL3, and IgL4) in the Swanson trout genome assembly were annotated in order to provide a measurement of the potential IgL repertoire. IgL1 and IgL3 gene segments are co-localized on chromosomes 21, 18, 15, and 7 while IgL2 and IgL4 were found on chromosomes 13 and 17, respectively. In total, 48 constant (C_L), 87 variable (V_L), and 59 joining (J_L) productive genes are described. Pairwise alignment of the V_L segments revealed that they belong to nine different families, three of which (kappa IV, V, and VI) are described for the first time in this study. V_L and C_L sequences on chromosome 15 and 21 and those on chromosomes 7 and 18 clustered together in phylogenetic analysis. PCR was used to examine IgL C_L and V_L genes in 9 lines of rainbow trout. IgL4 in the Hot Creek and Golden trout lines was missing 42 nucleotides resulting in a loss of 14 amino acids. The sigma IV variable family was completely absent from the Swanson, Arlee, Hot Creek, and wild type lines and silenced in the Skamania line with the addition of 176 bp mini-satellite insert. Similarly, the Whale Rock, Arlee, and wild type lines were all found to encode two sigma II products, a functional 252 bp product and a larger 425 bp product that contained a 172 bp insert. Results from this study indicate that there are genomic differences in IgL repertoire between different lines of trout that could affect humoral immune responses post vaccination and during disease.

Genome and Phylogenetic Analysis of Genes Involved in the Immune System of Solea senegalensis - Potential Applications in Aquaculture

Global aquaculture production continues to increase rapidly. One of the most important species of marine fish currently cultivated in Southern Europe is Solea senegalensis, reaching more than 300 Tn in 2017. In the present work, 14 Bacterial Artificial Chromosome (BAC) clones containing candidate genes involved in the immune system (b2m, il10, tlr3, tap1, tnfα, tlr8, trim25, lysg, irf5, hmgb2, calr, trim16, and mx), were examined and compared with other species using multicolor Fluorescence in situ Hybridization (mFISH), massive sequencing and bioinformatic analysis to determine the genomic surroundings and syntenic chromosomal conservation of the genomic region contained in each BAC clone. The mFISH showed that the groups of genes hmgb2-trim25-irf5-b2m; tlr3-lysg; tnfα-tap1, and il10-mx-trim16 were co-localized on the same chromosomes. Synteny results suggested that the studied BACs are placed in a smaller number of chromosomes in S. senegalensis that in other species. Phylogenetic analyses suggested that the evolutionary rate of immune system genes studied is similar among the taxa studied, given that the clustering obtained was in accordance with the accepted phylogenetic relationships among these species. This study contributes to a better understanding of the structure and function of the immune system of the Senegalese sole, which is essential for the development of new technologies and products to improve fish health and productivity.

Differential gene expression following TLR stimulation in rag1-/- mutant zebrafish tissues and morphological descriptions of lymphocyte-like cell populations

In the absence of lymphocytes, rag1-/- mutant zebrafish develop protective immunity to bacteria. In mammals, induction of protection by innate immunity can be mediated by macrophages or natural killer (NK) cells. To elucidate potential responsive cell populations, we morphologically characterized lymphocyte-like cells (LLCs) from liver, spleen and kidney hematopoietic tissues. In fish, these cells include NK cells and Non-specific cytotoxic cells (NCCs). We also evaluated the transcriptional expression response of select genes that are important indicators of NK and macrophage activation after exposure to specific TLR ligands. The LLC cell populations could be discriminated by size and further discriminated by the presence of cytoplasmic granules. Expression levels of mx, tnfα, ifnγ, t-bet and nitr9 demonstrated dynamic changes in response to intra-coelomically administered β glucan (a TLR2/6 ligand), Poly I:C (a TLR3 ligand) and resiquimod (R848) (a TLR7/8 ligand). Following TLR 2/6 stimulation, there was a greater than 100 fold increase in ifnγ in liver, kidney and spleen and moderate increases in tnfα in liver and kidney. TLR3 stimulation caused broad up regulation of mx, down-regulation of tnfα in kidney and spleen tissues and up regulation of nitr9 in the kidney. Following TLR 7/8 stimulation, there was a greater than 100 fold increase in ifnγ in liver and kidney and t-bet in liver. Our gene expression findings suggest that LLCs and macrophages are stimulated following β glucan exposure. Poly I:C causes type I interferon response and mild induction of LLC in the kidney and R-848 exposure causes the strongest LLC stimulation. Overall, the strongest NK like gene expression occurred in the liver. These differential effects of TLR ligands in rag1-/- mutant zebrafish shows strong NK cell-like gene expression responses, especially in the liver, and provides tools to evaluate the basis for protective immunity mediated by the innate immune cells of fish.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Construction of the BAC Library of Small Abalone (Haliotis diversicolor) for Gene Screening and Genome Characterization

The small abalone (Haliotis diversicolor) is one of the most important aquaculture species in East Asia. To facilitate gene cloning and characterization, genome analysis, and genetic breeding of it, we constructed a large-insert bacterial artificial chromosome (BAC) library, which is an important genetic tool for advanced genetics and genomics research. The small abalone BAC library includes 92,610 clones with an average insert size of 120 Kb, equivalent to approximately 7.6× of the small abalone genome. We set up three-dimensional pools and super pools of 18,432 BAC clones for target gene screening using PCR method. To assess the approach, we screened 12 target genes in these 18,432 BAC clones and identified 16 positive BAC clones. Eight positive BAC clones were then sequenced and assembled with the next generation sequencing platform. The assembled contigs representing these 8 BAC clones spanned 928 Kb of the small abalone genome, providing the first batch of genome sequences for genome evaluation and characterization. The average GC content of small abalone genome was estimated as 40.33%. A total of 21 protein-coding genes, including 7 target genes, were annotated into the 8 BACs, which proved the feasibility of PCR screening approach with three-dimensional pools in small abalone BAC library. One hundred fifty microsatellite loci were also identified from the sequences for marker development in the future. The BAC library and clone pools provided valuable resources and tools for genetic breeding and conservation of H. diversicolor.

Transcriptome assembly, gene annotation and tissue gene expression atlas of the rainbow trout

Efforts to obtain a comprehensive genome sequence for rainbow trout are ongoing and will be complemented by transcriptome information that will enhance genome assembly and annotation. Previously, transcriptome reference sequences were reported using data from different sources. Although the previous work added a great wealth of sequences, a complete and well-annotated transcriptome is still needed. In addition, gene expression in different tissues was not completely addressed in the previous studies. In this study, non-normalized cDNA libraries were sequenced from 13 different tissues of a single doubled haploid rainbow trout from the same source used for the rainbow trout genome sequence. A total of ~1.167 billion paired-end reads were de novo assembled using the Trinity RNA-Seq assembler yielding 474,524 contigs > 500 base-pairs. Of them, 287,593 had homologies to the NCBI non-redundant protein database. The longest contig of each cluster was selected as a reference, yielding 44,990 representative contigs. A total of 4,146 contigs (9.2%), including 710 full-length sequences, did not match any mRNA sequences in the current rainbow trout genome reference. Mapping reads to the reference genome identified an additional 11,843 transcripts not annotated in the genome. A digital gene expression atlas revealed 7,678 housekeeping and 4,021 tissue-specific genes. Expression of about 16,000–32,000 genes (35–71% of the identified genes) accounted for basic and specialized functions of each tissue. White muscle and stomach had the least complex transcriptomes, with high percentages of their total mRNA contributed by a small number of genes. Brain, testis and intestine, in contrast, had complex transcriptomes, with a large numbers of genes involved in their expression patterns. This study provides comprehensive de novo transcriptome information that is suitable for functional and comparative genomics studies in rainbow trout, including annotation of the genome.

Endocrinology: advances through omics and related technologies

The rapid development of new omics technologies to measure changes at genetic, transcriptomic, proteomic, and metabolomics levels together with the evolution of methods to analyze and integrate the data at a systems level are revolutionizing the study of biological processes. Here we discuss how new approaches using omics technologies have expanded our knowledge especially in nontraditional models. Our increasing knowledge of these interactions and evolutionary pathway conservation facilitates the use of nontraditional species, both invertebrate and vertebrate, as new model species for biological and endocrinology research. The increasing availability of technology to create organisms overexpressing key genes in endocrine function allows manipulation of complex regulatory networks such as growth hormone (GH) in transgenic fish where disregulation of GH production to produce larger fish has also permitted exploration of the role that GH plays in testis development, suggesting that it does so through interactions with insulin-like growth factors. The availability of omics tools to monitor changes at nearly any level in any organism, manipulate gene expression and behavior, and integrate data across biological levels, provides novel opportunities to explore endocrine function across many species and understand the complex roles that key genes play in different aspects of the endocrine function.

The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates

Vertebrate evolution has been shaped by several rounds of whole-genome duplications (WGDs) that are often suggested to be associated with adaptive radiations and evolutionary innovations. Due to an additional round of WGD, the rainbow trout genome offers a unique opportunity to investigate the early evolutionary fate of a duplicated vertebrate genome. Here we show that after 100 million years of evolution the two ancestral subgenomes have remained extremely collinear, despite the loss of half of the duplicated protein-coding genes, mostly through pseudogenization. In striking contrast is the fate of miRNA genes that have almost all been retained as duplicated copies. The slow and stepwise rediploidization process characterized here challenges the current hypothesis that WGD is followed by massive and rapid genomic reorganizations and gene deletions.

Although whole-genome duplications (WGDs) are rare events, they have an important role in shaping vertebrate evolution. Here, the authors sequence the rainbow trout genome and show that rediploidization after WGD occurs in a slow and stepwise manner.

A first generation BAC-based physical map of the half-smooth tongue sole (Cynoglossus semilaevis) genome

Background

Half-smooth tongue sole (Cynoglossus semilaevis Günther) has been exploited as a commercially important cultured marine flatfish, and female grows 2–3 times faster than male. Genetic studies, especially on the chromosomal sex-determining system of this species, have been carried out in the last decade. Although the genome of half-smooth tongue sole was relatively small (626.9 Mb), there are still some difficulties in the high-quality assembly of the next generation genome sequencing reads without the assistance of a physical map, especially for the W chromosome of this fish due to abundance of repetitive sequences. The objective of this study is to construct a bacterial artificial chromosome (BAC)-based physical map for half-smooth tongue sole with the method of high information content fingerprinting (HICF).

Results

A physical map of half-smooth tongue sole was constructed with 30, 294 valid fingerprints (7.5 × genome coverage) with a tolerance of 4 and an initial cutoff of 1e-60. A total of 29,709 clones were assembled into 1,485 contigs with an average length of 539 kb and a N50 length of 664 kb. There were 394 contigs longer than the N50 length, and these contigs will be a useful resource for future integration with linkage map and whole genome sequence assembly. The estimated physical length of the assembled contigs was 797 Mb, representing approximately 1.27 coverage of the half-smooth tongue sole genome. The largest contig contained 410 BAC clones with a physical length of 3.48 Mb. Almost all of the 676 BAC clones (99.9%) in the 21 randomly selected contigs were positively validated by PCR assays, thereby confirming the reliability of the assembly.

Conclusions

A first generation BAC-based physical map of half-smooth tongue sole was constructed with high reliability. The map will promote genetic improvement programs of this fish, especially integration of physical and genetic maps, fine-mappings of important gene and/or QTL, comparative and evolutionary genomics studies, as well as whole genome sequence assembly.

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.

Assessment of genetic correlation between bacterial cold water disease resistance and spleen index in a domesticated population of rainbow trout: identification of QTL on chromosome Omy19

Selective breeding of animals for increased disease resistance is an effective strategy to reduce mortality in aquaculture. However, implementation of selective breeding programs is limited by an incomplete understanding of host resistance traits. We previously reported results of a rainbow trout selection program that demonstrated increased survival following challenge with Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD). Mechanistic study of disease resistance identified a positive phenotypic correlation between post-challenge survival and spleen somatic-index (SI). Herein, we investigated the hypothesis of a genetic correlation between the two traits influenced by colocalizing QTL. We evaluated the inheritance and calculated the genetic correlation in five year-classes of odd- and even-year breeding lines. A total of 322 pedigreed families (n = 25,369 fish) were measured for disease resistance, and 251 families (n = 5,645 fish) were evaluated for SI. Spleen index was moderately heritable in both even-year (h(2)  = 0.56±0.18) and odd-year (h(2)  = 0.60±0.15) lines. A significant genetic correlation between SI and BCWD resistance was observed in the even-year line (rg  = 0.45±0.20, P = 0.03) but not in the odd-year line (rg  = 0.16±0.12, P = 0.19). Complex segregation analyses of the even-year line provided evidence of genes with major effect on SI, and a genome scan of a single family, 2008132, detected three significant QTL on chromosomes Omy19, 16 and 5, in addition to ten suggestive QTL. A separate chromosome scan for disease resistance in family 2008132 identified a significant BCWD QTL on Omy19 that was associated with time to death and percent survival. In family 2008132, Omy19 microsatellite alleles that associated with higher disease resistance also associated with increased spleen size raising the hypothesis that closely linked QTL contribute to the correlation between these traits. To our knowledge, this is the first estimation of spleen size heritability and evidence for genetic linkage with specific disease resistance in a teleost fish.

A preliminary genetic map in Solea senegalensis (Pleuronectiformes, Soleidae) using BAC-FISH and next-generation sequencing

This article presents the first physical mapping carried out in the Senegalese sole (Solea senegalensis), an important marine fish species of Southern Europe. Eight probes were designated to pick up genes of interest in aquaculture (candidate genes) from a bacterial artificial chromosome (BAC) library using a method of rapid screening based on a 4-dimension PCR. Seven known and 3 unknown clones were isolated and labeled. The 10 BAC clones were used as probes to map the karyotype of the species by fluorescence in situ hybridization (FISH). Nine out of the 10 clones were localized in only 1 chromosome pair, whereas the remaining one hybridized on 2 chromosome pairs. The 2-color FISH experiments showed colocation of 4 probes in 2 chromosome pairs. In addition, 2-color FISH was carried out both with 5S rDNA and the BAC containing the lysozyme gene published previously. This first genetic map of the Senegalese sole represents a starting point for future studies of the sole genome. In addition, 7 out of the 10 BAC clones were sequenced using next-generation sequencing, and bioinformatic characterization of the sequences was carried out. Hence the anchoring of the sequences to specific chromosomes or chromosome arms is now possible, leading to an initial scaffold of the Senegalese sole genome.

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.

Characterization of the OmyY1 Region on the Rainbow Trout Y Chromosome

We characterized the male-specific region on the Y chromosome of rainbow trout, which contains both sdY (the sex-determining gene) and the male-specific genetic marker, OmyY1. Several clones containing the OmyY1 marker were screened from a BAC library from a YY clonal line and found to be part of an 800 kb BAC contig. Using fluorescence in situ hybridization (FISH), these clones were localized to the end of the short arm of the Y chromosome in rainbow trout, with an additional signal on the end of the X chromosome in many cells. We sequenced a minimum tiling path of these clones using Illumina and 454 pyrosequencing. The region is rich in transposons and rDNA, but also appears to contain several single-copy protein-coding genes. Most of these genes are also found on the X chromosome; and in several cases sex-specific SNPs in these genes were identified between the male (YY) and female (XX) homozygous clonal lines. Additional genes were identified by hybridization of the BACs to the cGRASP salmonid 4x44K oligo microarray. By BLASTn evaluations using hypothetical transcripts of OmyY1-linked candidate genes as query against several EST databases, we conclude at least 12 of these candidate genes are likely functional, and expressed.

Meiotic maps of sockeye salmon derived from massively parallel DNA sequencing

Background

Meiotic maps are a key tool for comparative genomics and association mapping studies. Next-generation sequencing and genotyping by sequencing are speeding the processes of SNP discovery and the development of new genetic tools, including meiotic maps for numerous species. Currently there are limited genetic resources for sockeye salmon, Oncorhynchus nerka. We develop the first dense meiotic map for sockeye salmon using a combination of novel SNPs found in restriction site associated DNA (RAD tags) and SNPs available from existing expressed sequence tag (EST) based assays.

Results

We discovered and genotyped putative SNPs in 3,430 RAD tags. We removed paralogous sequence variants leaving 1,672 SNPs; these were combined with 53 EST-based SNP genotypes for linkage mapping. The map contained 29 male and female linkage groups, consistent with the haploid chromosome number expected for sockeye salmon. The female map contains 1,057 loci spanning 4,896 cM, and the male map contains 1,118 loci spanning 4,220 cM. Regions of conservation with rainbow trout and synteny between the RAD based rainbow trout map and the sockeye salmon map were established.

Conclusions

Using RAD sequencing and EST-based SNP assays we successfully generated the first high density linkage map for sockeye salmon.

A second generation integrated map of the rainbow trout (Oncorhynchus mykiss) genome: analysis of conserved synteny with model fish genomes

DNA fingerprints and end sequences from bacterial artificial chromosomes (BACs) from two new libraries were generated to improve the first generation integrated physical and genetic map of the rainbow trout (Oncorhynchus mykiss) genome. The current version of the physical map is composed of 167,989 clones of which 158,670 are assembled into contigs and 9,319 are singletons. The number of contigs was reduced from 4,173 to 3,220. End sequencing of clones from the new libraries generated a total of 11,958 high quality sequence reads. The end sequences were used to develop 238 new microsatellites of which 42 were added to the genetic map. Conserved synteny between the rainbow trout genome and model fish genomes was analyzed using 188,443 BAC end sequence (BES) reads. The fractions of BES reads with significant BLASTN hits against the zebrafish, medaka, and stickleback genomes were 8.8%, 9.7%, and 10.5%, respectively, while the fractions of significant BLASTX hits against the zebrafish, medaka, and stickleback protein databases were 6.2%, 5.8%, and 5.5%, respectively. The overall number of unique regions of conserved synteny identified through grouping of the rainbow trout BES into fingerprinting contigs was 2,259, 2,229, and 2,203 for stickleback, medaka, and zebrafish, respectively. These numbers are approximately three to five times greater than those we have previously identified using BAC paired ends. Clustering of the conserved synteny analysis results by linkage groups as derived from the integrated physical and genetic map revealed that despite the low sequence homology, large blocks of macrosynteny are conserved between chromosome arms of rainbow trout and the model fish species.

Genome physical mapping of polyploids: a BIBAC physical map of cultivated tetraploid cotton, Gossypium hirsutum L

Polyploids account for approximately 70% of flowering plants, including many field, horticulture and forage crops. Cottons are a world-leading fiber and important oilseed crop, and a model species for study of plant polyploidization, cellulose biosynthesis and cell wall biogenesis. This study has addressed the concerns of physical mapping of polyploids with BACs and/or BIBACs by constructing a physical map of the tetraploid cotton, Gossypium hirsutum L. The physical map consists of 3,450 BIBAC contigs with an N50 contig size of 863 kb, collectively spanning 2,244 Mb. We sorted the map contigs according to their origin of subgenome, showing that we assembled physical maps for the A- and D-subgenomes of the tetraploid cotton, separately. We also identified the BIBACs in the map minimal tilling path, which consists of 15,277 clones. Moreover, we have marked the physical map with nearly 10,000 BIBAC ends (BESs), making one BES in approximately 250 kb. This physical map provides a line of evidence and a strategy for physical mapping of polyploids, and a platform for advanced research of the tetraploid cotton genome, particularly fine mapping and cloning the cotton agronomic genes and QTLs, and sequencing and assembling the cotton genome using the modern next-generation sequencing technology.

A synthetic rainbow trout linkage map provides new insights into the salmonid whole genome duplication and the conservation of synteny among teleosts

Background

Rainbow trout is an economically important fish and a suitable experimental organism in many fields of biology including genome evolution, owing to the occurrence of a salmonid specific whole-genome duplication (4^th WGD). Rainbow trout is among some of the most studied teleosts and has benefited from substantial efforts to develop genomic resources (e.g., linkage maps. Here, we first generated a synthetic map by merging segregation data files derived from three independent linkage maps. Then, we used it to evaluate genome conservation between rainbow trout and three teleost models, medaka, stickleback and zebrafish and to further investigate the extent of the 4^th WGD in trout genome.

Results

The INRA linkage map was updated by adding 211 new markers. After standardization of marker names, consistency of marker assignment to linkage groups and marker orders was checked across the three different data sets and only loci showing consistent location over all or almost all of the data sets were kept. This resulted in a synthetic map consisting of 2226 markers and 29 linkage groups spanning over 3600 cM. Blastn searches against medaka, stickleback, and zebrafish genomic databases resulted in 778, 824 and 730 significant hits respectively while blastx searches yielded 505, 513 and 510 significant hits. Homology search results revealed that, for most rainbow trout chromosomes, large syntenic regions encompassing nearly whole chromosome arms have been conserved between rainbow trout and its closest models, medaka and stickleback. Large conserved syntenies were also found between the genomes of rainbow trout and the reconstructed teleost ancestor. These syntenies consolidated the known homeologous affinities between rainbow trout chromosomes due to the 4^th WGD and suggested new ones.

Conclusions

The synthetic map constructed herein further highlights the stability of the teleost genome over long evolutionary time scales. This map can be easily extended by incorporating new data sets and should help future rainbow trout whole genome sequence assembly. Finally, the persistence of large conserved syntenies across teleosts should facilitate the identification of candidate genes through comparative mapping, even if the occurrence of intra-chromosomal micro-rearrangement may hinder the accurate prediction their genomic location.

A BAC-based physical map of Zhikong scallop (Chlamys farreri Jones et Preston)

Zhikong scallop (Chlamys farreri) is one of the most economically important aquaculture species in China. Physical maps are crucial tools for genome sequencing, gene mapping and cloning, genetic improvement and selective breeding. In this study, we have developed a genome-wide, BAC-based physical map for the species. A total of 81,408 clones from two BAC libraries of the scallop were fingerprinted using an ABI 3130xl Genetic Analyzer and a fingerprinting kit developed in our laboratory. After data processing, 63,641 (∼5.8× genome coverage) fingerprints were validated and used in the physical map assembly. A total of 3,696 contigs were assembled for the physical map. Each contig contained an average of 10.0 clones, with an average physical size of 490 kb. The combined total physical size of all contigs was 1.81 Gb, equivalent to approximately 1.5 fold of the scallop haploid genome. A total of 10,587 BAC end sequences (BESs) and 167 markers were integrated into the physical map. We evaluated the physical map by overgo hybridization, BAC-FISH (fluorescence in situ hybridization), contig BAC pool screening and source BAC library screening. The results have provided evidence of the high reliability of the contig physical map. This is the first physical map in mollusc; therefore, it provides an important platform for advanced research of genomics and genetics, and mapping of genes and QTL of economical importance, thus facilitating the genetic improvement and selective breeding of the scallop and other marine molluscs.

Generation of the first BAC-based physical map of the common carp genome

Background

Common carp (Cyprinus carpio), a member of Cyprinidae, is the third most important aquaculture species in the world with an annual global production of 3.4 million metric tons, accounting for nearly 14% of the all freshwater aquaculture production in the world. Apparently genomic resources are needed for this species in order to study its performance and production traits. In spite of much progress, no physical maps have been available for common carp. The objective of this project was to generate a BAC-based physical map using fluorescent restriction fingerprinting.

Result

The first generation of common carp physical map was constructed using four- color High Information Content Fingerprinting (HICF). A total of 72,158 BAC clones were analyzed that generated 67,493 valid fingerprints (5.5 × genome coverage). These BAC clones were assembled into 3,696 contigs with the average length of 476 kb and a N50 length of 688 kb, representing approximately 1.76 Gb of the common carp genome. The largest contig contained 171 BAC clones with the physical length of 3.12 Mb. There are 761 contigs longer than the N50, and these contigs should be the most useful resource for future integrations with linkage map and whole genome sequence assembly. The common carp physical map is available at http://genomics.cafs.ac.cn/fpc/WebAGCoL/Carp/WebFPC/.

Conclusion

The reported common carp physical map is the first physical map of the common carp genome. It should be a valuable genome resource facilitating whole genome sequence assembly and characterization of position-based genes important for aquaculture traits.

Enhancing genome investigations in the mosquito Culex quinquefasciatus via BAC library construction and characterization

BACKGROUND

Culex quinquefasciatus (Say) is a major species in the Culex pipiens complex and an important vector for several human pathogens including West Nile virus and parasitic filarial nematodes causing lymphatic filariasis. It is common throughout tropical and subtropical regions and is among the most geographically widespread mosquito species. Although the complete genome sequence is now available, additional genomic tools are needed to improve the sequence assembly.

FINDINGS

We constructed a bacterial artificial chromosome (BAC) library using the pIndigoBAC536 vector and HindIII partially digested DNA isolated from Cx. quinquefasciatus pupae, Johannesburg strain (NDJ). Insert size was estimated by NotI digestion and pulsed-field gel electrophoresis of 82 randomly selected clones. To estimate genome coverage, each 384-well plate was pooled for screening with 29 simple sequence repeat (SSR) and five gene markers. The NDJ library consists of 55,296 clones arrayed in 144 384-well microplates. Fragment insert size ranged from 50 to 190 kb in length (mean = 106 kb). Based on a mean insert size of 106 kb and a genome size of 579 Mbp, the BAC library provides ~10.1-fold coverage of the Cx. quinquefasciatus genome. PCR screening of BAC DNA plate pools for SSR loci from the genetic linkage map and for four genes associated with reproductive diapause in Culex pipiens resulted in a mean of 9.0 positive plate pools per locus.

CONCLUSION

The NDJ library represents an excellent resource for genome assembly enhancement and characterization in Culex pipiens complex mosquitoes.

A genetically anchored physical framework for Theobroma cacao cv. Matina 1-6

Background

The fermented dried seeds of Theobroma cacao (cacao tree) are the main ingredient in chocolate. World cocoa production was estimated to be 3 million tons in 2010 with an annual estimated average growth rate of 2.2%. The cacao bean production industry is currently under threat from a rise in fungal diseases including black pod, frosty pod, and witches' broom. In order to address these issues, genome-sequencing efforts have been initiated recently to facilitate identification of genetic markers and genes that could be utilized to accelerate the release of robust T. cacao cultivars. However, problems inherent with assembly and resolution of distal regions of complex eukaryotic genomes, such as gaps, chimeric joins, and unresolvable repeat-induced compressions, have been unavoidably encountered with the sequencing strategies selected.

Results

Here, we describe the construction of a BAC-based integrated genetic-physical map of the T. cacao cultivar Matina 1-6 which is designed to augment and enhance these sequencing efforts. Three BAC libraries, each comprised of 10× coverage, were constructed and fingerprinted. 230 genetic markers from a high-resolution genetic recombination map and 96 Arabidopsis-derived conserved ortholog set (COS) II markers were anchored using pooled overgo hybridization. A dense tile path consisting of 29,383 BACs was selected and end-sequenced. The physical map consists of 154 contigs and 4,268 singletons. Forty-nine contigs are genetically anchored and ordered to chromosomes for a total span of 307.2 Mbp. The unanchored contigs (105) span 67.4 Mbp and therefore the estimated genome size of T. cacao is 374.6 Mbp. A comparative analysis with A. thaliana, V. vinifera, and P. trichocarpa suggests that comparisons of the genome assemblies of these distantly related species could provide insights into genome structure, evolutionary history, conservation of functional sites, and improvements in physical map assembly. A comparison between the two T. cacao cultivars Matina 1-6 and Criollo indicates a high degree of collinearity in their genomes, yet rearrangements were also observed.

Conclusions

The results presented in this study are a stand-alone resource for functional exploitation and enhancement of Theobroma cacao but are also expected to complement and augment ongoing genome-sequencing efforts. This resource will serve as a template for refinement of the T. cacao genome through gap-filling, targeted re-sequencing, and resolution of repetitive DNA arrays.

Construction and characterization of the BAC library for common carp Cyprinus carpio L. and establishment of microsynteny with zebrafish Danio rerio

A bacterial artificial chromosome (BAC) library of common carp Cyprinus carpio L. was constructed as a part of ongoing common carp genome project, which is aiming assembly of common carp genome. The library, containing a total of 92,160 BAC clones with an average insert size of 141 kb, was constructed into the restriction site of Hind III on BAC vector CopyControl pCC1BAC, covering 7.7 X haploid genome equivalents. Three dimension pools and superpools of the BAC library were established and 23 positive clones of 14 targets were identified from one-fifth of the BAC library. Pilot project of BAC end sequencing was conducted on 2,688 BAC ends from 1,344 clones and harvested 2,522 high-quality Q20 sequences with average length of 677 bp. The sequencing success rate was 93.8% and pair-end success rate was 92.3%. A total of 212 microsyntenies had been established between common carp and zebrafish genomes as a trial for genome-wide comparative genomics in these two closely related species.

Analysis of BAC-end sequences in rainbow trout: content characterization and assessment of synteny between trout and other fish genomes

BACKGROUND

Rainbow trout (Oncorhynchus mykiss) are cultivated worldwide for aquaculture production and are widely used as a model species to gain knowledge of many aspects of fish biology. The common ancestor of the salmonids experienced a whole genome duplication event, making extant salmonids such as the rainbow trout an excellent model for studying the evolution of tetraploidization and re-diploidization in vertebrates. However, the lack of a reference genome sequence hampers research progress for both academic and applied purposes. In order to enrich the genomic tools already available in this species and provide further insight on the complexity of its genome, we sequenced a large number of rainbow trout BAC-end sequences (BES) and characterized their contents.

RESULTS

A total of 176,485 high quality BES, were generated, representing approximately 4% of the trout genome. BES analyses identified 6,848 simple sequence repeats (SSRs), of which 3,854 had high quality flanking sequences for PCR primers design. The first rainbow trout repeat elements database (INRA RT rep1.0) containing 735 putative repeat elements was developed, and identified almost 59.5% of the BES database in base-pairs as repetitive sequence. Approximately 55% of the BES reads (97,846) had more than 100 base pairs of contiguous non-repetitive sequences. The fractions of the 97,846 non-repetitive trout BES reads that had significant BLASTN hits against the zebrafish, medaka and stickleback genome databases were 15%, 16.2% and 17.9%, respectively, while the fractions of the non-repetitive BES reads that had significant BLASTX hits against the zebrafish, medaka, and stickleback protein databases were 10.7%, 9.5% and 9.5%, respectively. Comparative genomics using paired BAC-ends revealed several regions of conserved synteny across all the fish species analyzed in this study.

CONCLUSIONS

The characterization of BES provided insights on the rainbow trout genome. The discovery of specific repeat elements will facilitate analyses of sequence content (e.g. for SNPs discovery and for transcriptome characterization) and future genome sequence assemblies. The numerous microsatellites will facilitate integration of the linkage and physical maps and serve as valuable resource for fine mapping QTL and positional cloning of genes affecting aquaculture production traits. Furthermore, comparative genomics through BES can be used for identifying positional candidate genes from QTL mapping studies, aid in future assembly of a reference genome sequence and elucidating sequence content and complexity in the rainbow trout genome.

Genomic insight into the common carp (Cyprinus carpio) genome by sequencing analysis of BAC-end sequences

Background

Common carp is one of the most important aquaculture teleost fish in the world. Common carp and other closely related Cyprinidae species provide over 30% aquaculture production in the world. However, common carp genomic resources are still relatively underdeveloped. BAC end sequences (BES) are important resources for genome research on BAC-anchored genetic marker development, linkage map and physical map integration, and whole genome sequence assembling and scaffolding.

Result

To develop such valuable resources in common carp (Cyprinus carpio), a total of 40,224 BAC clones were sequenced on both ends, generating 65,720 clean BES with an average read length of 647 bp after sequence processing, representing 42,522,168 bp or 2.5% of common carp genome. The first survey of common carp genome was conducted with various bioinformatics tools. The common carp genome contains over 17.3% of repetitive elements with GC content of 36.8% and 518 transposon ORFs. To identify and develop BAC-anchored microsatellite markers, a total of 13,581 microsatellites were detected from 10,355 BES. The coding region of 7,127 genes were recognized from 9,443 BES on 7,453 BACs, with 1,990 BACs have genes on both ends. To evaluate the similarity to the genome of closely related zebrafish, BES of common carp were aligned against zebrafish genome. A total of 39,335 BES of common carp have conserved homologs on zebrafish genome which demonstrated the high similarity between zebrafish and common carp genomes, indicating the feasibility of comparative mapping between zebrafish and common carp once we have physical map of common carp.

Conclusion

BAC end sequences are great resources for the first genome wide survey of common carp. The repetitive DNA was estimated to be approximate 28% of common carp genome, indicating the higher complexity of the genome. Comparative analysis had mapped around 40,000 BES to zebrafish genome and established over 3,100 microsyntenies, covering over 50% of the zebrafish genome. BES of common carp are tremendous tools for comparative mapping between the two closely related species, zebrafish and common carp, which should facilitate both structural and functional genome analysis in common carp.

A first generation integrated map of the rainbow trout genome

Background

Rainbow trout (Oncorhynchus mykiss) are the most-widely cultivated cold freshwater fish in the world and an important model species for many research areas. Coupling great interest in this species as a research model with the need for genetic improvement of aquaculture production efficiency traits justifies the continued development of genomics research resources. Many quantitative trait loci (QTL) have been identified for production and life-history traits in rainbow trout. An integrated physical and genetic map is needed to facilitate fine mapping of QTL and the selection of positional candidate genes for incorporation in marker-assisted selection (MAS) programs for improving rainbow trout aquaculture production.

Results

The first generation integrated map of the rainbow trout genome is composed of 238 BAC contigs anchored to chromosomes of the genetic map. It covers more than 10% of the genome across segments from all 29 chromosomes. Anchoring of 203 contigs to chromosomes of the National Center for Cool and Cold Water Aquaculture (NCCCWA) genetic map was achieved through mapping of 288 genetic markers derived from BAC end sequences (BES), screening of the BAC library with previously mapped markers and matching of SNPs with BES reads. In addition, 35 contigs were anchored to linkage groups of the INRA (French National Institute of Agricultural Research) genetic map through markers that were not informative for linkage analysis in the NCCCWA mapping panel. The ratio of physical to genetic linkage distances varied substantially among chromosomes and BAC contigs with an average of 3,033 Kb/cM.

Conclusions

The integrated map described here provides a framework for a robust composite genome map for rainbow trout. This resource is needed for genomic analyses in this research model and economically important species and will facilitate comparative genome mapping with other salmonids and with model fish species. This resource will also facilitate efforts to assemble a whole-genome reference sequence for rainbow trout.

Characterization of the rainbow trout transcriptome using Sanger and 454-pyrosequencing approaches

Background

Rainbow trout are important fish for aquaculture and recreational fisheries and serves as a model species for research investigations associated with carcinogenesis, comparative immunology, toxicology and evolutionary biology. However, to date there is no genome reference sequence to facilitate the development of molecular technologies that utilize high-throughput characterizations of gene expression and genetic variation. Alternatively, transcriptome sequencing is a rapid and efficient means for gene discovery and genetic marker development. Although a large number (258,973) of EST sequences are publicly available, the nature of rainbow trout duplicated genome hinders assembly and complicates annotation.

Results

High-throughput deep sequencing of the Swanson rainbow trout doubled-haploid transcriptome using 454-pyrosequencing technology yielded ~1.3 million reads with an average length of 344 bp, a total of 447 million bases. De novo assembly of the sequences yielded 151,847 Tentative Consensus (TC) sequences (average length of 662 bp) and 224,391 singletons. A combination assembly of both the 454-pyrosequencing ESTs and the pre-existing sequences resulted in 161,818 TCs (average length of 758 bp) and 261,071 singletons. Gene Ontology analysis of the combination assembly showed high similarities to transcriptomes of other fish species with known genome sequences.

Conclusion

The 454 library significantly increased the suite of ESTs available for rainbow trout, allowing improved assembly and annotation of the transcriptome. Furthermore, the 454 sequencing enables functional genome research in rainbow trout, providing a wealth of sequence data to serve as a reference transcriptome for future studies including identification of paralogous sequences and/or allelic variation, digital gene expression and proteomic research.

A first generation BAC-based physical map of the Asian seabass (Lates calcarifer)

Background

The Asian seabass (Lates calcarifer) is an important marine foodfish species in Southeast Asia and Australia. Genetic improvement of this species has been achieved to some extent through selective breeding programs since 1990s. Several genomic tools such as DNA markers, a linkage map, cDNA and BAC libraries have been developed to assist selective breeding. A physical map is still lacking, although it is essential for positional cloning of genes located in quantitative trait loci (QTL) and assembly of whole genome sequences.

Methodology/Principal Findings

A genome-wide physical map of the Asian seabass was constructed by restriction fingerprinting of 38,208 BAC clones with SNaPshot HICF FPC technique. A total of 30,454 were assembled into 2,865 contigs. The physical length of the assembled contigs summed up to 665 Mb. Analyses of some contigs using different methods demonstrated the reliability of the assembly.

Conclusions/Significance

The present physical map is the first physical map for Asian seabass. This physical map will facilitate the fine mapping of QTL for economically important traits and the positional cloning of genes located in QTL. It will also be useful for the whole genome sequencing and assembly. Detailed information about BAC-contigs and BAC clones are available upon request.

Generation of a BAC-based physical map of the melon genome

BACKGROUND

Cucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has high intra-specific genetic variation, morphologic diversity and a small genome size (450 Mb), which make this species suitable for a great variety of molecular and genetic studies that can lead to the development of tools for breeding varieties of the species. A number of genetic and genomic resources have already been developed, such as several genetic maps and BAC genomic libraries. These tools are essential for the construction of a physical map, a valuable resource for map-based cloning, comparative genomics and assembly of whole genome sequencing data. However, no physical map of any Cucurbitaceae has yet been developed. A project has recently been started to sequence the complete melon genome following a whole-genome shotgun strategy, which makes use of massive sequencing data. A BAC-based melon physical map will be a useful tool to help assemble and refine the draft genome data that is being produced.

RESULTS

A melon physical map was constructed using a 5.7 x BAC library and a genetic map previously developed in our laboratories. High-information-content fingerprinting (HICF) was carried out on 23,040 BAC clones, digesting with five restriction enzymes and SNaPshot labeling, followed by contig assembly with FPC software. The physical map has 1,355 contigs and 441 singletons, with an estimated physical length of 407 Mb (0.9 x coverage of the genome) and the longest contig being 3.2 Mb. The anchoring of 845 BAC clones to 178 genetic markers (100 RFLPs, 76 SNPs and 2 SSRs) also allowed the genetic positioning of 183 physical map contigs/singletons, representing 55 Mb (12%) of the melon genome, to individual chromosomal loci. The melon FPC database is available for download at http://melonomics.upv.es/static/files/public/physical_map/.

CONCLUSIONS

Here we report the construction of the first physical map of a Cucurbitaceae species described so far. The physical map was integrated with the genetic map so that a number of physical contigs, representing 12% of the melon genome, could be anchored to known genetic positions. The data presented is already helping to improve the quality of the melon genomic sequence available as a result of a project currently being carried out in Spain, adopting a whole genome shotgun approach based on 454 sequencing data.