Microsatellite Content in 397 Nuclear Exons and Their Flanking Regions in the Fern Family Ophioglossaceae

Microsatellites or SSRs are small tandem repeats that are 1-6 bp long. They are usually highly polymorphic and form important portions of genomes. They have been extensively analyzed in humans, animals and model plants; however, information from non-flowering plants is generally lacking. Here, we examined 29 samples of Ophioglossaceae ferns, mainly from the genera Botrychium and Sceptridium. We analyzed the SSR distribution, density and composition in almost 400 nuclear exons and their flanking regions. We detected 45 SSRs in exons and 1475 SSRs in the flanking regions. In the exons, only di-, tri- and tetranucleotides were found, and all of them were 12 bp long. The annotation of the exons containing SSRs showed that they were related to various processes, such as metabolism, catalysis, transportation or plant growth. The flanking regions contained SSRs from all categories, with the most numerous being dinucleotides, followed by tetranucleotides. More than one-third of all the SSRs in the flanking regions were 12 bp long. The SSR densities in the exons were very low, ranging from 0 to 0.07 SSRs/kb, while those in the flanking regions ranged from 0.24 to 0.81 SSRs/kb; and those in the combined dataset ranged from 0.2 to 0.81 SSRs/kb. The majority of the detected SSRs in the flanking regions were polymorphic and present at the same loci across two or more samples but differing in the number of repeats. The SSRs detected here may serve as a basis for further population genetic, phylogenetic or evolutionary genetic studies, as well as for further studies focusing on SSRs in the genomes and their roles in adaptation, evolution and diseases.

Genetic Diversity Maximization as a Strategy for Resilient Forest Ecosystems: A Case Study on Norway Spruce

Norway spruce, economically and ecologically one of the most important European forest tree species, rapidly declines due to massive bark beetle outbreaks across many countries. As a prerequisite of ecosystem stability facing climate changes of uncertain predictions, the reforestation management promoting locally adapted resources of broad genetic diversity should be prioritized, especially in nature conservation areas. In our case study carried out in the national park, Krkonoše Mountains (the Giant Mountains, the Czech Republic), we demonstrated a tree breeding strategy aiming at maximizing genetic diversity. More than four hundred unique Norway spruce accessions were genotyped on 15 microsatellite loci (Ne = 5.764, I = 1.713 and He = 0.685). Two core collection selection approaches were proposed to establish a new deployment population providing local gene sources of high genetic diversity. Namely, the Core Hunter selection algorithm, with average entry-to-nearest-entry distance (EN) optimization, was applied to identify the most diverse core collection set with the highest genetic diversity parameters obtained for 57 selected individuals (Ne = 6.507, I = 1.807, and He = 0.731). The latter core collection method proposed is innovative, based on choosing appropriate genotypes from a clustered heatmap. For simplicity, we demonstrated the principle of selection strategy on a reduced dataset. It is vital to promote panmixia of a newly established production population from a core collection to complete the conservation breeding effort. Thus, we demonstrated the utilization of the Optimum Neighborhood Algorithm (ONA) deployment that outperformed other deployment algorithms, especially in the case of balanced clone representation and uneven shapes of planting plots. We believe that the case study presented can be generalized and considered as a guideline for analogical tree breeding intentions.

Development and characterization of Novel EST-based single-copy genic microsatellite DNA markers in white spruce and black spruce

Due mainly to large genome size and prevalence of repetitive sequences in the nuclear genome of spruce (Picea Mill.), it is very difficult to develop single-copy genomic microsatellite markers. We have developed and characterized 25 polymorphic, single-copy genic microsatellites from white spruce (Picea glauca (Moench) Voss) EST sequences and determined their informativeness in white spruce and black spruce (Picea mariana (Mill.) B.S.P.) and inheritance in black spruce. White spruce EST sequences from NCBI dbEST were searched for the presence of microsatellite repeats. Forty-seven sequences containing dinucleotide, trinucleotide, tetranucleotide and compound repeats were selected to develop primers. Twenty-five of the designed primer pairs yielded scorable amplicons, with single-locus patterns, and were characterized in 20 individuals each of white spruce and black spruce. All 25 microsatellites were polymorphic in white spruce and 24 in black spruce. The number of alleles at a locus ranged from two to 18, with a mean of 8.8 in white spruce, and from one to 17, with a mean of 7.6 in black spruce. The expected heterozygosity/polymorphic information content ranged from 0.10 to 0.92, with a mean of 0.67 in white spruce, and from 0 to 0.93, with a mean of 0.59 in black spruce. Microsatellites with dinucleotide and compound repeats were more informative than those with trinucleotide and tetranucleotide repeats. Eighteen microsatellite markers polymorphic between the parents of a black spruce controlled cross inherited in a single-locus Mendelian fashion. The microsatellite markers developed can be applied for various genetics, genomics, breeding, and conservation studies and applications.

Genetic Structure of Norway Spruce Ecotypes Studied by SSR Markers

Norway spruce is a widespread and economically highly important tree species in Central Europe which occurs there in different morphotypic forms (also known as ecotypes). Previously established common garden experiments indicated that the morphological differentiation is most likely genetically determined. The genetic structure of Norway spruce morphological variants might be an indicator (marker) of specific sustainability in forest ecosystems. In this study, we investigated 436 individuals from autochthonous populations belonging to three different ecotypes. The main aim was to evaluate a level of genetic intra and interpopulation diversity among the low, medium and high-elevation ecotypes using both expressed sequence tag simple sequence repeats (EST – SSR) and genomic SSR markers. Sixteen highly polymorphic microsatellite loci folded in two newly designed multiplexes were used to depicture the genetic structure of targeted trees. Important allele frequency parameters, such as the mean expected (0.722, SE = 0.061) and observed (0.585, SE = 0.062) heterozygosity and mean effective number of alleles (Ne = 5.943, SE = 1.279), were estimated. The low genetic differentiation detected among different ecotypes (Fst = 0.008) was further discussed and clarified.

Use of Nuclear Microsatellite Loci for Evaluating Genetic Diversity of Selected Populations of Picea abies (L.) Karsten in the Czech Republic

DNA polymorphism at nine nuclear microsatellites of nine selected naturally-regenerated Norway spruce populations growing mainly within gene conservation units in different parts of the Czech Republic was studied. To verify the genetic quality of the selected gene conservation unit, we analyzed nine Norway spruce subpopulations from gene conservation unit GZ 102–Orlické hory. Genetic parameters can be used in state administrative decision making on including stands into gene conservation units. The level of genetic diversity within 17 investigated Czech Norway spruce units was relatively high. Mean values for the number of different alleles ranged from 12.2 (population SM 08) to 16.2 (subpopulation SM T4). The values of observed heterozygosity (Ho) ranged from 0.65 to 0.80 and expected heterozygosity (He) from 0.74 to 0.81. Pairwise population FST values ranging from 0.006 to 0.027 indicated low genetic differentiation between units, and values of Nei’s genetic distance among Norway spruce units ranged from 0.046 to 0.168, thus structuring of the investigated Norway spruce units was confirmed. Closer genetic similarity was seen in subpopulations from the gene conservation unit in Orlické hory than in the studied populations from other genetic conservation units. Additionally, the populations SM 01 and SM 05, both of Hurst ecotypes, were the closest to one another and the populations of mountain and alpine ecotypes were assembled into another group.

The molecular ecologist's guide to expressed sequence tags

Genomics and bioinformatics have great potential to help address numerous topics in ecology and evolution. Expressed sequence tags (ESTs) can bridge genomics and molecular ecology because they can provide a means of accessing the gene space of almost any organism. We review how ESTs have been used in molecular ecology research in the last several years by providing sequence data for the design of molecular markers, genome-wide studies of gene expression and selection, the identification of candidate genes underlying adaptation, and the basis for studies of gene family and genome evolution. Given the tremendous recent advances in inexpensive sequencing technologies, we predict that molecular ecologists will increasingly be developing and using EST collections in the years to come. With this in mind, we close our review by discussing aspects of EST resource development of particular relevance for molecular ecologists.

Comparative genome mapping among Picea glauca, P. mariana x P. rubens and P. abies, and correspondence with other Pinaceae

A composite linkage map was constructed from four individual maps for the conifer Picea glauca (Moench) Voss, from anonymous and gene-specific markfers (714 AFLPs, 38 SSRs, and 53 ESTPs). A total of 12 linkage groups were delineated with an average marker density of 2.7 cM. Macro-synteny and macro-colinearity comparisons with two other composite linkage maps developed for the species complex P. mariana (Mill.) B.S.P. x P. rubens Sarg., and for P. abies (L.) Karst. revealed an identical number of linkage groups and a remarkable conservation of the gene content and gene order of linkage groups over the million years since the split between these taxa. Identical gene order among taxa was observed for 10 of the 12 assembled composite linkage groups. The discovery of one breakdown in synteny between P. glauca and the other two taxa indicated the occurrence of an inter-chromosomal rearrangement involving an insertional translocation. Analysis of marker colinearity also revealed a putative segmental duplication. The combined information from these three Picea genomes validated and improved large-scale genome comparisons at the inter-generic level in the family Pinaceae by allowing for the identification of 11 homoeologous linkage groups between Picea and Pinus, and nine such groups between Picea and Pseudotsuga menziesii. Notably, the analysis of synteny among the three genera revealed a putative case of chromosomal fission and an inter-chromosomal rearrangement in the genome of P. menziesii. Both of these changes are inter-connected, indicating much instability in this part of the P. menziesii genome. Overall, the macro-structure of the Pinaceae genome was well conserved, which is notable given the Cretaceous origin of its main lineages.

A composite linkage map from two crosses for the species complex Picea mariana x Picea rubens and analysis of synteny with other Pinaceae

Four individual linkage maps were constructed from two crosses for the species complex Picea mariana (Mill.) B.S.P. x Picea rubens Sarg in order to integrate their information into a composite map and to compare with other Pinaceae. For all individual linkage maps, 12 major linkage groups were recovered with 306 markers per map on average. Before building the composite linkage map, the common male parent between the two crosses made it possible to construct a reference linkage map to validate the relative position of homologous markers. The final composite map had a length of 2,319 cM (Haldane) and contained a total of 1,124 positioned markers, including 1,014 AFLPs, 3 RAPDs, 53 SSRs, and 54 ESTPs, assembled into 12 major linkage groups. Marker density of the composite map was statistically homogenous and was much higher (one marker every 2.1 cM) than that of the individual linkage maps (one marker every 5.7 to 7.1 cM). Synteny was well conserved between individual, reference, and composite linkage maps and 94% of homologous markers were colinear between the reference and composite maps. The combined information from the two crosses increased by about 24% the number of anchor markers compared to the information from any single cross. With a total number of 107 anchor markers (SSRs and ESTPs), the composite linkage map is a useful starting point for large-scale genome comparisons at the intergeneric level in the Pinaceae. Comparisons of this map with those in Pinus and Pseudotsuga allowed the identification of one breakdown in synteny where one linkage group homologous to both Picea and Pinus corresponded to two linkage groups in Pseudotsuga. Implications for the evolution of the Pinaceae genome are discussed.

Proteome analysis of early somatic embryogenesis in Picea glauca

Forestry is a valuable natural resource for many countries. Rapid production of large quantities of genetically improved and uniform seedlings for restocking harvested lands is a key component of sustainable forest management programs. Clonal propagation through somatic embryogenesis has the potential to meet this need in conifers and can offer the added benefit of ensuring consistent seedling quality. Although in commercial use, mass production of conifers through somatic embryogenesis is relatively new and there are numerous biological unknowns regarding this complex developmental pathway. To aid in unravelling the embryo developmental process, two-dimensional electrophoresis was employed to quantitatively assess the expression levels of proteins across four stages of somatic embryo maturation in white spruce (0, 7, 21 and 35 days post abscisic acid treatment). Forty-eight differentially expressed proteins have been identified, which display a significant change in abundance as early as day 7 of embryo development. These proteins are involved in a variety of cellular processes, many of which have not previously been associated with embryo development. The identification of these proteins was greatly assisted by the availability of a substantial expressed sequence tag (EST) resource developed for white, sitka and interior spruce. The combined use of these spruce ESTs in conjunction with GenBank accessions for other plants improved the rate of protein identification from 38% to 62% when compared with GenBank alone using automated, high-throughput techniques. This underscores the utility of EST resources in a proteomic study of any species for which a genome sequence is unavailable.

Robust simple sequence repeat markers for spruce (Picea spp.) from expressed sequence tags

Traditionally, simple sequence repeat (SSR)markers have been developed from libraries of genomic DNA. However, the large, repetitive nature of conifer genomes makes development of robust, single-copy SSR markers from genomic DNA difficult. Expressed sequence tags (ESTs), or sequences of messenger RNA, offer the opportunity to exploit single, low-copy, conserved sequence motifs for SSR development. From a 20,275-unigene spruce EST set, we identified 44 candidate EST-SSR markers. Of these, 25 amplified and were polymorphic in white, Sitka, and black spruce; 20 amplified in all 3 spruce species tested; the remaining five amplified in all except one species. In addition, 101 previously described spruce SSRs (mostly developed from genomic DNA), were tested. Of these, 17 amplified across white,Sitka, and black spruce. The 25 EST-SSRs had approximately 9% less heterozygosity than the 17 genomic-derived SSRs (mean H=0.65 vs 0.72), but appeared to have less null alleles, as evidenced by much lower apparent inbreeding (mean F=0.046 vs 0.126). These robust SSRs are of particular use in comparative studies,and as the EST-SSRs are within the expressed portion of the genome, they are more likely to be associated with a particular gene of interest, improving their utility for quantitative trait loci mapping and allowing detection of selective sweeps at specific genes.

Inheritance and diversity of simple sequence repeat (SSR) microsatellite markers in various families of Picea abies

A large number of sequence-specific SSRs were screened by using electrophoresis on metaphore agarose gels with the bands visualized by ethidium bromide staining. Many SSRs appeared as codominant and many as dominant markers, with presence or absence of bands. A simple Mendelian inheritance pattern for most codominant and dominant SSR loci was found. For many codominant SSR markers, null alleles were detected. The proportion of dominant microsatellites detected in this study (close to 50 %) was much higher than that commonly reported in many other studies. A high proportion of dominant markers together with a high frequency of codominant markers with null alleles may represent two important limitations for the use of microsatellites in different studies. On the other hand, many polymorphic codominant SSR microsatellite markers were found to be highly repeatable, and can be used for population studies, seed certification, quality control of controlled crosses, paternity analysis, pollen contamination, and mapping of QTL in related families. In this paper, we report on the inheritance pattern and diversity of codominant and dominant SSR microsatellites in seven families of Picea abies sharing a common mother.

Single-copy, species-transferable microsatellite markers developed from loblolly pine ESTs

Microsatellites, or simple sequence repeats (SSRs), are usually regarded as the "markers of choice" in population genetics research because they exhibit high variability. The development cost of these markers is usually high. In addition, microsatellite primers developed for one species often do not cross-amplify in related species, requiring separate development for each species. However, microsatellites found in expressed sequence tags (ESTs) might better cross-amplify as they reside in or near conserved coding DNA. In this study, we identified 14 Pinus taeda (loblolly pine) EST-SSRs from public EST databases and tested for their cross-species transferability to P. contorta ssp. latifolia, P. ponderosa, and P. sylvestris. As part of our development of a P. contorta microsatellite set, we also compared their transferability to that of 99 traditional microsatellite markers developed in P. taeda and tested on P. contorta ssp. latifolia. Compared to traditional microsatellites, EST-SSRs had higher transfer rates across pine species; however, the level of polymorphism of microsatellites derived from ESTs was lower. Sequence analyses revealed that the frequencies of insertions/deletions and base substitutions were lower in EST-SSRs than in other types of microsatellites, confirming that EST-SSRs are more conserved than traditional SSRs. Our results also provide a battery of 23 polymorphic, robust microsatellite primer pairs for lodgepole pine.

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.

Rapid microsatellite analysis using discontinuous polyacrylamide gel electrophoresis

A method for screening large numbers of samples for microsatellites using discontinuous, non-denaturing polyacrylamide gels and rapid fluorescent gel staining is described. Disc electrophoresis on slab gels provides high-resolution of PCR products. It is useful for collecting population data once microsatellite loci have been characterized.

Highly informative nature of inter simple sequence repeat (ISSR) sequences amplified using tri- and tetra-nucleotide primers from DNA of cauliflower (Brassica oleracea var. botrytis L.)

Inter simple sequence repeat (ISSR) sequences as molecular markers can lead to the detection of polymorphism and also be a new approach to the study of SSR distribution and frequency. In this study, ISSR amplification with nonanchored primer was performed in closely related cauliflower lines. Fourty-four different amplified fragments were sequenced. Sequences of PCR products are delimited by the expected motifs and number of repeats, which validates the ISSR nonanchored primer amplification technique. DNA and amino acids homology search between internal sequences and databases (i) show that the majority of the internal regions of ISSR had homologies with known sequences, mainly with genes coding for proteins implicated in DNA interaction or gene expression, which reflected the significance of amplified ISSR sequences and (ii) display long and numerous homologies with the Arabidopsis thaliana genome. ISSR amplifications revealed a high conservation of these sequences between Arabidopsis thaliana and Brassica oleracea var. botrytis. Thirty-four of the 44 ISSRs had one or several perfect or imperfect internal microsatellites. Such distribution indicates the presence in genomes of highly concentrated regions of SSR, or "SSR hot spots." Among the four nonanchored primers used in this study, trinucleotide repeats, and especially (CAA)5, were the most powerful primers for ISSR amplifications regarding the number of amplified bands, level of polymorphism, and their nature.

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.

Microsatellite isolation and characterization in sunflower (Helianthus annuus L.)

Development of microsatellite markers for sunflower (Helianthus annuus L.) was performed to estimate their frequency, nature (structure), levels of polymorphism, usefulness for genotype identification, and calculation of genetic relationships between inbred lines representing the species diversity. Isolation was performed from a small-insert genomic library followed by hybridization screening using oligonucleotide probes containing different nucleotide arrays. In this work, 503 unique microsatellite clones were sequenced and 271 PCR primer sequences bordering the microsatellite repeat were designed. For polymorphism assessment, 16 H. annuus germplasm accessions were checked and 170 of the primers tested were shown to be polymorphic for the selected lines. The polymorphic microsatellites produced an average of 3.5 alleles/locus and an average polymorphism information content (PIC) of 0.55. The most frequently found motifs within polymorphic simple-sequence repeats (SSRs) were: (GA)n, (GT)n, (AT)n, followed by trinucleotides (ATT)n, (TGG)n, and (ATC)n, and the tetranucleotide (CATA)n. Most of the 170 SSRs obtained showed important differences in the 16 reference inbred lines used for their characterization. In this work, 20 of the most informative SSRs destined to sunflower genotyping and legal fingerprinting purposes are fully described.

Diversity, origin, and distribution of retrotransposons (gypsy and copia) in conifers

We examined the diversity, evolution, and genomic organization of retroelements in a wide range of gymnosperms. In total, 165 fragments of the reverse transcriptase (RT) gene domain were sequenced from PCR products using newly designed primers for gypsy-like retrotransposons and well-known primers for copia-like retrotransposons; representatives of long interspersed nuclear element (LINE) retroposons were also found. Gypsy and copia-like retroelements are a major component of the gymnosperm genome, and in situ hybridization showed that individual element families were widespread across the chromosomes, consistent with dispersion and amplification via an RNA intermediate. Most of the retroelement families were widely distributed among the gymnosperms, including species with wide taxonomic separation from the Northern and Southern Hemispheres. When the gymnosperm sequences were analyzed together with retroelements from other species, the monophyletic origin of plant copia, gypsy, and LINE groups was well supported, with an additional clade including badnaviral and other, probably virus-related, plant sequences as well as animal and fungal gypsy elements. Plant retroelements showed high diversity within the phylogenetic trees of both copia and gypsy RT domains, with, for example, retroelement sequences from Arabidopsis thaliana being present in many supported groupings. No primary branches divided major taxonomic clades such as angiosperms, monocotyledons, gymnosperms, or conifers or (based on smaller samples) ferns, Gnetales, or Sphenopsida (Equisetum), suggesting that much of the existing diversity was present early in plant evolution, or perhaps that horizontal transfer of sequences has occurred. Within the phylogenetic trees for both gypsy and copia, two clearly monophyletic gymnosperm/conifer clades were revealed, providing evidence against recent horizontal transfer. The results put the evolution of the large and relatively conserved genome structure of gymnosperms into the context of the diversity of other groups of plants.