Families of clustered microsatellites in a conifer genome

Development of nuclear microsatellite loci for Pinus albicaulis Engelm. (Pinaceae), a conifer of conservation concern

Pinus albicaulis (whitebark pine) is a widely-distributed but rapidly declining high elevation western North American tree and a candidate for listing under the U.S. Endangered Species Act. Our objectives were to develop reliable nuclear microsatellite markers that can be used to assess within-population genetic diversity as well as seed and pollen migration dynamics, and to validate markers using two geographically proximal P. albicaulis populations. We identified 1,667 microsatellite-containing sequences from shotgun DNA libraries of P. albicaulis. Primer pairs were designed for 308 unique microsatellite-containing loci, and these were evaluated for PCR amplification success and segregation in a panel of diploid needle tissue. DNA was extracted with an SDS protocol, and primers were screened through gel electrophoresis. Microsatellites were genotyped through fluorescent primer fragment analysis. Ten novel and 13 transferred loci were found to be reproducible in analyses based on 20 foliage samples from each of two locations: Henderson Mountain, Custer Gallatin National Forest, Montana, and Mt. Washburn, Yellowstone National Park, Wyoming (USA). Transferred loci had higher numbers of alleles and expected heterozygosities than novel loci, but also revealed evidence for a higher frequency of null alleles. Eight of the 13 transferred loci deviated significantly from Hardy-Weinberg Equilibrium, and showed large positive F_IS values that were likely inflated by null alleles. Mantel’s tests of transferred and novel markers showed no correlation between genetic and geographic distances within or among the two sampled populations. AMOVA suggests that 91% of genetic variability occurs within populations and 9% between the two populations. Studies assessing genetic diversity using these microsatellite loci can help guide future management and restoration activities for P. albicaulis.

Insights into the loblolly pine genome: characterization of BAC and fosmid sequences

Despite their prevalence and importance, the genome sequences of loblolly pine, Norway spruce, and white spruce, three ecologically and economically important conifer species, are just becoming available to the research community. Following the completion of these large assemblies, annotation efforts will be undertaken to characterize the reference sequences. Accurate annotation of these ancient genomes would be aided by a comprehensive repeat library; however, few studies have generated enough sequence to fully evaluate and catalog their non-genic content. In this paper, two sets of loblolly pine genomic sequence, 103 previously assembled BACs and 90,954 newly sequenced and assembled fosmid scaffolds, were analyzed. Together, this sequence represents 280 Mbp (roughly 1% of the loblolly pine genome) and one of the most comprehensive studies of repetitive elements and genes in a gymnosperm species. A combination of homology and de novo methodologies were applied to identify both conserved and novel repeats. Similarity analysis estimated a repetitive content of 27% that included both full and partial elements. When combined with the de novo investigation, the estimate increased to almost 86%. Over 60% of the repetitive sequence consists of full or partial LTR (long terminal repeat) retrotransposons. Through de novo approaches, 6,270 novel, full-length transposable element families and 9,415 sub-families were identified. Among those 6,270 families, 82% were annotated as single-copy. Several of the novel, high-copy families are described here, with the largest, PtPiedmont, comprising 133 full-length copies. In addition to repeats, analysis of the coding region reported 23 full-length eukaryotic orthologous proteins (KOGS) and another 29 novel or orthologous genes. These discoveries, along with other genomic resources, will be used to annotate conifer genomes and address long-standing questions about gymnosperm evolution.

Development of SSR and gene-targeted markers for construction of a framework linkage map of Catharanthus roseus

BACKGROUND AND AIMS

Catharanthus roseus is a plant of great medicinal importance, yet inadequate knowledge of its genome structure and the unavailability of genomic resources have been major impediments in the development of improved varieties. The aims of this study were to develop co-dominant sequence-tagged microsatellite sites (STMS) and gene-targeted markers (GTMs) and utilize them for the construction of a framework intraspecific linkage map of C. roseus.

METHODS

For simple sequence repeat (SSR) isolation, a genomic library enriched for (GA)(n) repeats was constructed from C. roseus 'Nirmal' (CrN1). In addition, GTMs were also designed from 12 genes of the TIA (terpenoid indole alkaloid) pathway - the medicinally most significant pathway in C. roseus. An F(2) mapping population was also generated by crossing two diverse accessions of C. roseus CrN1 (Nirmal)×CrN82 (Kew).

KEY RESULTS

A new set of 314 STMS markers and 64 GTMs were developed in this study. A segregating F(2) mapping population consisting of 111 F(2) individuals was generated. For generating the linkage map, a set of 423 co-dominant markers (378 newly developed and 45 published earlier) were screened for polymorphism between the parental genotypes, of which 134 were identified to be polymorphic. A total of 114 markers were mapped on eight linkage groups that spanned a 632·7 cM region of the genome with an average marker distance of 5·55 cM. Further, the mechanism of hypervariability at the gene-targeted loci was investigated at the sequence level.

CONCLUSIONS

For the first time, a large array of STMS markers and GTMs was generated in the model medicinal plant C. roseus. Moreover, the first microsatellite marker-based linkage map was described in this study. Together, these will serve as a foundation for future genomics studies related to quantitative trait loci analysis and molecular breeding in C. roseus.

Advancing the STMS genomic resources for defining new locations on the intraspecific genetic linkage map of chickpea (Cicer arietinum L.)

Background

Chickpea (Cicer arietinum L.) is an economically important cool season grain legume crop that is valued for its nutritive seeds having high protein content. However, several biotic and abiotic stresses and the low genetic variability in the chickpea genome have continuously hindered the chickpea molecular breeding programs. STMS (Sequence Tagged Microsatellite Sites) markers which are preferred for the construction of saturated linkage maps in several crop species, have also emerged as the most efficient and reliable source for detecting allelic diversity in chickpea. However, the number of STMS markers reported in chickpea is still limited and moreover exhibit low rates of both inter and intraspecific polymorphism, thereby limiting the positions of the SSR markers especially on the intraspecific linkage maps of chickpea. Hence, this study was undertaken with the aim of developing additional STMS markers and utilizing them for advancing the genetic linkage map of chickpea which would have applications in QTL identification, MAS and for de novo assembly of high throughput whole genome sequence data.

Results

A microsatellite enriched library of chickpea (enriched for (GT/CA)_n and (GA/CT)_n repeats) was constructed from which 387 putative microsatellite containing clones were identified. From these, 254 STMS primers were designed of which 181 were developed as functional markers. An intraspecific mapping population of chickpea, [ICCV-2 (single podded) × JG-62 (double podded)] and comprising of 126 RILs, was genotyped for mapping. Of the 522 chickpea STMS markers (including the double-podding trait, screened for parental polymorphism, 226 (43.3%) were polymorphic in the parents and were used to genotype the RILs. At a LOD score of 3.5, eight linkage groups defining the position of 138 markers were obtained that spanned 630.9 cM with an average marker density of 4.57 cM. Further, based on the common loci present between the current map and the previously published chickpea intraspecific map, integration of maps was performed which revealed improvement of marker density and saturation of the region in the vicinity of sfl (double-podding) gene thereby bringing about an advancement of the current map.

Conclusion

An arsenal of 181 new chickpea STMS markers was reported. The developed intraspecific linkage map defined map positions of 138 markers which included 101 new locations.Map integration with a previously published map was carried out which revealed an advanced map with improved density. This study is a major contribution towards providing advanced genomic resources which will facilitate chickpea geneticists and molecular breeders in developing superior genotypes with improved traits.

The Pinus taeda genome is characterized by diverse and highly diverged repetitive sequences

BACKGROUND

In today's age of genomic discovery, no attempt has been made to comprehensively sequence a gymnosperm genome. The largest genus in the coniferous family Pinaceae is Pinus, whose 110-120 species have extremely large genomes (c. 20-40 Gb, 2N = 24). The size and complexity of these genomes have prompted much speculation as to the feasibility of completing a conifer genome sequence. Conifer genomes are reputed to be highly repetitive, but there is little information available on the nature and identity of repetitive units in gymnosperms. The pines have extensive genetic resources, with approximately 329000 ESTs from eleven species and genetic maps in eight species, including a dense genetic map of the twelve linkage groups in Pinus taeda.

RESULTS

We present here the Sanger sequence and annotation of ten P. taeda BAC clones and Genome Analyzer II whole genome shotgun (WGS) sequences representing 7.5% of the genome. Computational annotation of ten BACs predicts three putative protein-coding genes and at least fifteen likely pseudogenes in nearly one megabase of sequence. We found three conifer-specific LTR retroelements in the BACs, and tentatively identified at least 15 others based on evidence from the distantly related angiosperms. Alignment of WGS sequences to the BACs indicates that 80% of BAC sequences have similar copies (> or = 75% nucleotide identity) elsewhere in the genome, but only 23% have identical copies (99% identity). The three most common repetitive elements in the genome were identified and, when combined, represent less than 5% of the genome.

CONCLUSIONS

This study indicates that the majority of repeats in the P. taeda genome are 'novel' and will therefore require additional BAC or genomic sequencing for accurate characterization. The pine genome contains a very large number of diverged and probably defunct repetitive elements. This study also provides new evidence that sequencing a pine genome using a WGS approach is a feasible goal.

Microsatellites for the genus Cucurbita and an SSR-based genetic linkage map of Cucurbita pepo L

Until recently, only a few microsatellites have been available for Cucurbita, thus their development is highly desirable. The Austrian oil-pumpkin variety Gleisdorfer Olkürbis (C. pepo subsp. pepo) and the C. moschata cultivar Soler (Puerto Rico) were used for SSR development. SSR-enriched partial genomic libraries were established and 2,400 clones were sequenced. Of these 1,058 (44%) contained an SSR at least four repeats long. Primers were designed for 532 SSRs; 500 primer pairs produced fragments of expected size. Of these, 405 (81%) amplified polymorphic fragments in a set of 12 genotypes: three C. moschata, one C. ecuadorensis, and eight C. pepo representing all eight cultivar groups. On an average, C. pepo and C. moschata produced 3.3 alleles per primer pair, showing high inter-species transferability. There were 187 SSR markers detecting polymorphism between the USA oil-pumpkin variety "Lady Godiva" (O5) and the Italian crookneck variety "Bianco Friulano" (CN), which are the parents of our previous F(2) mapping population. It has been used to construct the first published C. pepo map, containing mainly RAPD and AFLP markers. Now the updated map comprises 178 SSRs, 244 AFLPs, 230 RAPDs, five SCARs, and two morphological traits (h and B). It contains 20 linkage groups with a map density of 2.9 cM. The observed genome coverage (Co) is 86.8%.

Evolution of a triplet repeat in a conifer

The opportunity to trace the evolution of a triplet repeat is rare, especially for seed-plant lineages with a well-defined fossil record. Microsatellite PtTX2133 sequences from 18 species in 2 conifer genera were used to calibrate the birth of a CAGn repeat, from its protomicrosatellite origins to its repeat expansion. Birth occurred in the hard-pine genome ~ 136 million years ago, or 14 million generations ago, then expanded as a polymorphic triplet repeat 136-100 million years before a major North American vicariance event. Calibration of the triplet-repeat birth and expansion is supported by the shared allelic lineages among Old and New World hard pines and the shared alleles solely among North American diploxylon or hard pines. Five CAGn repeat units appeared to be the expansion threshold for Old and New World diploxylon pines. Haploxylon pine species worldwide did not undergo birth and repeat expansion, remaining monomorphic, with a single imperfect 198-bp allele. A sister genus, Picea, had only a region of cryptic simplicity, preceding a proto-microsatellite region. The polymorphic triplet repeat in hard pines is older than some long-lived microsatellites reported for reptiles, yet younger than those reported for insects. Some cautionary points are raised about phylogenetic applications for this long-lived microsatellite.

Development of simple sequence repeat (SSR) markers and construction of an SSR-based linkage map in Italian ryegrass (Lolium multiflorum Lam.)

In order to develop simple sequence repeat (SSR) markers in Italian ryegrass, we constructed a genomic library enriched for (CA)n-containing SSR repeats. A total of 1,544 clones were sequenced, of which 1,044 (67.6%) contained SSR motifs, and 395 unique clones were chosen for primer design. Three hundred and fifty-seven of these clones amplified products of the expected size in both parents of a two-way pseudo-testcross F(1) mapping population, and 260 primer pairs detected genetic polymorphism in the F(1) population. Genetic loci detected by a total of 218 primer pairs were assigned to locations on seven linkage groups, representing the seven chromosomes of the haploid Italian ryegrass karyotype. The SSR markers covered 887.8 cM of the female map and 795.8 cM of the male map. The average distance between two flanking SSR markers was 3.2 cM. The SSR markers developed in this study will be useful in cultivar discrimination, linkage analysis, and marker-assisted selection of Italian ryegrass and closely related species.

An enhanced microsatellite map of diploid Fragaria

A total of 45 microsatellites (SSRs) were developed for mapping in Fragaria. They included 31 newly isolated codominant genomic SSRs from F. nubicola and a further 14 SSRs, derived from an expressed sequence tagged library (EST-SSRs) of the cultivated strawberry, F. x ananassa. These, and an additional 64 previously characterised but unmapped SSRs and EST-SSRs, were scored in the diploid Fragaria interspecific F2 mapping population (FVxFN) derived from a cross between F. vesca 815 and F. nubicola 601. The cosegregation data of these 109 SSRs, and of 73 previously mapped molecular markers, were used to elaborate an enhanced linkage map. The map is composed of 182 molecular markers (175 microsatellites, six gene specific markers and one sequence-characterised amplified region) and spans 424 cM over seven linkage groups. The average marker spacing is 2.3 cM/marker and the map now contains just eight gaps longer than 10 cM. The transferability of the new SSR markers to the cultivated strawberry was demonstrated using eight cultivars. Because of the transferable nature of these markers, the map produced will provide a useful reference framework for the development of linkage maps of the cultivated strawberry and for the development of other key resources for Fragaria such as a physical map. In addition, the map now provides a framework upon which to place transferable markers, such as genes of known function, for comparative mapping purposes within Rosaceae.

Lepidopteran microsatellite DNA: redundant but promising

The development of microsatellite DNA markers for use with butterflies and moths is extremely difficult for no apparent reason. New findings demonstrate that many lepidopteran microsatellite sequences exist in multiple copies in the genome, and have similar or almost identical flanking regions. These findings provide a compelling explanation for the low efficiency of microsatellite isolation in Lepidoptera, and might also shed light on the evolutionary dynamics of microsatellite sequences in these insects and other organisms.

Cross-amplification and sequence variation of microsatellite loci in Eurasian hard pines

Microsatellite transfer across coniferous species is a valued methodology because de novo development for each species is costly and there are many species with only a limited commodity value. Cross-species amplification of orthologous microsatellite regions provides valuable information on mutational and evolutionary processes affecting these loci. We tested 19 nuclear microsatellite markers from Pinus taeda L. (subsection Australes) and three from P. sylvestris L. (subsection Pinus) on seven Eurasian hard pine species ( P. uncinata Ram., P. sylvestris L., P. nigra Arn., P. pinaster Ait., P. halepensis Mill., P. pinea L. and P. canariensis Sm.). Transfer rates to species in subsection Pinus (36-59%) were slightly higher than those to subsections Pineae and Pinaster (32-45%). Half of the trans-specific microsatellites were found to be polymorphic over evolutionary times of approximately 100 million years (ten million generations). Sequencing of three trans-specific microsatellites showed conserved repeat and flanking regions. Both a decrease in the number of perfect repeats in the non-focal species and a polarity for mutation, the latter defined as a higher substitution rate in the flanking sequence regions close to the repeat motifs, were observed in the trans-specific microsatellites. The transfer of microsatellites among hard pine species proved to be useful for obtaining highly polymorphic markers in a wide range of species, thereby providing new tools for population and quantitative genetic studies.

A full saturated linkage map of Picea abies including AFLP, SSR, ESTP, 5S rDNA and morphological markers

Based on an F(1) progeny of 73 individuals, two parental maps were constructed according to the double pseudo-test cross strategy. The paternal map contained 16 linkage groups for a total genetic length of 1,792 cM. The maternal map covered 1,920 cM, and consisted of 12 linkage groups. These parental maps were then integrated using 66 intercross markers. The resulting consensus map covered 2,035 cM and included 755 markers (661 AFLPs, 74 SSRs, 18 ESTPs, the 5S rDNA and the early cone formation trait) on 12 linkage groups, reflecting the haploid number of chromosomes of Picea abies. The average spacing between two adjacent markers was 2.6 cM. The presence of 39 of the SSR and/or ESTP markers from this consensus map on other published maps of different Picea and Pinus species allowed us to establish partial linkage group homologies across three P. abies maps (up to five common markers per linkage group). This first saturated linkage map of P. abies could be therefore used as a support for developing comparative genome mapping in conifers.

Detection and physical mapping of the 18S-5.8S-26S rDNA and the pKFJ660 probe with microsatellite sequences derived from the rice blast fungus (Magnaporthe grisea) in conifer species

Sequences homologous to the pKFJ660 probe, a fragment of DNA derived from the rice blast fungus (Magnaporthe grisea) carrying TC/AG repeat microsatellite sequences and 30 bp direct repeats were identified in the genome of Picea (spruce) and Pinus (pine) species by fluorescence in situ hybridization (FISH) and slot blot analyses. Slot blot analysis using the pKFJ660 probe revealed hybridization signals with genomic DNAs from various pine and spruce species. Further analyses indicated that the copy number of the (AG)30 motif was higher than 5 x 10(4) per plant genome for all plant samples tested, but the copy number of the sequences homologous to the whole pKFJ660 probe varies considerably among the 25 plant species tested. In situ hybridization of metaphase chromosomes from Pinus resinosa, P. banksiana and P. strobus showed the presence of sequences homologous to this probe on several chromosomes in a dispersed pattern. Major signals were observed on a few chromosomes indicating that some of these sequences are clustered in specific genomic locations. The locations of these repeats were compared to those of 18S-5.8S-26S rDNA in pine species. Chromosomal distribution of 18S-5.8S-26S rDNA varied among the three pine species (P. resinosa, P. banksiana and P. strobus) studied. Ribosomal DNA (rDNA) sites were identified on 14 to 20 chromosomes in these pine species.

No clustering for linkage map based on low-copy and undermethylated microsatellites

Clustering has been reported for conifer genetic maps based on hypomethylated or low-copy molecular markers, resulting in uneven marker distribution. To test this, a framework genetic map was constructed from three types of microsatellites: low-copy, undermethylated, and genomic. These Pinus taeda L. microsatellites were mapped using a three-generation pedigree with 118 progeny. The microsatellites were highly informative; of the 32 markers in intercross configuration, 29 were segregating for three or four alleles in the progeny. The sex-averaged map placed 51 of the 95 markers in 15 linkage groups at LOD > 4.0. No clustering or uneven distribution across the genome was observed. The three types of P. taeda microsatellites were randomly dispersed within each linkage group. The 51 microsatellites covered a map distance of 795 cM, an average distance of 21.8 cM between markers, roughly half of the estimated total map length. The minimum and maximum distances between any two bins was 4.4 and 45.3 cM, respectively. These microsatellites provided anchor points for framework mapping for polymorphism in P. taeda and other closely related hard pines.Key words: simple sequence repeats, genome duplication, conifers, Pinus taeda L. hypomethylation, low-copy kinetic component.

Comparing EST-based genetic maps between Pinus sylvestris and Pinus taeda

A genetic map of Pinus sylvestris was constructed using ESTP (expressed sequence tag polymorphism) markers and other gene-based markers, AFLP markers and microsatellites. Part of the ESTP markers (40) were developed and mapped earlier in Pinus taeda, and additional markers were generated based on P. sylvestris sequences or sequences from other pine species. The mapping in P. sylvestris was based on 94 F(1) progeny from a cross between plus-tree parents E635C and E1101. AFLP framework maps for the parent trees were first constructed. The ESTP and other gene sequence-based markers were added to the framework maps, as well as five published microsatellite loci. The separate maps were then integrated with the aid of AFLPs segregating in both trees (dominant segregation ratios 3:1) as well as gene markers and microsatellites segregating in both parent trees (segregation ratios 1:1:1:1 or 1:2:1). The integrated map consisted of 12 groups corresponding to the P. taeda linkage groups, and additionally three and six smaller groups for E1101 and E635C, respectively. The number of framework AFLP markers in the integrated map is altogether 194 and the number of gene markers 61. The total length of the integrated map was 1,314 cM. The set of markers developed for P. sylvestris was also added to existing maps of two P. taeda pedigrees. Starting with a mapped marker from one pedigree in the source species resulted in a mapped marker in a pedigree of the other species in more than 40% of the cases, with about equal success in both directions. The maps of the two species are largely colinear, even if the species have diverged more than 70 MYA. Most cases of different locations were probably due to problems in identifying the orthologous members of gene families. These data provide a first ESTP-containing map of P. sylvestris, which can also be used for comparing this species to additional species mapped with the same markers.

Undermethylated DNA as a source of microsatellites from a conifer genome

Developing microsatellites from the large, highly duplicated conifer genome requires special tools. To improve the efficiency of developing Pinus taeda L. microsatellites, undermethylated (UM) DNA fragments were used to construct a microsatellite-enriched copy library. A methylation-sensitive restriction enzyme, McrBC, was used to enrich for UM DNA before library construction. Digested DNA fragments larger than 9 kb were then excised and digested with RsaI and used to construct nine dinucleotide and trinucleotide libraries. A total of 1016 microsatellite-positive clones were detected among 11 904 clones and 620 of these were unique. Of 245 primer sets that produced a PCR product, 113 could be developed as UM microsatellite markers and 70 were polymorphic. Inheritance and marker informativeness were tested for a random sample of 36 polymorphic markers using a three-generation outbred pedigree. Thirty-one microsatellites (86%) had single-locus inheritance despite the highly duplicated nature of the P. taeda genome. Nineteen UM microsatellites had highly informative intercross mating type configurations. Allele number and frequency were estimated for eleven UM microsatellites using a population survey. Allele numbers for these UM microsatellites ranged from 3 to 12 with an average of 5.7 alleles/locus. Frequencies for the 63 alleles were mostly in the low-common range; only 14 of the 63 were in the rare allele (q < 0.05) class. Enriching for UM DNA was an efficient method for developing polymorphic microsatellites from a large plant genome.