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.     

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.     

A deer (subfamily Cervinae) genetic linkage map and the evolution of ruminant genomes

Comparative maps between ruminant species and humans are increasingly important tools for the discovery of genes underlying economically important traits. In this article we present a primary linkage map of the deer genome derived from an interspecies hybrid between red deer (Cervus elaphus) and Père David's deer (Elaphurus davidianus). The map is approximately 2500 cM long and contains >600 markers including both evolutionary conserved type I markers and highly polymorphic type II markers (microsatellites). Comparative mapping by annotation and sequence similarity (COMPASS) was demonstrated to be a useful tool for mapping bovine and ovine ESTs in deer. Using marker order as a phylogenetic character and comparative map information from human, mouse, deer, cattle, and sheep, we reconstructed the karyotype of the ancestral Pecoran mammal and identified the chromosome rearrangements that have occurred in the sheep, cattle, and deer lineages. The deer map and interspecies hybrid pedigrees described here are a valuable resource for (1) predicting the location of orthologs to human genes in ruminants, (2) mapping QTL in farmed and wild deer populations, and (3) ruminant phylogenetic studies.     

A comprehensive radiation hybrid map of the bovine genome comprising 5593 loci

A bovine whole genome 7000-rad radiation hybrid (RH) panel, SUNbRH(7000-rad), was constructed to build a high-resolution RH map. The Shirakawa-USDA linkage map served as a scaffold to construct a framework map of 3216 microsatellites on which 2377 ESTs were ordered. The resulting RH map provided essentially complete coverage across the genome, with 1 cR7000 corresponding to 114 kb, and a cattle-human comparative map of 1716 bovine genes and sequences annotated in the human genome, which covered 79 and 72% of the bovine and human genomes, respectively. We then integrated the bovine RH and comparative maps with BAC fingerprint information in to construct a detailed, BAC-based physical map covering a reported 40-cM quantitative trait locus region for intramuscular fat or "marbling" on BTA 4. In summary, the new, high-resolution SUNbRH7000-rad, comparative, Shirakawa-USDA linkage, and BAC fingerprint maps provide a set of genomic tools for fine mapping regions of interest in cattle.     

Improved resolution of the comparative horse-human map: investigating markers with in silico and linkage mapping approaches

Genetic maps are extremely important tools for tracing the genes that govern economically significant traits, and microsatellites are a significant component of these. In this study, we isolated 2346 novel horse microsatellites as resources for the construction of high-density horse genetic maps. Of these 2346 markers, 339 (14.5%) horse sequences showed sequence homology to DNA sequences in the human genome, demonstrating that microsatellites as type II markers are valuable resources for developing linkage maps and that they have a potential equal to that of type I markers for developing comparative maps. Of the 339 markers, 206 (60.8%) were assigned to horse chromosomes using the Animal Health Trust (AHT) full-sib reference family, and 195 (94.6%) of these localized to the expected syntenic locations on the human genome. These results confirmed the high level of accuracy of in silico mapping. Thus, the 339 markers that exhibited homology to the human genome increased the density of markers on the horse-human comparative map. The resulting comparative map will facilitate the use of horse microsatellites as genetic markers for the identification of quantitative trait loci (QTL) that have been mapped on the human genome. In addition, although the in silico and linkage mapping data did not agree for the other 11 (5.4%) of the assigned 206 markers, these may represent new putative regions of horse-human synteny.     

Physical and linkage mapping of mammary-derived expressed sequence tags in cattle

This study describes the physical and linkage mapping of 42 gene-associated markers developed from mammary gland-derived expressed sequence tags to the cattle genome. Of the markers, 25 were placed on the USDA reference linkage map and 37 were positioned on the Roslin 3000-rad radiation hybrid (RH) map, with 20 assignments shared between the maps. Although no novel regions of conserved synteny between the cattle and the human genomes were identified, the coverage was extended for bovine chromosomes 3, 7, 15, and 29 compared with previously published comparative maps between human and bovine genomes. Overall, these data improve the resolution of the human-bovine comparative maps and will assist future efforts to integrate bovine RH and linkage map data.     

Bioinformatic Mining of Type I Microsatellites from Expressed Sequence Tags of Channel Catfish (<Emphasis Type="Italic">Ictalurus punctatus</Emphasis>)

Gene-derived markers are pivotal to the analysis of genome structure, organization, and evolution and necessary for comparative genomics. However, gene-derived markers are relatively difficult to develop. This project utilized the genomic resources of channel catfish expressed sequence tags (ESTs) to identify simple sequence repeats (SSRs), or microsatellites. It took the advantage of ESTs for the establishment of gene identities, and of microsatellites for the acquisition of high polymorphism. When microsatellites are tagged to genes, the microsatellites can then be used as gene markers. A bioinformatic analysis of 43,033 ESTs identified 4855 ESTs containing microsatellites. Cluster analysis indicated that 1312 of these ESTs fell into 569 contigs, and the remaining 3534 ESTs were singletons. A total of 4103 unique microsatellite-containing genes were identified. The dinucleotide CA/TG and GA/TC pairs were the most abundant microsatellites. AT-rich microsatellite types were predominant among trinucleotide and tetranucleotide microsatellites, consistent with our earlier estimation that the catfish genome is highly AT-rich. Our preliminary results indicated that the majority of the identified microsatellites were polymorphic and, therefore, useful for genetic linkage mapping of catfish. Mapping of these gene-derived markers is under way, which will set the foundation for comparative genome analysis in catfish.     

A high-resolution comparative RH map of the proximal part of bovine chromosome 1

Current comparative maps between human chromosome 21 and the proximal part of cattle chromosome 1 are insufficient to define chromosomal rearrangements because of the low density of mapped genes in the bovine genome. The recently completed sequence of human chromosome 21 facilitates the detailed comparative analysis of corresponding segments on BTA1. In this study eight bovine bacterial artificial chromosome (BAC) clones containing bovine orthologues of human chromosome 21 genes, i.e. GRIK1, CLDN8, TIAM1, HUNK, SYNJ1, OLIG2, IL10RB, and KCNE2 were physically assigned by fluorescence in situ hybridization (FISH) to BTA1q12.1-q12.2. Sequence tagged site (STS) markers derived from these clones were mapped on the 3000 rad Roslin/Cambridge bovine radiation hybrid (RH) panel. In addition to these eight novel markers, 17 known markers from previously published BTA1 linkage or RH maps were also mapped on the Roslin/Cambridge bovine RH panel resulting in an integrated map with 25 markers of 355.4 cR(3000) length. The human-cattle genome comparison revealed the existence of three chromosomal breakpoints and two probable inversions in this region.     

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.     

Use of bovine EST data and human genomic sequences to map 100 gene-specific bovine markers

A system to use bovine EST data in conjunction with human genomic sequence to improve the bovine linkage map over the entire genome or on specific chromosomes was evaluated. Bovine EST sequence was used to provide primer sequences corresponding to bovine genes, while human genomic sequence directed primer design to flank introns and produce amplicons of appropriate size for efficient direct sequencing. The sequence tagged sites (STS) produced in this way from the four sires of the MARC reference families were examined for single nucleotide polymorphisms (SNPs) that could be used to map the corresponding genes. With this approach, along with a primer/extension mass spectrometry SNP genotyping assay, 100 ESTs were placed on the bovine genetic linkage map. The first 70 were chosen at random from bovine EST–human genomic comparisons. An additional 30 ESTs were successfully mapped to bovine Chromosome 19 (BTA19), and comparison of the resulting BTA19 map to the position of the corresponding human orthologs on the HSA17 draft sequences revealed differences in the spacing and order of genes. Over 80% of successful amplicons contained SNPs, indicating that this is an efficient approach to generating EST-associated genetic markers. We have demonstrated the feasibility of constructing a linkage map based on SNPs associated with ESTs and the plausibility of utilizing EST, comparative mapping information, and human sequence data to target regions of the bovine genome for SNP marker development.     

Refinement of bovine chromosome 2 linkage map near the mh locus reveals rearrangements between the bovine and human genomes

The locus responsible for the appearance of muscular hypertrophy (mh) in double muscled cattle breeds has recently been shown to encode a secreted growth factor designated myostatin (MSTN). This conclusion was based in part on the placement of MSTN in the interval to which mh had been mapped on bovine chromosome 2 (BTA2). During the mapping phase of the study, numerous yeast artificial chromosome (YAC) clones were isolated that contained genetic markers closely linked to mh. Other YACs and cosmids were identified that contained genes selected from human chromosome 2q (HSA2q), with the goal of defining the position of breakpoints in conserved synteny between the bovine and human comparative maps, thereby permitting accurate selection of positional candidate genes. An efficient subcloning procedure was developed to obtain microsatellites (ms) from YAC clones, to increase the number of informative meioses in herds segregating for mh. The same procedure was used to place the human orthologues of engrailed-1 (EN1), interleukin 1 beta (IL1B), and paired-box-containing 8 (PAX8) genes on the cattle map to further define the positions of breakpoints in conserved synteny and gene order. Twenty-three of 28 ms identified from YAC subclone libraries were informative in the mapping families. Seven mapped to the centromeric end of BTA2, which contains the mh locus, improving marker density and informativeness. The two MSTN and four EN1 gene-associated ms markers developed from YACs, map to positions 1·5 and 61·6 cm in the BTA2 linkage group, respectively. In addition, ms markers developed from cosmids containing either IL1B or PAX8, map to positions 56·6 and 56·9 cm in the BTA11 linkage group, respectively. These linkage data confirm the location and orientation of orthologous segments of HSA2q that were previously indistinguishable on the bovine map, and demonstrates the presence of microrearrangements of gene order (segments <10 cm ) and conserved synteny between the human and bovine genomes.     

An expanded comparative map of bovine chromosome 27 targeting dairy form QTL regions

At present, the density of genes on the bovine maps is extremely limited and current resolution of the human-bovine comparative map is insufficient for selection of candidate genes controlling many economic traits of interest in dairy cattle. This study describes the chromosomal mapping of 10 selected gene-associated markers to bovine linkage and radiation hybrid maps to improve the breakpoint resolution in the human-bovine comparative map near two previously identified quantitative trait loci for the linear type trait, dairy form. Two regions of conserved synteny not previously described are reported between the telomeric region of bovine chromosome 27 (BTA27) and human chromosome 3 (HSA3) p24 region and between the HSA4q34.1 region and BTA8. These data increase the number of genes positioned on the bovine gene maps, refine the human-bovine comparative map, and should improve the efficiency of candidate gene selection for the dairy form trait in cattle.     

Directed cosmid isolation of bovine markers for physical assignment by fish

Abstract

A directed PCR‐based iterative screening protocol was developed to isolate cosmids containing microsatellite markers linked to chromosomal regions of interest, such as those near the ends of linkage groups or quantitative trait loci. This method was optimized for large‐scale screening of total genomic libraries and used to purify bovine cosmids that anchor the ends of the bovine linkage group 28 (BTA28). Cosmids containing ms markers RBP3 and BMS2060 were purified for fluorescence in situ hybridization and assigned to 28ql8‐ql9 and 28ql2, respectively. These assignments indicated that approximately 73% of BTA28 (90% excluding the centromere) is covered by the current linkage map. Since this method is applicable to any target gene sequence suitable for PCR amplification, it may be extended to comparative mapping of genes.     

Comparative mapping of bovine chromosome 13 by fluorescence in situ hybridization

We present chromosomal fluorescence in situ hybridization (FISH) results that both extend the HSA20/BTA13 comparative map as well as cytogenetically anchor two microsatellite markers. A bovine bacterial artificial chromosome (BAC) library was screened for conserved genes (type I loci) previously assigned to HSA10 or HSA20 and BTA13, and for microsatellites selected from two published BTA13 linkage maps. Clones from six out of nine comparative loci and both microsatellites were found represented in the BAC library. These BAC clones were used as probes in single colour FISH to determine the chromosome band position of each locus. As predicted by the human/bovine comparative map, all type I loci mapped to BTA13. Because single colour FISH analysis revealed that the loci were clustered within the distal half of BTA13, dual colour FISH was used to confirm the locus order. Established order was centromere- PRNP-(SODIL/AVP/OXT)-(BL42/GNAS1)-HCK-CSSM30 . The findings confirm the presence of a conserved HSA20 homologous synteny group on BTA13 distal of a HSA10 homologous segment.     

Development of transcript-associated microsatellite markers for diversity and linkage mapping studies in hop (<Emphasis Type="Italic">Humulus lupulus</Emphasis> L.)

Data mining of gene sequences available from various projects dealing with the development of expressed sequence tags (ESTs) can contribute to the discovery of new microsatellite markers. Our aim was to develop new microsatellite markers in hop isolated from an enriched cDNA library and from coding GenBank sequences and to test their suitability in hop diversity studies and for construction of a linkage map. In a set of 614 coding GenBank sequences, 72 containing microsatellites were found (11.7%); the most frequent were trinucleotide repeats (54.0%) followed by dinucleotide repeats (34.5%). Additionally, 11 sequences containing microsatellites were isolated from an enriched cDNA library. A total of 34 primer pairs were designed, 29 based on GenBank sequences and five on sequences from the cDNA enriched library. Twenty-seven (79.4%) coding microsatellites were successfully amplified and used in diversity and linkage mapping studies. Eleven primer pairs amplified 12 coding microsatellite loci suitable for mapping and were placed on female and male linkage maps. We were able to extend previous simple sequence repeat (SSR) female, male and integral maps by 38.8, 25.8 and 40.0 cM, respectively. In the diversity study, 36 diverse hop genotypes were analyzed. Twenty-four coding microsatellites were polymorphic, 17 showing co-dominant behavior and 7 primer pairs amplifying three or more bands in some hop genotypes. Altogether, 143 microsatellite DNA fragments were amplified and they revealed a clear separation of hop genotypes according to geographical region, use or breeding history. In addition, a discussion and comparison of results with other plant coding/EST SSR studies is presented. Our results showed that these microsatellite markers can enhance hop diversity and linkage mapping studies and are a comparable marker system to non-coding SSRs.     

Rapid evolution of mammalian X-linked testis microRNAs

Background

Microsatellites in cDNA are useful as molecular markers because they represent transcribed genes and can be used as anchor markers for linkage and comparative mapping, as well as for studying genome evolution. Microsatellites in cDNA can be detected in existing ESTs by data mining. However, in most fish species, no ESTs are available or the number of ESTs is limited, although fishes represent half of the vertebrates on the earth. We developed a simple and efficient method for isolation of microsatellites from cDNA in fish.

Results

The method included normalization of 150 ng cDNA using 0.5 U duplex-specific nuclease (DSN) at 65°C for 30 min, enrichment of microsatellites using biotinylated oligonucleotides and magnetic field, and directional cloning of cDNA into a vector. We tested this method to enrich CA- and GA-microsatellites from cDNA of Asian seabass, and demonstrated that enrichment of microsatellites from normalized cDNA could increased the efficiency of microsatellite isolation over 30 times as compared to direct sequencing of clones from cDNA libraries. One hundred and thirty-nine (36.2%) out of 384 clones from normalized cDNA contained microsatellites. Unique microsatellite sequences accounted for 23.6% (91/384) of sequenced clones. Sixty microsatellites isolated from cDNA were characterized, and 41 were polymorphic. The average allele number of the 41 microsatellites was 4.85 ± 0.54, while the expected heterozygosity was 0.56 ± 0.03. All the isolated microsatellites inherited in a Mendelian pattern.

Conclusion

Normalization of cDNA substantially increased the efficiency of enrichment of microsatellites from cDNA. The described method for isolation of microsatellites from cDNA has the potential to be applied to a wide range of fish species. The microsatellites isolated from cDNA could be useful for linkage and comparative mapping, as well as for studying genome evolution.