Putative SNP discovery in interspecific hybrids of catfish by comparative EST analysis

In this study, we identified putative SNP markers within genes by comparative analysis of expressed sequence tags (ESTs). Comparison of 849 ESTs from blue catfish (Ictalurus furcatus) with >11,000 ESTs from channel catfish (I. punctatus) deposited in GenBank resulted in the identification of 1020 putative SNPs within 161 genes, of which 145 were nuclear genes of known function. The observed frequency of SNPs within ESTs of the two closely related catfish species was 1.32 SNP per 100 bp. The majority of identified SNPs differed between the two species and, therefore, these SNPs are useful for mapping genes in channel catfish x blue catfish interspecific resource families. The SNPs that differed within species were also observed; these can be applied to genome scans in channel catfish resource families.     

Development of Polymorphic EST Markers Suitable for Genetic Linkage Mapping of Catfish

Expressed sequence tag (EST) markers are important for gene mapping and for marker-assisted selection (MAS). To develop EST markers for use in catfish gene mapping, 100 randomly picked complementary DNAs from the channel catfish (Ictalurus punctatus) pituitary library were sequenced. The EST sequences were used to design primers to amplify channel catfish and blue catfish (I. furcatus) genomic DNAs. Polymerase chain reaction products of the ESTs were analyzed to determine length polymorphism between the channel catfish and blue catfish. Eleven polymorphic EST markers were identified. Five of the 11 EST markers were from known genes and the other six were from unidentified ESTs. Seven ESTs were found to be associated with microsatellite sequences. Analysis of channel catfish gene sequences indicated highly biased codon usage, with 16 codons being preferably used. These codons were more preferably used in highly expressed ribosomal protein genes and in highly expressed pituitary hormone genes. G/C-rich codons are less used in channel catfish than those in other vertebrates suggesting AT-richness of the channel catfish genome. Received June 29, 1998; accepted March 29, 1999.     

Characterization of EST-derived microsatellites for gene mapping and evolutionary genomics in turbot

The detection of microsatellite sequences within expressed sequence tags (ESTs) connects potential markers with specific genes, generating type I markers. We have developed and mapped by linkage analysis a set of EST-derived microsatellites in the turbot, Scophthalmus maximus. One hundred and ninety-one microsatellites were identified from 9256 turbot ESTs. Primer design was possible with 98 microsatellites. After genotyping 25 wild turbot and the parents of two reference families for linkage analysis, 43 EST-derived microsatellites were selected because they met technical and polymorphism criteria. A final set of 31 EST-derived microsatellites could be mapped to 17 linkage groups of the turbot consensus map based on 242 anonymous microsatellites. Twenty-four microsatellite-containing ESTs were functionally annotated, confirming them as type I markers. Nineteen were mapped in the turbot consensus map. These EST-derived microsatellites constitute useful tools for genome scanning of turbot populations, marker-assisted selection programmes and comparative mapping.     

Development of polymorphic expressed sequence tag-derived microsatellites for the extension of the genetic linkage map of the black tiger shrimp (Penaeus monodon)

In this study, microsatellite markers were developed for the genetic linkage mapping and breeding program of the black tiger shrimp Penaeus monodon. A total of 997 unique microsatellite-containing expressed sequence tags (ESTs) were identified from 10 100 EST sequences in the P. monodon EST database. AT-rich microsatellite types were predominant in the EST sequences. Homology searching by the blastn and blastx programs revealed that these 997 ESTs represented 8.6% known gene products, 27.8% hypothetical proteins and 63.6% unknown gene products. Characterization of 50 markers on a panel of 35-48 unrelated shrimp indicated an average number of alleles of 12.6 and an average polymorphic information content of 0.723. These EST microsatellite markers along with 208 other markers (185 amplified fragment length polymorphisms, one exon-primed intron-crossing, six single strand conformation polymorphisms, one single nucleotide polymorphism, 13 non-EST-associated microsatellites and two EST-associated microsatellites) were analysed across the international P. monodon mapping family. A total of 144 new markers were added to the P. monodon maps, including 36 of the microsatellite-containing ESTs. The current P. monodon male and female linkage maps have 47 and 36 linkage groups respectively with coverage across half the P. monodon genome.     

An in silico mining for simple sequence repeats from expressed sequence tags of zebrafish, medaka, Fundulus, and Xiphophorus

Teleost fish genome projects involving model species are resulting in a rapid accumulation of genomic and expressed DNA sequences in public databases. The expressed sequence tags (ESTs) collected in the databases can be mined for the analysis of both structural and functional genomics. In this study, we in silico analyzed 49,430 unigenes representing a total of 692,654 ESTs from four model fish for their potential use in developing simple sequence repeats (SSRs), or microsatellites. After bioinformatical mining, a total of 3,018 EST derived SSRs (EST-SSRs) were identified for 2,335 SSR containing ESTs (SSR-ESTs). The frequency of identified SSR-ESTs ranged from 1.5% for Xiphophorus to 7.3% for zebrafish. The dinucleotide repeat motif is the most abundant SSR, accounting for 47%, 52%, 64%, and 78% for medaka, Fundulus, zebrafish, and Xiphophorus, respectively. Simulation analysis suggests that a majority of these EST-SSRs have sufficient flanking sequences for polymerase chain reaction (PCR) primer design. Comparative DNA sequence analyses of SSR-ESTs identified several cross-species SSRs and sequences that may be used as cross-reference genes in comparative studies. For example, the flanking sequences of one SSR (CTG)n within the pituitary tumor-transforming gene (PTTG) 1 interacting protein (PTTGIP), showed conservation spanning the medaka, Fundulus, human, and mouse genomes. This study provides a large body of information on EST-SSRs that can be useful for the development of polymorphic markers, gene mapping, and comparative genome analysis. Functional analysis of these SSR-ESTs may reveal their role in metabolism and gene evolution of these model species.     

Construction of Genetic Linkage Maps and Comparative Genome Analysis of Catfish Using Gene-Associated Markers

A genetic linkage map of the channel catfish genome (N = 29) was constructed using EST-based microsatellite and single nucleotide polymorphism (SNP) markers in an interspecific reference family. A total of 413 microsatellites and 125 SNP markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 allowed mapping of 331 markers (259 microsatellites and 72 SNPs) to 29 linkage groups. Each linkage group contained 3–18 markers. The largest linkage group contained 18 markers and spanned 131.2 cM, while the smallest linkage group contained 14 markers and spanned only 7.9 cM. The linkage map covered a genetic distance of 1811 cM with an average marker interval of 6.0 cM. Sex-specific maps were also constructed; the recombination rate for females was 1.6 times higher than that for males. Putative conserved syntenies between catfish and zebrafish, medaka, and Tetraodon were established, but the overall levels of genome rearrangements were high among the teleost genomes. This study represents a first-generation linkage map constructed by using EST-derived microsatellites and SNPs, laying a framework for large-scale comparative genome analysis in catfish. The conserved syntenies identified here between the catfish and the three model fish species should facilitate structural genome analysis and evolutionary studies, but more importantly should facilitate functional inference of catfish genes. Given that determination of gene functions is difficult in nonmodel species such as catfish, functional genome analysis will have to rely heavily on the establishment of orthologies from model species.     

Characterization of a BAC Library from Channel Catfish <Emphasis Type="Italic">Ictalurus punctatus</Emphasis>: Indications of High Levels of Chromosomal Reshuffling Among Teleost Genomes

The CHORI-212 bacterial artificial chromosome (BAC) library was constructed by cloning EcoRI/EcoRI partially digested DNA into the pTARBAC2.1 vector. The library has an average insert size of 161 kb, and provides 10.6-fold coverage of the channel catfish haploid genome. Screening of 32 genes using overgo or cDNA probes indicated that this library had a good representation of the genome as all tested genes existed in the library. We previously reported sequencing of approximately 25,000 BAC ends that generated 20,366 high-quality BAC end sequences (BES) and identified a large number of sequences similar to known genes using BLASTX searches. In this work, particular attention was given to identification of BAC mate pairs with known genes from both ends. When identified, comparative genome analysis was conducted to determine syntenic regions of the catfish genome with the genomes of zebrafish and Tetraodon. Of the 141 mate pairs with known genes from channel catfish, conserved syntenies were identified in 34 (24.1%), with 30 conserved in the zebrafish genome and 14 conserved in the Tetraodon genome. Additional analysis of three of the 34 conserved syntenic groups by direct sequencing indicated conserved gene contents in all three species. This indicates that comparative genome analysis may provide shortcuts to genome analysis in catfish, especially for short genomic regions once the conserved syntenies are identified. Shaolin Wang and Peng Xu contributed equally to the article.     

Characterization of AT-rich microsatellites in common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.)

Polymorphism of microsatellite markers is often associated with the simple sequence repeat motif targeted. AT-rich microsatellites tend to be highly variable and this appears to be notable, especially in legume genomes. To analyze the value of AT-rich microsatellites for common bean (Phaseolus vulgaris L.), we developed a total of 85 new microsatellite markers, 74 of which targeted ATA or other AT-rich motif loci and 11 of which were made for GA, CA or CAC motif loci. We evaluated the loci for the level of allelic diversity in comparison to previously characterized microsatellites using a panel of 18 standard genotypes and genetically mapped any loci polymorphic in the DOR364 × G19833 population. The majority of the microsatellites produced single bands and detected single loci, however, 15 of the AT-rich microsatellites produced multiple or double banding patterns; while only one of the GA or CA-rich microsatellites did. The polymorphism information content (PIC) values averaged 0.892 and 0.600 for the AT and ATA motif microsatellites, respectively, but only 0.140 for the CA-rich microsatellites. GA microsatellites, which had a large average number of repeats, had high to intermediate PIC, averaging 0.706. A total of 45 loci could be genetically mapped and distribution of the loci across the genome was skewed towards non-distal locations with a greater prevalence of loci on linkage groups b02, b09 and b11. AT-rich microsatellites were found to be a useful source of polymorphic markers for mapping and diversity assessment in common bean that appears to uncover higher diversity than other types of simple sequence repeat markers.     

Isolation and Enrichment of Abundant Microsatellites from a Channel Catfish (Ictalurus punctatus) Brain cDNA Library

Efforts to construct a genetic linkage map of channel catfish have involved identification of random genomic microsatellite markers, as well as anchored Type I loci (expressed genes) from channel catfish. To identify Type I markers we constructed a directional cDNA library from brain tissue to obtain expressed catfish sequences that could be used for single nucleotide polymorphism (SNP) marker development. These cDNA sequences surprisingly contained a high proportion of microsatellites (about 14%) in noncoding regions of expressed sequence tags (ESTs), many of which were not associated with known sequences. To further identify cDNAs with microsatellites and reduce the number of sequencing reactions needed for marker development, we enriched this library for repeat sequences and sequenced clones from both directions. A total of 1644 clones from seven repeat-enriched captures (CA, GT, CT, GA, MTT, TAG, and TAC) were sequenced from both ends, and 795 nonredundant clones were assembled. Thirty-seven percent of the clones contained microsatellites in the trimmed sequence. After assembly in the TIGR Catfish Gene Index (CfGI), 154 contigs matched known vertebrate genes and 92 contigs contained microsatellites. When BLAST-matched orthologues were available for similarity alignments, 28% of these contigs contained repeats in the 5'-UTR, 72% contained repeats in the 3'-UTR, and 8% contained repeats at both ends. Using biotinylated repeat oligonucleotides coupled with streptavidin-coated magnetic beads, and rapid, single-pass hybridization, we were able to enrich our plasmid library greater than two-fold for repeat sequences and increase the ability to link these ESTs with known sequences greater than six-fold.     

Isolation and enrichment of abundant microsatellites from a channel catfish (Ictalurus punctatus) brain cDNA library

Efforts to construct a genetic linkage map of channel catfish have involved identification of random genomic microsatellite markers, as well as anchored Type I loci (expressed genes) from channel catfish. To identify Type I markers we constructed a directional cDNA library from brain tissue to obtain expressed catfish sequences that could be used for single nucleotide polymorphism (SNP) marker development. These cDNA sequences surprisingly contained a high proportion of microsatellites (about 14%) in noncoding regions of expressed sequence tags (ESTs), many of which were not associated with known sequences. To further identify cDNAs with microsatellites and reduce the number of sequencing reactions needed for marker development, we enriched this library for repeat sequences and sequenced clones from both directions. A total of 1644 clones from seven repeat-enriched captures (CA, GT, CT, GA, MTT, TAG, and TAC) were sequenced from both ends, and 795 nonredundant clones were assembled. Thirty-seven percent of the clones contained microsatellites in the trimmed sequence. After assembly in the TIGR Catfish Gene Index (CfGI), 154 contigs matched known vertebrate genes and 92 contigs contained microsatellites. When BLAST-matched orthologues were available for similarity alignments, 28% of these contigs contained repeats in the 5'-UTR, 72% contained repeats in the 3'-UTR, and 8% contained repeats at both ends. Using biotinylated repeat oligonucleotides coupled with streptavidin-coated magnetic beads, and rapid; single-pass hybridization, we were able to enrich our plasmid library greater than two-fold for repeat sequences and increase the ability to link these ESTs with known sequences greater than six-fold.     

Development of EST derived SSRs and SNPs as a genomic resource in Indian catfish, <Emphasis Type="Italic">Clarias batrachus</Emphasis>

Clarias batrachus, an Indian catfish species, is endemic to the Indian subcontinent and potential cultivable species. The genomic resources in C. batrachus in the form of ESTs containing microsatellite repeats (EST-SSR) and single nucleotide polymorphisms (SNPs) that are associated with the expressed genes from spleen were mined. From a total of 1,937 ESTs generated, 1,698 unique sequences were obtained, out of which 221 EST-SSRs were identified and 54% could be functionally annotated by similarity searches. A total of 23 contigs containing 3 or more ESTs were found to contain 31 SNP loci, out of which 8 ESTs showed similarity to genes of known function and 1 for hypothetical protein. Nine ESTs with SSRs and/or SNPs identified in this study were reported to be associated with diseases in human and animals. These identified loci can be developed into markers in C. batrachus, which can be useful in linkage mapping, comparative genomics studies and for its genetic improvement programmes.     

Development of gene-based markers and construction of an integrated linkage map in eggplant by using Solanum orthologous (SOL) gene sets

We constructed an integrated DNA marker linkage map of eggplant (Solanum melongena L.) using DNA marker segregation data sets obtained from two independent intraspecific F(2) populations. The linkage map consisted of 12 linkage groups and encompassed 1,285.5 cM in total. We mapped 952 DNA markers, including 313 genomic SSR markers developed by random sequencing of simple sequence repeat (SSR)-enriched genomic libraries, and 623 single-nucleotide polymorphisms (SNP) and insertion/deletion polymorphisms (InDels) found in eggplant-expressed sequence tags (ESTs) and related genomic sequences [introns and untranslated regions (UTRs)]. Because of their co-dominant inheritance and their highly polymorphic and multi-allelic nature, the SSR markers may be more versatile than the SNP and InDel markers for map-based genetic analysis of any traits of interest using segregating populations derived from any intraspecific crosses of practical breeding materials. However, we found that the distribution of microsatellites in the genome was biased to some extent, and therefore a considerable part of the eggplant genome was first detected when gene-derived SNP and InDel markers were mapped. Of the 623 SNP and InDel markers mapped onto the eggplant integrated map, 469 were derived from eggplant unigenes contained within Solanum orthologous (SOL) gene sets (i.e., sets of orthologous unigenes from eggplant, tomato, and potato). Out of the 469 markers, 326 could also be mapped onto the tomato map. These common markers will be informative landmarks for the transfer of tomato's more saturated genomic information to eggplant and will also provide comparative information on the genome organization of the two solanaceous species. The data are available from the DNA marker database of vegetables, VegMarks (http://vegmarks.nivot.affrc.go.jp).     

Characterization and radiation hybrid mapping of expressed sequence tags from the canine brain

Abstract Maps of the canine genome are now developing rapidly. Most of the markers on the current integrated canine radiation hybrid/genetic linkage/cytogenetic map are highly polymorphic microsatellite (type II) markers that are very useful for mapping disease loci. However, there is still an urgent need for the mapping of gene-based (type I) markers that are required for comparative mapping, as well as identifying candidate genes for disease loci that have been genetically mapped. We constructed an adult brain cDNA library as a resource to increase the number of gene-based markers on the canine genome map. Eighty-one percent of the 2700 sequenced expressed sequence tags (ESTs) represented unique sequences. The canine brain ESTs were compared with sequences in public databases to identify putative canine orthologs of human genes. One hundred nine of the canine ESTs were mapped on the latest canine radiation hybrid (RH) panel to determine the location of the respective canine gene. The addition of these new gene-based markers revealed three conserved segments (CS) between human and canine genomes previously detected by fluorescence in situ hybridization (FISH), but not by RH mapping. In addition, five new CS between dog and human were identified that had not been detected previously by RH mapping or FISH. This work has increased the number of gene-based markers on the canine RH map by approximately 30% and indicates the benefit to be gained by increasing the gene content of the current canine comparative map.     

A gene map of the rat derived from linkage analysis and related regions in the mouse and human genomes

We report the localization by linkage analysis in the rat genome of 148 new markers derived from 128 distinct known gene sequences, ESTs, and anonymous sequences selected in GenBank database on the basis of the presence of a repeated element. The composite linkage map of the rat contributed by our group integrates mapping information on a total of 370 different known genes, ESTs, and anonymous mouse or human sequences, and provides a valuable tool for comparative genome analysis. 206 and 254 homologous loci were identified in the mouse and human genomes respectively. Our linkage map, which combines both anonymous markers and gene markers, should facilitate the advancement of genetic studies for a wide variety of rat models characterized for complete phenotypes. The comparative genome mapping should define genetic regions in human likely to be homologous to susceptibility loci identified in rat and provide useful information for the identification of new potential candidates for genetic disorders. Received: 2 January 1999 / Accepted: 7 March 1999     

A first generation microsatellite- and SNP-based linkage map of Jatropha

Jatropha curcas is a potential plant species for biodiesel production. However, its seed yield is too low for profitable production of biodiesel. To improve the productivity, genetic improvement through breeding is essential. A linkage map is an important component in molecular breeding. We established a first-generation linkage map using a mapping panel containing two backcross populations with 93 progeny. We mapped 506 markers (216 microsatellites and 290 SNPs from ESTs) onto 11 linkage groups. The total length of the map was 1440.9 cM with an average marker space of 2.8 cM. Blasting of 222 Jatropha ESTs containing polymorphic SSR or SNP markers against EST-databases revealed that 91.0%, 86.5% and 79.2% of Jatropha ESTs were homologous to counterparts in castor bean, poplar and Arabidopsis respectively. Mapping 192 orthologous markers to the assembled whole genome sequence of Arabidopsis thaliana identified 38 syntenic blocks and revealed that small linkage blocks were well conserved, but often shuffled. The first generation linkage map and the data of comparative mapping could lay a solid foundation for QTL mapping of agronomic traits, marker-assisted breeding and cloning genes responsible for phenotypic variation.