Utility of EST-derived SSRs as population genetics markers in a beetle

Microsatellite, or simple sequence repeat (SSR), loci can be identified by mining expressed sequence tag (EST) databases, and where these are available, marker development time and expense can be decreased considerably over conventional strategies of probing the entire genome. However, it is unclear whether they provide information on population structure similar to that generated by anonymous genomic SSRs. We performed comparative population genetic analyses between EST-derived SSRs (EST-SSRs) and anonymous SSRs developed from genomic DNA for the same set of populations of the insect Diabrotica virgifera, a beetle in the family Chrysomelidae. Compared with noncoding, nontranscribed regions, EST-SSRs were generally less polymorphic but had reduced occurrence of null alleles and greater cross-species amplification. Neutrality tests suggested the loci were not under positive selection. Across all populations and all loci, the genomic and EST-SSRs performed similarly in estimating genetic diversity, F(IS), F(ST), population assignment and exclusion tests, and detection of distinct populations. These findings, therefore, indicate that the EST-SSRs examined can be used with confidence in future genetic studies of Diabrotica populations and suggest that EST libraries can be added as a valuable source of markers for population genetics studies in insects and other animals.     

SSR and morphological trait based population structure analysis of 130 diverse flax (Linum usitatissimum L.) accessions

A total of 130 flax accessions of diverse morphotypes and worldwide origin were assessed for genetic diversity and population structure using 11 morphological traits and microsatellite markers (15 gSSRs and 7 EST–SSRs). Analysis performed after classifying these accessions on the basis of plant height, branching pattern, seed size, Indian/foreign origin into six categories called sub-populations viz. fibre type exotic, fibre type indigenous, intermediate type exotic, intermediate type indigenous, linseed type exotic and linseed type indigenous. The study assessed different diversity indices, AMOVA, population structure and included a principal coordinate analysis based on different marker systems. The highest diversity was exhibited by gSSR markers (SI = 0.46; He = 0.31; P = 85.11). AMOVA based on all markers explained significant difference among fibre type, intermediate type and linseed type populations of flax. In terms of variation explained by different markers, EST-SSR markers (12%) better differentiated flax populations compared to morphological (9%) and gSSR (6%) markers at P = 0.01. The maximum Nei's unbiased genetic distance (D = 0.11) was observed between fibre type and linseed type exotic sub-populations based on EST-SSR markers. The combined structure analysis by using all markers grouped Indian fibre type accessions (63.4%) in a separate cluster along with the Indian intermediate type (48.7%), whereas Indian accessions (82.16%) of linseed type constituted an independent cluster. These findings were supported by the results of the principal coordinate analysis. Morphological markers employed in the study found complementary with microsatellite based markers in deciphering genetic diversity and population structure of the flax germplasm.     

Pros and cons of using genomic SSRs and EST-SSRs for resolving phylogeny of the genus Gossypium

The genus Gossypium is comprised of 50 diverse cotton species representing eight different genomes (A through G and K), however, phylogenetic relationship using various DNA marker types such as RAPD and SSRs was determined on limited number of cotton species. In this report, we have demonstrated the application of genomic SSRs (gSSRs) and EST-SSRs, and after combining both the data sets, for resolving the phylogenies of 36 cotton species including seven races. Out of the 100 primer pairs surveyed (50 for gSSRs and 50 for EST-SSRs), 75 produced scorable amplification products in all species. Out of these, 73 were found to be polymorphic and amplified 135 alleles ranging from 1 to 5 alleles per SSR marker (average 2.87 alleles per marker). The gSSRs amplified higher number of alleles (72) compared to the EST-SSRs (63). In total 22 highly informative SSRs with PIC values ≥0.5 were identified. Genomic SSRs containing di-while EST-SSRs containing tri-nucleotide repeats exhibited high polymorphism compared to the other nucleotide repeats containing gSSRs/EST-SSRs. Number of tandem repeats and polymorphism were positively correlated. Neither the type of chromosome nor the location of the SSRs showed association with the polymorphism. Gossypium herbaceum var. africanum (Watt) Hutch. ex and Ghose and Gossypium robinsonii F. Muell. were found the most genetically diverse, while among races of Gossypium hirsutum L. “yucatanense” and G. hirsutum “punctatum” were found genetically diverse. Of the three data sets, clustering analysis based on EST-SSRs and combined data sets, revealed parallel results reported in earlier studies. This study further confirmed that Gossypium darwinii Watt has close relationship with Gossypium barbadense L. Moreover, Gossypium raimondii Ulbr. and G. herbaceum/Gossypium arboreum L. are close living relatives of the ancestor allotetraploid species. Our studies suggest that for resolving phylogenetic relationship among the various plant species EST-SSRs could be a better choice. This information can be instrumental in transferring novel alleles or loci from the wild species into the cultivated cotton species which would set a stage for cultivating genetically diverse cultivars—a way to achieve sustainable cotton production in changing climate.     

Signatures of adaptation and genetic structure among the mainland populations of Pinus radiata (D. Don) inferred from SNP loci

Insights into the relative contributions of locus specific and genome-wide effects on population genetic diversity can be gained through separation of their resulting genetic signals. Here we explore patterns of adaptive and neutral genetic diversity in the disjunct natural populations of Pinus radiata (D. Don) from mainland California. A first-generation common garden of 447 individuals revealed significant differentiation of wood phenotypes among populations (P _ST), possibly reflecting local adaptation in response to environment. We subsequently screened all trees for genetic diversity at 149 candidate gene single nucleotide polymorphism (SNP) loci for signatures of adaptation. Ten loci were identified as being possible targets of diversifying selection following F _ST outlier tests. Multivariate canonical correlation performed on a data set of 444 individuals identified significant covariance between environment, adaptive phenotypes and outlier SNP diversity, lending support to the case for local adaptation suggested from F _ST and P _ST tests. Covariation among discrete sets of outlier SNPs and adaptive phenotypes (inferred from multivariate loadings) with environment are supported by existing studies of candidate gene function and genotype–phenotype association. Canonical analyses failed to detect significant correlations between environment and 139 non-outlier SNP loci, which were applied to estimate neutral patterns of genetic differentiation among populations (F _ST 4.3 %). Using this data set, significant hierarchical structure was detected, indicating three populations on the mainland. The hierarchical relationships based on neutral SNP markers (and SSR) were in contrast with those inferred from putatively adaptive loci, potentially highlighting the independent action of selection and demography in shaping genetic structure in this species.     

Establishment of a microsatellite panel covering the sugi (Cryptomeria japonica) genome,and its application for localization of a male-sterile gene (ms-2)

A microsatellite (simple sequence repeat; SSR) panel for Cryptomeria japonica was established, using both newly developed and previously reported markers, to construct a frame of linkage map and facilitate localization of important genes in this species. In this study, 32 new expressed sequence tag SSRs (EST-SSRs) and 12 new genomic SSRs (gSSRs) were developed. Their average polymorphism information content (PIC) values were 0.549 and 0.776, respectively. The markers were mapped onto a high-density linkage map. The SSR panel that was established to cover the genome consisted of 46 gSSRs and 47 EST-SSRs. The number of SSR markers in each linkage group, the average map distance between loci within a linkage group, and the average PIC values in each linkage group ranged from 6 to 13, 6.77 to 19.88 and 0.475 to 0.712, respectively. The utility of the SSR panel was tested by using it to localize a male-sterile gene, ms-2. The ms-2 locus was successfully localized on the linkage group 5 using 33 SSR markers (three SSRs per linkage group) which were selected from the SSR panel based on the existence of polymorphisms and the absence of null alleles in the mapping population for ms-2. A partial linkage map surrounding the ms-2 locus was then constructed using a further 57 single nucleotide polymorphisms and three SSRs, to facilitate future development of markers tightly linked to the ms-2 locus for use in marker-assisted selection. The SSR panel covering the C. japonica genome will allow researchers to localize important genes efficiently.     

Genetic Variation and Breeding Signature in Mass Selection Lines of the Pacific Oyster (Crassostrea gigas) Assessed by SNP Markers

In breeding industries, a challenging problem is how to keep genetic diversity over generations. To investigate genetic variation and identify breeding signatures in mass selected lines of Pacific oyster (Crassostrea gigas), three sixth-generation selected lines and four wild populations were assessed using 103 single nucleotide polymorphism (SNP) markers. The genetic diversity data indicated that the selected lines exhibited a significant reduction in the observed heterozygosity and observed number of alleles per locus compared with the wild populations (P≤0.05), indicating the selected lines tended to lose genetic diversity contrasted with the wild populations. The unweighted pair-group method with arithmetic mean (UPGMA) analysis showed that the wild populations and selected lines were not separated into two groups. Using four outlier tests, a total of 17 loci were found under selection at two levels. The global outlier detection suggested that 4 common outlier loci were subject to selection using both the hierarchical island model and Bayesian likelihood approaches. At regional level, 3 SNPs were detected as outlier using at least two outlier tests and one outlier SNP (CgSNP309) was overlapped in the two wild-selected population comparisons. The candidate outlier SNPs provide valuable resources for future association studies in C. gigas.     

Genetic structure of wild emmer wheat populations as reflected by transcribed versus anonymous SSR markers.

Simple sequence repeat (SSR) markers have become a major tool in population genetic analyses. The anonymous genomic SSRs (gSSRs) have been recently supplemented with expressed sequence tag (EST) derived SSRs (eSSRs), which represent the transcribed regions of the genome. In the present study, we used 8 populations of wild emmer wheat (Triticum turgidum subsp. dicoccoides) to compare the usefulness of the two types of SSR markers in assessing allelic diversity and population structure. gSSRs revealed significantly higher diversity than eSSRs in terms of average number of alleles (14.92 vs. 7.4, respectively), polymorphic information content (0.87 vs. 0.68, respectively), and gene diversity (He; 0.55 vs. 0.38, respectively). Despite the overall differences in the level of diversity, Mantel tests for correlations between eSSR and gSSR pairwise genetic distances were found to be significant for each population as well as for all accessions jointly (RM=0.54, p=0.01). Various genetic structure analyses (AMOVA, PCoA, STRUCTURE, unrooted UPGMA tree) revealed a better capacity of eSSRs to distinguish between populations, while gSSRs showed a higher proportion of intrapopulation (among accessions) diversity. We conclude that eSSR and gSSR markers should be employed in conjunction to obtain a high inter- and intra-specific (or inter- and intra-varietal) distinctness.     

Development and Characterization of Polymorphic EST-SSR and Genomic SSR Markers for Tibetan Annual Wild Barley

Tibetan annual wild barley is rich in genetic variation. This study was aimed at the exploitation of new SSRs for the genetic diversity and phylogenetic analysis of wild barley by data mining. We developed 49 novel EST-SSRs and confirmed 20 genomic SSRs for 80 Tibetan annual wild barley and 16 cultivated barley accessions. A total of 213 alleles were generated from 69 loci with an average of 3.14 alleles per locus. The trimeric repeats were the most abundant motifs (40.82%) among the EST-SSRs, while the majority of the genomic SSRs were di-nuleotide repeats. The polymorphic information content (PIC) ranged from 0.08 to 0.75 with a mean of 0.46. Besides this, the expected heterozygosity (He) ranged from 0.0854 to 0.7842 with an average of 0.5279. Overall, the polymorphism of genomic SSRs was higher than that of EST-SSRs. Furthermore, the number of alleles and the PIC of wild barley were both higher than that of cultivated barley, being 3.12 vs 2.59 and 0.44 vs 0.37. Indicating more polymorphism existed in the Tibetan wild barley than in cultivated barley. The 96 accessions were divided into eight subpopulations based on 69 SSR markers, and the cultivated genotypes can be clearly separated from wild barleys. A total of 47 SSR-containing EST unigenes showed significant similarities to the known genes. These EST-SSR markers have potential for application in germplasm appraisal, genetic diversity and population structure analysis, facilitating marker-assisted breeding and crop improvement in barley.     

Comparative analysis of population genetic structure in Athyrium distentifolium (Pteridophyta) using AFLPs and SSRs from anonymous and transcribed gene regions

To examine the performance and information content of different marker systems, comparative assessment of population genetic diversity was undertaken in nine populations of Athyrium distentifolium using nine genomic and 10 expressed sequence tag (EST) microsatellite (SSR) loci, and 265 amplified fragment length polymorphism (AFLP) loci from two primer combinations. In range-wide comparisons (European vs. North American populations), the EST-SSR loci showed more reliable amplification and produced more easily scorable bands than genomic simple sequence repeats (SSRs). Genomic SSRs showed significantly higher levels of allelic diversity than EST-SSRs, but there was a significant correlation in the rank order of population diversities revealed by both marker types. When AFLPs, genomic SSRs, and EST-SSRs are considered, comparisons of different population diversity metrics/markers revealed a mixture of significant and nonsignificant rank-order correlations. However, no hard incongruence was detected (in no pairwise comparison of populations did different marker systems or metrics detect opposingly significant different amounts of variation). Comparable population pairwise estimates of F(ST) were obtained for all marker types, but whilst absolute values for genomic and EST-SSRs were very similar (F(ST) = 0.355 and 0.342, respectively), differentiation was consistently higher for AFLPs in pairwise and global comparisons (global AFLP F(ST) = 0.496). The two AFLP primer combinations outperformed 18 SSR loci in assignment tests and discriminatory power in phenetic cluster analyses. The results from marker comparisons on A. distentifolium are discussed in the context of the few other studies on natural plant populations comparing microsatellite and AFLP variability.     

Development and analysis of EST-SSRs for flax (Linum usitatissimum L.)

A set of 146,611 expressed sequence tags (ESTs) were generated from 10 flax cDNA libraries. After assembly, a total of 11,166 contigs and 11,896 singletons were mined for the presence of putative simple sequence repeats (SSRs) and yielded 806 (3.5%) non-redundant sequences which contained 851 putative SSRs. This is equivalent to one EST-SSR per 16.5 kb of sequence. Trinucleotide motifs were the most abundant (76.9%), followed by dinucleotides (13.9%). Tetra-, penta- and hexanucleotide motifs represented <10% of the SSRs identified. A total of 83 SSR motifs were identified. Motif (TTC/GAA)n was the most abundant (10.2%) followed by (CTT/AAG)n (8.7%), (TCT/AGA)n (8.6%), (CT/AG)n (6.7%) and (TC/GA)n (5.3%). A total of 662 primer pairs were designed, of which 610 primer pairs yielded amplicons in a set of 23 flax accessions. Polymorphism between the accessions was found for 248 primer pairs which detected a total of 275 EST-SSR loci. Two to seven alleles were detected per marker. The polymorphism information content value for these markers ranged from 0.08 to 0.82 and averaged 0.35. The 635 alleles detected by the 275 polymorphic EST-SSRs were used to study the genetic relationship of 23 flax accessions. Four major clusters and two singletons were observed. Sub-clusters within the main clusters correlated with the pedigree relationships amongst accessions. The EST-SSRs developed herein represent the first large-scale development of SSR markers in flax. They have potential to be used for the development of genetic and physical maps, quantitative trait loci mapping, genetic diversity studies, association mapping and fingerprinting cultivars for example. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genetic structure of Quercus rubra L. and Quercus ellipsoidalis E. J. Hill populations at gene-based EST-SSR and nuclear SSR markers

Sympatric hybridizing oak species provide a model system for studying local adaptation. Disjunct populations of Quercus rubra L. and Quercus ellipsoidalis E. J. Hill at the northern edge of their distribution may harbor important reservoirs of adaptive genetic variation. Genic (expressed sequence tag- simple sequence repeat?=?EST-SSR) and non-genic nuclear microsatellite (nuclear SSR?=?nSSR) markers were used to estimate neutral and potentially adaptive genetic variation in these two supposedly interfertile oak species showing different adaptations to drought. Eleven populations of putative Q. rubra and Q. ellipsoidalis located in the Western Upper Peninsula of Michigan were characterized using seven EST-SSRs and eight nSSRs. Bayesian cluster analysis revealed two distinct groups corresponding to each species with evidence of low levels of potential introgression. A comparison of the genetic structure of adult trees and seedlings revealed no evidence for selection against hybrids. Overall, similar levels of genetic variation and differentiation between populations and species were found at both EST-SSRs and nSSRs indicating that most EST-SSRs chosen reflect neutral variation. Two loci, 3A05 (nSSR) and GOT021 (EST-SSR, putative histidine kinase 4-like), were identified as putative outlier loci between species showing largely reduced variation in Q. ellipsoidalis. Future analyses of an increased number of EST-SSRs located in functional genes will allow the identification of genes involved in the reproductive isolation between both species.     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.     

Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

De novo sequencing and development of EST-SSR for Gevuina avellana (Proteaceae), a nut tree from South America

?Premise of the study: The aim of this study was to develop molecular tools to investigate the genetic structure and diversity of natural populations of Gevuina avellana (Molina, Proteaceae), the Chilean hazelnut. Simple sequence repeats (SSRs) derived from expressed sequence tags (EST-SSRs) were developed, optimized, and characterized. ?Methods and Results: The microsatellite-containing sequences were selected from a cDNA library developed from the nut. The eight marker loci showed two to seven alleles in 60 unrelated trees, from two different natural populations. The observed heterozygosity (H(o)) was 0.44, ranging from 0.07 to 0.92 for different loci. When the multilocus genotype of the eight EST-SSRs was considered, all trees could be differentiated. ?Conclusions: The newly identified EST-SSRs in G. avellana and the multi-pooling strategy involving eight markers will facilitate future studies of clone identification, genetic structure, and diversity.     

Development and Characterization of Genomic and Gene-Based Microsatellite Markers in North American Red Oak Species

Oaks (Quercus: Fagaceae) are ecological and economic keystones of many forested ecosystems but effective genetic management strategies are hindered by high levels of phenotypic plasticity within species and frequent hybridization among them. These same features, however, make oak communities suited for the study of speciation, hybridization, and genetic adaptation. Efforts to develop new and to adapt existing genomic resources to less-studied members of this genus should not only improve oak conservation and management but also aid the study of fundamental evolutionary processes. Here, we present a suite of 27 highly polymorphic simple sequence repeat (SSR) markers tested in four North American red oak (Quercus section Lobatae) species: Q. rubra, Q. ellipsoidalis, Q. coccinea, and Q. velutina. Five markers are genomic SSRs (gSSRs) — four novel and one previously transferred from Q. petraea — and 22 are gene-based SSRs derived from Q. robur and Q. petraea expressed sequence tags (EST-SSRs). Overall, levels of polymorphism detected with these primer pairs were high, with gene diversity (H _e) averaging 0.66 across all loci in natural populations. In addition, we show that EST-SSR markers may have the potential to detect divergent selection at stress-resistance candidate genes among closely related oak species.     

SSR-based population structure, molecular diversity and linkage disequilibrium analysis of a collection of flax (Linum usitatissimum L.) varying for mucilage seed-coat content

Flax seed mucilage (SM) presents specific biological activities useful for the food and pharmaceutical industries. Understanding the population structure, genetic diversity and linkage disequilibrium (LD) of germplasm varying for mucilage content is pivotal for the identification of genes and quantitative trait loci underlying mucilage variation by association mapping (AM). In this study, 150 microsatellite loci were used to assess the population structure, genetic diversity and LD of a set of 60 flax cultivars/accessions capturing the breadth of SM variation in flax germplasm. STRUCTURE analysis and similarity-based methods revealed the presence of three populations reflecting mainly their geographic origins (South Asia, South America and North America), and the impact of germplasm exchange within and between North American flax breeding programs. Analysis of molecular variance showed that 78.32% of the genetic variation resided within populations and 21.68% among populations. The phi-statistic (??_st) value of 0.22 confirmed the presence of a strong population structure. A total of 408 alleles were detected, with the South American population capturing the highest overall diversity. However, the genetic diversity was narrow, as indicated by the small number of alleles per locus (2.72) and gene diversity (mean?=?0.34). LD was significant between 3.9% (r ²) and 36.2% (D??) of the loci pairs (FDR?<?0.05). The mean r ² and D?? were 0.26 and 0.53, respectively. The results suggest that the collection could be useful in AM studies aimed at the discovery of genes/alleles involved in SM; however a greater diversity may be required to improve the AM resolution.     

Assessment of genetic diversity and population structure of Chinese wild almond, Amygdalus nana, using EST- and genomic SSRs

Amygdalus nana L., commonly known as wild almond, is an endangered wild relative of cultivated almond, which has great potential in almond crop breeding. In this study, we used microsatellite (SSR) loci derived from both expressed sequence tag (EST) and anonymous genomic sequence to explore the genetic diversity and population structure of A. nana in Xinjiang of China. Seven natural populations were collected across the whole distribution of A. nana in China, including populations from both inside (four populations) and outside (three populations) the established protected areas. A total of 22 and 19 alleles were detected from the seven pairs of EST and genomic SSR loci, respectively. Generally, the genomic SSRs showed lower levels of variation than EST-SSRs, which may partially due to the higher cross-species transferability in EST-SSRs than in genomic SSRs. The population-level genetic diversity (A = 1.84, P = 50.00%, H_o = 0.3491, H_E = 0.2271) was lower than cultivated almond and several wild fruit species with similar breeding system. Most of the genetic variation (82.16%) was partitioned within populations. In particular, the population collected from Tacheng County (outside the protected areas) had the highest levels of genetic diversity and had significantly different genetic constitution from other populations.