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

Population genetics and plant growth experiments as prerequisite for conservation measures of the rare European aquatic plant Luronium natans (Alismataceae)

Information provided by population genetic studies is often necessary to effectively protect endangered species. In general, such data is scarce for aquatic plants and this holds also for Luronium natans, an aquatic macrophyte endemic to northwestern and western Europe. It is threatened across its whole distribution range due to human influences, in particular due to eutrophication and intensive fish farming. In spite of habitat protection populations continue to decline and re-introductions are one possibility to prevent the species' extinction. Therefore, insights in genetic diversity and relatedness of source populations is warranted. Thus, we performed Amplified Fragment-Length Polymorphism (AFLP) on two large populations in Saxony, Germany (Großenhainer Pflege and Niederspree), complemented with numerous additional occurrences from Europe. In addition, we conducted experiments on plant growth to assess optimal conditions for ex-situ cultivation taking water temperature, water level and substrate into account. We revealed considerably high levels of genetic diversity within populations (Shannon Indices ranged from 0.367 to 0.416) implying that populations are not restricted to clonal growth only but reproduce also by open-pollinated flowers. Remarkably, the two geographically close Saxon populations were genetically distant to each other but subpopulations within a locality were completely intermingled. Concerning optimal cultivation conditions, longest roots were obtained at temperatures >14°C and saturated, but not submerging water levels. Thus, our findings advocate for a re-introduction scheme from nearby source populations and provide detailed information on successful ex-situ cultivation.

Identification and Validation of Marketing Weight-Related SNP Markers Using SLAF Sequencing in Male Yangzhou Geese

Growth performance is a complex economic trait for avian production. The swan goose (Anser cygnoides) has never been exploited genetically like chickens or other waterfowl species such as ducks. Traditional phenotypic selection is still the main method for genetic improvement of geese body weight. In this study, specific locus amplified fragment sequencing (SLAF-seq) with bulked segregant analysis (BSA) was conducted for discovering and genotyping single nucleotide polymorphisms (SNPs) associated with marketing weight trait in male geese. A total of 149,045 SNPs were obtained from 427,093 SLAF tags with an average sequencing depth of 44.97-fold and a Q30 value of 93.26%. After SNPs' filtering, a total of 12,917 SNPs were included in the study. The 31 highest significant SNPs-which had different allelic frequencies-were further validated by individual-based AS-PCR genotyping in two populations. The association between 10 novel SNPs and the marketing weight of male geese was confirmed. The 10 significant SNPs were involved in linear regression model analysis, which confirmed single-SNP associations and revealed three types of SNP networks for marketing weight. The 10 significant SNPs were located within or close to 10 novel genes, which were identified. The qPCR analysis showed significant difference between genotypes of each SNP in seven genes. Developed SLAF-seq and identified genes will enrich growth performance studies, promoting molecular breeding applications to boost the marketing weight of Chinese geese.

Genetic Distinctiveness but Low Diversity Characterizes Rear-Edge Thuja standishii (Gordon) Carr. (Cupressaceae) Populations in Southwest Japan

Rear-edge populations are of significant scientific interest because they can contain allelic variation not found in core-range populations. However, such populations can differ in their level of genetic diversity and divergence reflecting variation in life-history traits, demographic histories and human impacts. Using 13 EST-microsatellites, we investigated the genetic diversity and differentiation of rear-edge populations of the Japanese endemic conifer Thuja standishii (Gordon) Carr. in southwest Japan from the core-range in northeast Japan. Range-wide genetic differentiation was moderate (Fst = 0.087), with northeast populations weakly differentiated (Fst = 0.047), but harboring high genetic diversity (average population-level Ar = 4.76 and Ho = 0.59). In contrast, rear-edge populations were genetically diverged (Fst = 0.168), but contained few unique alleles with lower genetic diversity (Ar = 3.73, Ho = 0.49). The divergence between rear-edge populations exceeding levels observed in the core-range and results from ABC analysis and species distribution modelling suggest that these populations are most likely relicts of the Last Glacial Maximum. However, despite long term persistence, low effective population size, low migration between populations and genetic drift have worked to promote the genetic differentiation of southwest Japan populations of T. standishii without the accumulation of unique alleles.

Construction of an SNP-based high-density linkage map for flax (Linum usitatissimum L.) using specific length amplified fragment sequencing (SLAF-seq) technology

Flax is an important crop for oil and fiber, however, no high-density genetic maps have been reported for this species. Specific length amplified fragment sequencing (SLAF-seq) is a high-resolution strategy for large scale de novo discovery and genotyping of single nucleotide polymorphisms. In this study, SLAF-seq was employed to develop SNP markers in an F2 population to construct a high-density genetic map for flax. In total, 196.29 million paired-end reads were obtained. The average sequencing depth was 25.08 in male parent, 32.17 in the female parent, and 9.64 in each F2 progeny. In total, 389,288 polymorphic SLAFs were detected, from which 260,380 polymorphic SNPs were developed. After filtering, 4,638 SNPs were found suitable for genetic map construction. The final genetic map included 4,145 SNP markers on 15 linkage groups and was 2,632.94 cM in length, with an average distance of 0.64 cM between adjacent markers. To our knowledge, this map is the densest SNP-based genetic map for flax. The SNP markers and genetic map reported in here will serve as a foundation for the fine mapping of quantitative trait loci (QTLs), map-based gene cloning and marker assisted selection (MAS) for flax.

Multiple markers, niche modelling, and bioregions analyses to evaluate the genetic diversity of a plant species complex

Background

The classification of closely related plants is not straightforward. These morphologically similar taxa frequently maintain their inter-hybridization potential and share ancestral polymorphisms as a consequence of their recent divergence. Under the biological species concept, they may thus not be considered separate species. The Petunia integrifolia complex is especially interesting because, in addition to the features mentioned above, its taxa share a pollinator, and their geographical ranges show multiple overlaps. Here, we combined plastid genome sequences, nuclear microsatellites, AFLP markers, ecological niche modelling, and bioregions analysis to investigate the genetic variability between the different taxa of the P. integrifolia complex in a comprehensive sample covering the entire geographical range of the complex.

Results

Results from molecular markers did not fully align with the current taxonomic classification. Niche modelling and bioregions analyses revealed that taxa were associated with different ecological constraints, indicating that the habitat plays an important role in preserving species boundaries. For three taxa, our analyses showed a mostly conserved, non-overlapping geographical distribution over time. However, for two taxa, niche modelling found an overlapping distribution over time; these taxa were also associated with the same bioregions.

Conclusions

cpDNA markers were better able to discriminate between Petunia taxa than SSRs and AFLPs. Overall, our results suggest that the P. integrifolia complex represents a continuum of individuals from distant and historically isolated populations, which share some morphological traits, but are established in four different evolutionary lineages.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1084-y) contains supplementary material, which is available to authorized users.

Development and characterization of EST-SSR markers via transcriptome sequencing in Brainea insignis (Aspleniaceae s.l.)

PREMISE OF THE STUDY

Brainea insignis (Aspleniaceae) is an endangered tree fern in China whose wild populations have been seriously damaged due to overexploitation. Expressed sequence tag-simple sequence repeat (EST-SSR) primers were developed to investigate its genetic diversity and provide resources for future conservation studies.

METHODS AND RESULTS

We obtained 72,897 unigenes of B. insignis using transcriptome sequencing and detected 15,006 SSRs in 12,058 unigenes. Based on these results, we designed 100 EST-SSR primer pairs and successfully amplified 52 of them in six individuals; 27 demonstrated polymorphisms after amplification against 72 individuals across three populations. Allele numbers ranged from three to 10, and the observed and expected heterozygosities ranged from 0.105 to 1.000 and from 0.523 to 0.865, respectively, in the tested populations. Most of these primers could be successfully amplified in two other fern species (Blechnumorientale and Chieniopteris harlandii).

CONCLUSIONS

These selected EST-SSRs are valuable for genetic diversity and conservation studies in B. insignis and other related fern species.

Molecular Resources from Transcriptomes in the Brassicaceae Family

The rapidly falling costs and the increasing availability of large DNA sequence data sets facilitate the fast and affordable mining of large molecular markers data sets for comprehensive evolutionary studies. The Brassicaceae (mustards) are an important species-rich family in the plant kingdom with taxa distributed worldwide and a complex evolutionary history. We performed Simple Sequence Repeats (SSRs) mining using de novo assembled transcriptomes from 19 species across the Brassicaceae in order to study SSR evolution and provide comprehensive sets of molecular markers for genetic studies within the family. Moreover, we selected the genus Cochlearia to test the transferability and polymorphism of these markers among species. Additionally, we annotated Cochlearia pyrenaica transcriptome in order to identify the position of each of the mined SSRs. While we introduce a new set of tools that will further enable evolutionary studies across the Brassicaceae, we also discuss some broader aspects of SSR evolution. Overall, we developed 2012 ready-to-use SSR markers with their respective primers in 19 Brassicaceae species and a high quality annotated transcriptome for C. pyrenaica. As indicated by our transferability test with the genus Cochlearia these SSRs are transferable to species within the genus increasing exponentially the number of targeted species. Also, our polymorphism results showed substantial levels of variability for these markers. Finally, despite its complex evolutionary history, SSR evolution across the Brassicaceae family is highly conserved and we found no deviation from patterns reported in other Angiosperms.

High-density genetic map construction and QTLs identification for plant height in white jute (Corchorus capsularis L.) using specific locus amplified fragment (SLAF) sequencing

Background

Genetic mapping and quantitative trait locus (QTL) detection are powerful methodologies in plant improvement and breeding. White jute (Corchorus capsularis L.) is an important industrial raw material fiber crop because of its elite characteristics. However, construction of a high-density genetic map and identification of QTLs has been limited in white jute due to a lack of sufficient molecular markers. The specific locus amplified fragment sequencing (SLAF-seq) strategy combines locus-specific amplification and high-throughput sequencing to carry out de novo single nuclear polymorphism (SNP) discovery and large-scale genotyping. In this study, SLAF-seq was employed to obtain sufficient markers to construct a high-density genetic map for white jute. Moreover, with the development of abundant markers, genetic dissection of fiber yield traits such as plant height was also possible. Here, we present QTLs associated with plant height that were identified using our newly constructed genetic linkage groups.

Results

An F₈ population consisting of 100 lines was developed. In total, 69,446 high-quality SLAFs were detected of which 5,074 SLAFs were polymorphic; 913 polymorphic markers were used for the construction of a genetic map. The average coverage for each SLAF marker was 43-fold in the parents, and 9.8-fold in each F₈ individual. A linkage map was constructed that contained 913 SLAFs on 11 linkage groups (LGs) covering 1621.4 cM with an average density of 1.61 cM per locus. Among the 11 LGs, LG1 was the largest with 210 markers, a length of 406.34 cM, and an average distance of 1.93 cM between adjacent markers. LG11 was the smallest with only 25 markers, a length of 29.66 cM, and an average distance of 1.19 cM between adjacent markers. ‘SNP_only’ markers accounted for 85.54% and were the predominant markers on the map. QTL mapping based on the F₈ phenotypes detected 11 plant height QTLs including one major effect QTL across two cultivation locations, with each QTL accounting for 4.14–15.63% of the phenotypic variance.

Conclusions

To our knowledge, the linkage map constructed here is the densest one available to date for white jute. This analysis also identified the first QTL in white jute. The results will provide an important platform for gene/QTL mapping, sequence assembly, genome comparisons, and marker-assisted selection breeding for white jute.

Level of genetic differentiation affects relative performances of expressed sequence tag and genomic SSRs

Microsatellites, also called simple sequence repeats (SSRs), are markers of choice to estimate relevant parameters for conservation genetics, such as migration rates, effective population size and kinship. Cross-amplification of SSRs is the simplest way to obtain sets of markers, and highly conserved SSRs have recently been developed from expressed sequence tags (EST) to improve SSR cross-species utility. As EST-SSRs are located in coding regions, the higher stability of their flanking regions reduces the frequency of null alleles and improves cross-species amplification. However, EST-SSRs have generally less allelic variability than genomic SSRs, potentially leading to differences in estimates of population genetic parameters such as genetic differentiation. To assess the potential of EST-SSRs in studies of within-species genetic diversity, we compared the relative performance of EST- and genomic SSRs following a multispecies approach on passerine birds. We tested whether patterns and levels of genetic diversity within and between populations assessed from EST- and from genomic SSRs are congruent, and we investigated how the relative efficiency of EST- and genomic SSRs is influenced by levels of differentiation. EST- and genomic SSRs ensured comparable inferences of population genetic structure in cases of strong genetic differentiation, and genomic SSRs performed slightly better than EST-SSRs when differentiation is moderate. However and interestingly, EST-SSRs had a higher power to detect weak genetic structure compared to genomic SSRs. Our study attests that EST-SSRs may be valuable molecular markers for conservation genetic studies in taxa such as birds, where the development of genomic SSRs is impeded by their low frequency.

Comparison of the effectiveness of ISJ and SSR markers and detection of outlier loci in conservation genetics of Pulsatilla patens populations

BACKGROUND

Research into the protection of rare and endangered plant species involves genetic analyses to determine their genetic variation and genetic structure. Various categories of genetic markers are used for this purpose. Microsatellites, also known as simple sequence repeats (SSR), are the most popular category of markers in population genetics research. In most cases, microsatellites account for a large part of the noncoding DNA and exert a neutral effect on the genome. Neutrality is a desirable feature in evaluations of genetic differences between populations, but it does not support analyses of a population's ability to adapt to a given environment or its evolutionary potential. Despite the numerous advantages of microsatellites, non-neutral markers may supply important information in conservation genetics research. They are used to evaluate adaptation to specific environmental conditions and a population's adaptive potential. The aim of this study was to compare the level of genetic variation in Pulsatilla patens populations revealed by neutral SSR markers and putatively adaptive ISJ markers (intron-exon splice junction).

METHODS

The experiment was conducted on 14 Polish populations of P. patens and three P. patens populations from the nearby region of Vitebsk in Belarus. A total of 345 individuals were examined. Analyses were performed with the use of eight SSR primers specific to P. patens and three ISJ primers.

RESULTS

SSR markers revealed a higher level of genetic variation than ISJ markers (He = 0.609, He = 0.145, respectively). An analysis of molecular variance (AMOVA) revealed that, the overall genetic diversity between the analyzed populations defined by parameters FST and Φ PT for SSR (20%) and Φ PT for ISJ (21%) markers was similar. Analysis conducted in the Structure program divided analyzed populations into two groups (SSR loci) and three groups (ISJ markers). Mantel test revealed correlations between the geographic distance and genetic diversity of Polish populations of P. patens for ISJ markers, but not for SSR markers.

CONCLUSIONS

The results of the present study suggest that ISJ markers can complement the analyses based on SSRs. However, neutral and adaptive markers should not be alternatively applied. Neutral microsatellite markers cannot depict the full range of genetic variation in a population because they do not enable to analyze functional variation. Although ISJ markers are less polymorphic, they can contribute to the reliability of analyses based on SSRs.

Identification and differentiation of Ricinus communis L. using SSR markers

The castor-oil plant (Ricinus communis L.), a member of the spurge family (Euphorbiaceae), is a versatile industrial oil crop that is cultivated in many tropical and subtropical regions of the world. Castor oil is of continuing importance to the global specialty chemical industry because it is the only commercial source of a hydroxylated fatty acid. Castor also has tremendous future potential as an industrial oilseed crop because of its high seed oil content, unique fatty acid composition, potentially high oil yields and ability to be grown under drought and saline conditions. Knowledge of genetic variability is important for breeding programs to provide the basis for developing desirable genotypes. The aim of this study was to assess genetic diversity within the set of 60 ricin genotypes using 10 SSR primers. Ten SSR primers revealed a total of 67 alleles ranging from 4 to 9 alleles per locus with a mean value of 6.70 alleles per locus. The PIC values ranged from 0.719 to 0.860 with an average value of 0.813 and the DI value ranged from 0.745 to 0.862 with an average value of 0.821. Probability of identity (PI) was low ranged from 0.004 to 0.018 with an average of 0.008. A dendrogram was constructed from a genetic distance matrix based on profiles of the 10 SSR loci using the unweighted pair-group method with the arithmetic average (UPGMA). According to analysis, the collection of 60 diverse accessions of castor bean was clustered into six clusters. We could not distinguish 2 genotypes grouped in cluster 1, RM-96 and RM-98, which are genetically the closest. Knowledge on the genetic diversity of castor can be used to future breeding programs of castor.

QTL mapping of cucumber fruit flesh thickness by SLAF-seq

Cucumber is an agriculturally and economically important vegetable crop worldwide. Fruit flesh thickness is an important trait for cucumber and also a central determinant of yield, yet little is known about the underlying mechanism of this trait. In this study, bulked segregant analysis (BSA) combined with specific length amplified fragment sequencing (SLAF-seq) was applied to finely map the gene that underlies fruit flesh thickness in cucumber. A 0.19-Mb-long quantitative trait locus on chromosome 2 controlling fruit flesh thickness (QTL fft2.1) was identified and further confirmed by simple sequence repeat (SSR) marker-based classical QTL mapping in 138 F₂ individuals. Gene prediction of this 0.19-Mb region identified 20 genes. Quantitative RT-PCR revealed higher expression levels of Csa2 M058670.1 (SET domain protein-lysine methyltransferase) in D8 (thick fruit flesh parent) compared with that in XUE1 (thin fruit flesh parent) during fruit development. Sequence alignment analysis of Csa2M058670.1 from thick and thin fruit flesh cucumber lines revealed a 4-bp deletion mutation in the promoter region of this candidate gene, which may result in the loss of Csa2M058670.1 activation in thin fruit flesh lines. The data presented herein suggest that Csa2M058670.1 is a possible candidate gene for controlling flesh thickness in cucumber.

Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

BACKGROUND

Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest.

METHODS

Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea.

RESULTS

We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics.

CONCLUSIONS

The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.

An empirical review: Characteristics of plant microsatellite markers that confer higher levels of genetic variation

During microsatellite marker development, researchers must choose from a pool of possible primer pairs to further test in their species of interest. In many cases, the goal is maximizing detectable levels of genetic variation. To guide researchers and determine which markers are associated with higher levels of genetic variation, we conducted a literature review based on 6782 genomic microsatellite markers published from 1997-2012. We examined relationships between heterozygosity (H e or H o) or allele number (A) with the following marker characteristics: repeat type, motif length, motif region, repeat frequency, and microsatellite size. Variation across taxonomic groups was also analyzed. There were significant differences between imperfect and perfect repeat types in A and H e. Dinucleotide motifs exhibited significantly higher A, H e, and H o than most other motifs. Repeat frequency and motif region were positively correlated with A, H e, and H o, but correlations with microsatellite size were minimal. Higher taxonomic groups were disproportionately represented in the literature and showed little consistency. In conclusion, researchers should carefully consider marker characteristics so they can be tailored to the desired application. If researchers aim to target high genetic variation, dinucleotide motif lengths with large repeat frequencies may be best.

Genetic differentiation and hybrid identification using microsatellite markers in closely related wild species

Identifying the genetic basis of speciation is critical for understanding the evolutionary history of closely related wild species. Recently diverged species facilitate the study of speciation because many genetic and morphological characteristics are still shared by the organisms under study. The Petunia genus grows in South American grasslands and comprises both recently diverged wild species and commercial species. In this work, we analysed two closely related species: Petunia exserta, which has a narrow endemic range and grows exclusively in rocky shelters, and Petunia axillaris, which is widely distributed and comprises three allopatric subspecies. Petunia axillaris ssp. axillaris and P. exserta occur in sympatry, and putative hybrids between them have been identified. Here, we analysed 14 expressed sequence tag-simple sequence repeats (EST-SSRs) in 126 wild individuals and 13 putative morphological hybrids with the goals of identifying differentially encoded alleles to characterize their natural genetic diversity, establishing a genetic profile for each taxon and to verify the presence of hybridization signal. Overall, 143 alleles were identified and all taxa contained private alleles. Four major groups were identified in clustering analyses, which indicated that there are genetic distinctions among the groups. The markers evaluated here will be useful in evolutionary studies involving these species and may help categorize individuals by species, thus enabling the identification of hybrids between both their putative taxa. The individuals with intermediate morphology presented private alleles of their both putative parental species, although they showed a level of genetic mixing that was comparable with some of the individuals with typical P. exserta morphology. The EST-SSR markers scattered throughout the Petunia genome are very efficient tools for characterizing the genetic diversity in wild taxa of this genus and aid in identifying interspecific hybrids based on the presence of private alleles. These properties indicate that these markers will be helpful tools in evolutionary studies.

European flint landraces grown in situ reveal adaptive introgression from modern maize

We have investigated the role of selection in the determination of the detected levels of introgression from modern maize hybrid varieties into maize landraces still cultivated in situ in Italy. We exploited the availability of a historical collection of landraces undertaken before the introduction and widespread use of modern maize, to analyse genomic changes that have occurred in these maize landraces over 50 years of co-existence with hybrid varieties. We have combined a previously published SSR dataset (n=21) with an AFLP loci dataset (n=168) to provide higher resolution power and to obtain a more detailed picture. We show that selection pressures for adaptation have favoured new alleles introduced by migration from hybrids. This shows the potential for analysis of historical introgression even over this short period of 50 years, for an understanding of the evolution of the genome and for the identification of its functionally important regions. Moreover, this demonstrates that landraces grown in situ represent almost unique populations for use for such studies when the focus is on the domesticated plant. This is due to their adaptation, which has arisen from their dynamic evolution under a continuously changing agro-ecological environment, and their capture of new alleles from hybridisation. We have also identified loci for which selection has inhibited introgression from modern germplasm and has enhanced the distinction between landraces and modern maize. These loci indicate that selection acted in the past, during the formation of the flint and dent gene pools. In particular, the locus showing the strongest signals of selection is a Misfit transposable element. Finally, molecular characterisation of the same samples with two different molecular markers has allowed us to compare their performances. Although the genetic-diversity and population-structure analyses provide the same global qualitative pattern, which thus provides the same inferences, there are differences related to their natures and characteristics.

Development of Cymbidium ensifolium genic-SSR markers and their utility in genetic diversity and population structure analysis in cymbidiums

Background

Cymbidium is a genus of 68 species in the orchid family, with extremely high ornamental value. Marker-assisted selection has proven to be an effective strategy in accelerating plant breeding for many plant species. Analysis of cymbidiums genetic background by molecular markers can be of great value in assisting parental selection and breeding strategy design, however, in plants such as cymbidiums limited genomic resources exist. In order to obtain efficient markers, we deep sequenced the C. ensifolium transcriptome to identify simple sequence repeats derived from gene regions (genic-SSR).

Result

The 7,936 genic-SSR markers were identified. A total of 80 genic-SSRs were selected, and primers were designed according to their flanking sequences. Of the 80 genic-SSR primer sets, 62 were amplified in C. ensifolium successfully, and 55 showed polymorphism when cross-tested among 9 Cymbidium species comprising 59 accessions. Unigenes containing the 62 genic-SSRs were searched against Non-redundant (Nr), Gene Ontology database (GO), eukaryotic orthologous groups (KOGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The search resulted in 53 matching Nr sequences, of which 39 had GO terms, 18 were assigned to KOGs, and 15 were annotated with KEGG. Genetic diversity and population structure were analyzed based on 55 polymorphic genic-SSR data among 59 accessions. The genetic distance averaged 0.3911, ranging from 0.016 to 0.618. The polymorphic index content (PIC) of 55 polymorphic markers averaged 0.407, ranging from 0.033 to 0.863. A model-based clustering analysis revealed that five genetic groups existed in the collection. Accessions from the same species were typically grouped together; however, C. goeringii accessions did not always form a separate cluster, suggesting that C. goeringii accessions were polyphyletic.

Conclusion

The genic-SSR identified in this study constitute a set of markers that can be applied across multiple Cymbidium species and used for the evaluation of genetic relationships as well as qualitative and quantitative trait mapping studies. Genic-SSR’s coupled with the functional annotations provided by the unigenes will aid in mapping candidate genes of specific function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0124-5) contains supplementary material, which is available to authorized users.

Developing conversed microsatellite markers and their implications in evolutionary analysis of the Bemisia tabaci complex

The study of population genetics among the Bemisia tabaci complex is limited due to the lack of conserved molecular markers. In this study, 358, 433 and 322 new polynucleotide microsatellites are separately identified from the transcriptome sequences of three cryptic species of the B. tabaci complex. The cross species transferability of 57 microsatellites was then experimentally validated. The results indicate that these markers are conserved and have high inter-taxon transferability. Thirteen markers were employed to assess the genetic relationships among six cryptic species of the B. tabaci complex. To our surprise, the inferred phylogeny was consistent with that of mitochondrial COI sequences, indicating that microsatellites have the potential to distinguish species of the B. tabaci complex. Our results demonstrate that development of microsatellites from transcriptome data is a fast and cost-effective approach. These markers can be used to analyze the population genetics and evolutionary patterns of the B. tabaci complex.

A high-density genetic map of cucumber derived from Specific Length Amplified Fragment sequencing (SLAF-seq)

High-density genetic map provides an essential framework for accurate and efficient genome assembly and QTL fine mapping. Construction of high-density genetic maps appears more feasible since the advent of next-generation sequencing (NGS), which eases SNP discovery and high-throughput genotyping of large population. In this research, a high-density genetic map of cucumber (Cucumis sativus L.) was successfully constructed across an F2 population by a recently developed Specific Length Amplified Fragment sequencing (SLAF-seq) method. In total, 18.69 GB of data containing 93,460,000 paired-end reads were obtained after preprocessing. The average sequencing depth was 44.92 in the D8 (female parent), 42.16 in the Jin5-508 (male parent), and 5.01 in each progeny. 79,092 high-quality SLAFs were detected, of which 6784 SLAFs were polymorphic, and 1892 of the polymorphic markers met the requirements for constructing genetic map. The genetic map spanned 845.87 cm with an average genetic distance of 0.45 cm. It is a reliable linkage map for fine mapping and molecular breeding of cucumber for its high marker density and well-ordered markers.

High-utility conserved avian microsatellite markers enable parentage and population studies across a wide range of species

BACKGROUND

Microsatellites are widely used for many genetic studies. In contrast to single nucleotide polymorphism (SNP) and genotyping-by-sequencing methods, they are readily typed in samples of low DNA quality/concentration (e.g. museum/non-invasive samples), and enable the quick, cheap identification of species, hybrids, clones and ploidy. Microsatellites also have the highest cross-species utility of all types of markers used for genotyping, but, despite this, when isolated from a single species, only a relatively small proportion will be of utility. Marker development of any type requires skill and time. The availability of sufficient "off-the-shelf" markers that are suitable for genotyping a wide range of species would not only save resources but also uniquely enable new comparisons of diversity among taxa at the same set of loci. No other marker types are capable of enabling this. We therefore developed a set of avian microsatellite markers with enhanced cross-species utility.

RESULTS

We selected highly-conserved sequences with a high number of repeat units in both of two genetically distant species. Twenty-four primer sets were designed from homologous sequences that possessed at least eight repeat units in both the zebra finch (Taeniopygia guttata) and chicken (Gallus gallus). Each primer sequence was a complete match to zebra finch and, after accounting for degenerate bases, at least 86% similar to chicken. We assessed primer-set utility by genotyping individuals belonging to eight passerine and four non-passerine species. The majority of the new Conserved Avian Microsatellite (CAM) markers amplified in all 12 species tested (on average, 94% in passerines and 95% in non-passerines). This new marker set is of especially high utility in passerines, with a mean 68% of loci polymorphic per species, compared with 42% in non-passerine species.

CONCLUSIONS

When combined with previously described conserved loci, this new set of conserved markers will not only reduce the necessity and expense of microsatellite isolation for a wide range of genetic studies, including avian parentage and population analyses, but will also now enable comparisons of genetic diversity among different species (and populations) at the same set of loci, with no or reduced bias. Finally, the approach used here can be applied to other taxa in which appropriate genome sequences are available.

Characterization of comparative genome-derived simple sequence repeats for acanthopterygian fishes

Simple sequence repeats (SSRs) have become one of the most popular molecular markers for population genetic studies. The application of SSR markers has often been limited to source species because SSR loci are too labile to be maintained in even closely related species. However, a few extremely conserved SSR loci have been reported. Here, we tested for the presence of conserved SSR loci in acanthopterygian fishes, which include over 14 000 species, by comparing the genome sequences of four acanthopterygian fishes. We also examined the comparative genome-derived SSRs (CG-SSRs) for their transferability across acanthopterygian fishes and their applicability to population genetic analysis. Forty-six SSR loci with conserved flanking regions were detected and examined for their transferability among seven nonacanthopterygian and 27 acanthopterygian fishes. The PCR amplification success rate in nonacanthopterygian fishes was low, ranging from 2.2% to 21.7%, except for Lophius litulon (Lophiiformes; 80.4%). Conversely, the rate in most acanthopterygian fishes exceeded 70.0%. Sequencing of these 46 loci revealed the presence of SSRs suitable for scoring while fragment analysis of 20 loci revealed polymorphisms in most of the acanthopterygian fishes. Population genetic analysis of Cottus pollux (Scorpaeniformes) and Sphaeramia orbicularis (Perciformes) using CG-SSRs showed that these populations did not deviate from linkage equilibrium or Hardy-Weinberg equilibrium. Furthermore, almost no loci showed evidence of null alleles, suggesting that CG-SSRs have strong resolving power for population genetic analysis. Our findings will facilitate the use of these markers in species in which markers remain to be identified.

Microsatellites for the marsh fritillary butterfly: de novo transcriptome sequencing, and a comparison with amplified fragment length polymorphism (AFLP) markers

BACKGROUND

Until recently the isolation of microsatellite markers from Lepidoptera has proved troublesome, expensive and time-consuming. Following on from a previous study of Edith's checkerspot butterfly, Euphydryas editha, we developed novel microsatellite markers for the vulnerable marsh fritillary butterfly, E. aurinia. Our goal was to optimize the process in order to reduce both time and cost relative to prevailing techniques. This was accomplished by using a combination of previously developed techniques: in silico mining of a de novo assembled transcriptome sequence, and genotyping the microsatellites found there using an economic method of fluorescently labelling primers.

PRINCIPAL FINDINGS

In total, we screened nine polymorphic microsatellite markers, two of which were previously published, and seven that were isolated de novo. These markers were able to amplify across geographically isolated populations throughout Continental Europe and the UK. Significant deviations from Hardy-Weinberg equilibrium were evident in some populations, most likely due to the presence of null alleles. However, we used an F(st) outlier approach to show that these markers are likely selectively neutral. Furthermore, using a set of 128 individuals from 11 populations, we demonstrate consistency in population differentiation estimates with previously developed amplified fragment length polymorphism (AFLP) markers (r = 0.68, p<0.001).

SIGNIFICANCE

Rapid development of microsatellite markers for difficult taxa such as Lepidoptera, and concordant results with other putatively neutral molecular markers, demonstrate the potential of de novo transcriptional sequencing for future studies of population structure and gene flow that are desperately needed for declining species across fragmented landscapes.

Genetic evidence for the uncoupling of local aquaculture activities and a population of an invasive species--a case study of Pacific oysters (Crassostrea gigas)

Human-mediated introduction of nonnative species into coastal areas via aquaculture is one of the main pathways that can lead to biological invasions. To develop strategies to counteract invasions, it is critical to determine whether populations establishing in the wild are self-sustaining or based on repeated introductions. Invasions by the Pacific oyster (Crassostrea gigas) have been associated with the growing oyster aquaculture industry worldwide. In this study, temporal genetic variability of farmed and wild oysters from the largest enclosed bay in Ireland was assessed to reconstruct the recent biological history of the feral populations using 7 anonymous microsatellites and 7 microsatellites linked to expressed sequence tags (ESTs). There was no evidence of EST-linked markers showing footprints of selection. Allelic richness was higher in feral than in aquaculture samples (P = 0.003, paired t-test). Significant deviations from Hardy-Weinberg equilibrium due to heterozygote deficiencies were detected for almost all loci and samples, most likely explained by the presence of null alleles. Relatively high genetic differentiation was found between aquaculture and feral oysters (largest pairwise multilocus F(ST) 0.074, P < 0.01) and between year classes of oysters from aquaculture (largest pairwise multilocus F(ST) 0.073, P < 0.01), which was also confirmed by the strong separation of aquaculture and wild samples using Bayesian clustering approaches. A 10-fold higher effective population size (N(e)) and a high number of private alleles in wild oysters suggest an established self-sustaining feral population. The wild oyster population studied appears demographically independent from the current aquaculture activities in the estuary and alternative scenarios of introduction pathways are discussed.

Genetic variation and structure of the endangered Lady Fern Athyrium viridescentipes based on ubiquitous genotyping

To clarify the genetic status and provide effective information for the conservation of Athyrium viridescentipes, a critically endangered fern species with only 103 individuals remaining in the wild, we conducted ubiquitous genotyping to determine the genotypes of all remnant individuals of the target species. We analyzed the genetic variation of the 103 known individuals in four populations by using 13 microsatellite loci. The genotypes of single spores from a sporophytic individual were also determined in order to reveal the breeding system of this species. The level of allelic variation in A. viridescentipes was significantly lower than that of closely related Athyrium species. The genetic composition of the four populations was rather similar. Sixty-nine individuals (67%) possessed an identical pattern in the allele combinations at 13 microsatellite loci. The mean pairwise F (ST) among four populations was 0.018. The segregated pattern of alleles, determined by single-spore genotyping, revealed that allelic recombination occurs through meiosis. The results indicate that this species contains a low level of genetic variation, has low population differentiation, and maintains populations by sexual reproduction. These findings could lead to more effective conservation programs, the selection of the most appropriate individuals for ex situ conservation efforts, and separate management of extant populations.

Construction of a high-density genetic map for grape using next generation restriction-site associated DNA sequencing

BACKGROUND

Genetic mapping and QTL detection are powerful methodologies in plant improvement and breeding. Construction of a high-density and high-quality genetic map would be of great benefit in the production of superior grapes to meet human demand. High throughput and low cost of the recently developed next generation sequencing (NGS) technology have resulted in its wide application in genome research. Sequencing restriction-site associated DNA (RAD) might be an efficient strategy to simplify genotyping. Combining NGS with RAD has proven to be powerful for single nucleotide polymorphism (SNP) marker development.

RESULTS

An F1 population of 100 individual plants was developed. In-silico digestion-site prediction was used to select an appropriate restriction enzyme for construction of a RAD sequencing library. Next generation RAD sequencing was applied to genotype the F1 population and its parents. Applying a cluster strategy for SNP modulation, a total of 1,814 high-quality SNP markers were developed: 1,121 of these were mapped to the female genetic map, 759 to the male map, and 1,646 to the integrated map. A comparison of the genetic maps to the published Vitis vinifera genome revealed both conservation and variations.

CONCLUSIONS

The applicability of next generation RAD sequencing for genotyping a grape F1 population was demonstrated, leading to the successful development of a genetic map with high density and quality using our designed SNP markers. Detailed analysis revealed that this newly developed genetic map can be used for a variety of genome investigations, such as QTL detection, sequence assembly and genome comparison.

Construction of a high-density genetic map for grape using next generation restriction-site associated DNA sequencing

BACKGROUND

Genetic mapping and QTL detection are powerful methodologies in plant improvement and breeding. Construction of a high-density and high-quality genetic map would be of great benefit in the production of superior grapes to meet human demand. High throughput and low cost of the recently developed next generation sequencing (NGS) technology have resulted in its wide application in genome research. Sequencing restriction-site associated DNA (RAD) might be an efficient strategy to simplify genotyping. Combining NGS with RAD has proven to be powerful for single nucleotide polymorphism (SNP) marker development.

RESULTS

An F1 population of 100 individual plants was developed. In-silico digestion-site prediction was used to select an appropriate restriction enzyme for construction of a RAD sequencing library. Next generation RAD sequencing was applied to genotype the F1 population and its parents. Applying a cluster strategy for SNP modulation, a total of 1,814 high-quality SNP markers were developed: 1,121 of these were mapped to the female genetic map, 759 to the male map, and 1,646 to the integrated map. A comparison of the genetic maps to the published Vitis vinifera genome revealed both conservation and variations.

CONCLUSIONS

The applicability of next generation RAD sequencing for genotyping a grape F1 population was demonstrated, leading to the successful development of a genetic map with high density and quality using our designed SNP markers. Detailed analysis revealed that this newly developed genetic map can be used for a variety of genome investigations, such as QTL detection, sequence assembly and genome comparison.

Population structure and linkage disequilibrium in Lupinus albus L. germplasm and its implication for association mapping

White lupin (Lupinus albus L.) has been around since 300 B.C. and is recognized for its ability to grow on poor soils and application as green manure in addition to seed harvest. The seed has very high levels of protein (33-47 %) and oil (6-13 %). It also has many secondary metabolites that are potentially of nutraceutical value to animals and humans. Despite such a great potential, lupins role in modern agriculture began only in the twentieth century. Although a large collection of Lupinus germplasm accessions is available worldwide, rarely have they been genetically characterized. Additionally, scarce genomic resources in terms of recombinant populations and genome information have been generated for L. albus. With the advancement in association mapping methods, the natural populations have the potential to replace the recombinant populations in gene mapping and marker-trait associations. Therefore, we studied the genetic similarity, population structure and marker-trait association in a USDA germplasm collection for their current and future application in this crop improvement. A total of 122 PI (Plant Inventory) lines were screened with 18 AFLP primer pairs that generated 2,277 fragments. A subset of 892 polymorphic markers with MAF >0.05 (minor allele frequency) were used for association mapping. The cluster analysis failed to group accessions on the basis of their passport information, and a weak structure and low linkage disequilibrium (LD) were observed indicating the usefulness of the collection for association mapping. Moreover, we were also able to identify two markers (a p value of 1.53 × 10(-4) and 2.3 × 10(-4)) that explained 22.69 and 20.5 % of seed weight variation determined using R (LR) (2) . The implications of lack of geographic clustering, population structure, low LD and the ability of AFLP to map seed weight trait using association mapping and the usefulness of the PI collections in breeding programs are discussed.

Comparative population genetic analysis of bocaccio rockfish Sebastes paucispinis using anonymous and gene-associated simple sequence repeat loci

Comparative population genetic analyses of traditional and emergent molecular markers aid in determining appropriate use of new technologies. The bocaccio rockfish Sebastes paucispinis is a high gene-flow marine species off the west coast of North America that experienced strong population decline over the past 3 decades. We used 18 anonymous and 13 gene-associated simple sequence repeat (SSR) loci (expressed sequence tag [EST]-SSRs) to characterize range-wide population structure with temporal replicates. No F(ST)-outliers were detected using the LOSITAN program, suggesting that neither balancing nor divergent selection affected the loci surveyed. Consistent hierarchical structuring of populations by geography or year class was not detected regardless of marker class. The EST-SSRs were less variable than the anonymous SSRs, but no correlation between F(ST) and variation or marker class was observed. General linear model analysis showed that low EST-SSR variation was attributable to low mean repeat number. Comparative genomic analysis with Gasterosteus aculeatus, Takifugu rubripes, and Oryzias latipes showed consistently lower repeat number in EST-SSRs than SSR loci that were not in ESTs. Purifying selection likely imposed functional constraints on EST-SSRs resulting in low repeat numbers that affected diversity estimates but did not affect the observed pattern of population structure.

Diverse spore rains and limited local exchange shape fern genetic diversity in a recently created habitat colonized by long-distance dispersal

BACKGROUND AND AIMS

Populations established by long-distance colonization are expected to show low levels of genetic variation per population, but strong genetic differentiation among populations. Whether isolated populations indeed show this genetic signature of isolation depends on the amount and diversity of diaspores arriving by long-distance dispersal, and time since colonization. For ferns, however, reliable estimates of long-distance dispersal rates remain largely unknown, and previous studies on fern population genetics often sampled older or non-isolated populations. Young populations in recent, disjunct habitats form a useful study system to improve our understanding of the genetic impact of long-distance dispersal.

METHODS

Microsatellite markers were used to analyse the amount and distribution of genetic diversity in young populations of four widespread calcicole ferns (Asplenium scolopendrium, diploid; Asplenium trichomanes subsp. quadrivalens, tetraploid; Polystichum setiferum, diploid; and Polystichum aculeatum, tetraploid), which are rare in The Netherlands but established multiple populations in a forest (the Kuinderbos) on recently reclaimed Dutch polder land following long-distance dispersal. Reference samples from populations throughout Europe were used to assess how much of the existing variation was already present in the Kuinderbos.

KEY RESULTS

A large part of the Dutch and European genetic diversity in all four species was already found in the Kuinderbos. This diversity was strongly partitioned among populations. Most populations showed low genetic variation and high inbreeding coefficients, and were assigned to single, unique gene pools in cluster analyses. Evidence for interpopulational gene flow was low, except for the most abundant species.

CONCLUSIONS

The results show that all four species, diploids as well as polyploids, were capable of frequent long-distance colonization via single-spore establishment. This indicates that even isolated habitats receive dense and diverse spore rains, including genotypes capable of self-fertilization. Limited gene flow may conserve the genetic signature of multiple long-distance colonization events for several decades.

Next-generation hybridization and introgression

Hybridization has a major role in evolution-from the introgression of important phenotypic traits between species, to the creation of new species through hybrid speciation. Molecular studies of hybridization aim to understand the class of hybrids and the frequency of introgression, detect the signature of ancient hybridization, and understand the behaviour of introgressed loci in their new genomic background. This often involves a large investment in the design and application of molecular markers, leading to a compromise between the depth and breadth of genomic data. New techniques designed to assay a large sub-section of the genome, in association with next-generation sequencing (NGS) technologies, will allow genome-wide hybridization and introgression studies in organisms with no prior sequence data. These detailed genotypic data will unite the breadth of sampling of loci characteristic of population genetics with the depth of sequence information associated with molecular phylogenetics. In this review, we assess the theoretical and methodological constraints that limit our understanding of natural hybridization, and promote the use of NGS for detecting hybridization and introgression between non-model organisms. We also make recommendations for the ways in which emerging techniques, such as pooled barcoded amplicon sequencing and restriction site-associated DNA tags, should be used to overcome current limitations, and enhance our understanding of this evolutionary significant process.

Wide variation in spatial genetic structure between natural populations of the European beech (Fagus sylvatica) and its implications for SGS comparability

Identification and quantification of spatial genetic structure (SGS) within populations remains a central element of understanding population structure at the local scale. Understanding such structure can inform on aspects of the species' biology, such as establishment patterns and gene dispersal distance, in addition to sampling design for genetic resource management and conservation. However, recent work has identified that variation in factors such as sampling methodology, population characteristics and marker system can all lead to significant variation in SGS estimates. Consequently, the extent to which estimates of SGS can be relied on to inform on the biology of a species or differentiate between experimental treatments is open to doubt. Following on from a recent report of unusually extensive SGS when assessed using amplified fragment length polymorphisms in the tree Fagus sylvatica, we explored whether this marker system led to similarly high estimates of SGS extent in other apparently similar populations of this species. In the three populations assessed, SGS extent was even stronger than this previously reported maximum, extending up to 360 m, an increase in up to 800% in comparison with the generally accepted maximum of 30-40 m based on the literature. Within this species, wide variation in SGS estimates exists, whether quantified as SGS intensity, extent or the Sp parameter. Consequently, we argue that greater standardization should be applied in sample design and SGS estimation and highlight five steps that can be taken to maximize the comparability between SGS estimates.

Cross-species transferability of SSR loci developed from transciptome sequencing in lodgepole pine

With the advent of next generation sequencing technologies, transcriptome level sequence collections are arising as prominent resources for the discovery of gene-based molecular markers. In a previous study more than 15,000 simple sequence repeats (SSRs) in expressed sequence tag (EST) sequences resulting from 454 pyrosequencing of Pinus contorta cDNA were identified. From these we developed PCR primers for approximately 4000 candidate SSRs. Here, we tested 184 of these SSRs for successful amplification across P. contorta and eight other pine species and examined patterns of polymorphism and allelic variability for a subset of these SSRs. Cross-species transferability was high, with high percentages of loci producing PCR products in all species tested. In addition, 50% of the loci we screened across panels of individuals from three of these species were polymorphic and allelically diverse. We examined levels of diversity in a subset of these SSRs by collecting genotypic data across several populations of Pinus ponderosa in northern Wyoming. Our results indicate the utility of mining pyrosequenced EST collections for gene-based SSRs and provide a source of molecular markers that should bolster evolutionary genetic investigations across the genus Pinus.

In silico mining and characterization of 12 EST-SSRs for the invasive slipper limpet Crepidula fornicata

In silico mining of an expressed sequence tags (ESTs) library was found to be efficient at isolating simple sequence repeats (SSRs) loci in the non-indigenous marine mollusc Crepidula fornicata. Twelve SSR loci were developed for routine genotyping. Cross-species amplification to 8 other Crepidula species showed that the 12 loci are highly specific for C. fornicata. Mendelian inheritance was shown for 11 of them (1 being monomorphic in the analyzed offspring array). The genetic diversity for 88 adults was found to be variable across the 12 loci (2-40 alleles, expected heterozygosity between 0.023 and 0.898) with a high overall exclusion probability of 0.99. The degree of genetic polymorphism found here is similar to that shown for 7 anonymous SSRs previously developed and here used on the same samples. This set of 12 specific loci is relevant to perform reliable population and relatedness analyses in Crepidula fornicata.

Mechanisms governing the responses to anthracnose pathogen in Juglans spp

Juglans nigra and Juglans regia are two highly economically important species for wood and fruit production that are susceptible to anthracnose caused by Gnomonia leptostyla. The identification of genotypes resistant to anthracnose could represent a valid alternative to agronomic and chemical management. In this study, we analyzed 72 walnut genotypes that showed a variety of resistance phenotypes in response to natural infection. According to the disease severity rating and microsatellite fingerprinting analysis, these genotypes were divided into three main groups: (40) J. nigra resistant, (1) J. nigra susceptible, and (31) J. regia susceptible. Data on leaf emergence rates and analysis of in vivo pathogenicity indicated that the incidence of anthracnose disease in the field might be partially conditioned by two key factors: the age and/or availability of susceptible leaves during the primary infection of fungus (avoidance by late flushing) and partial host resistance. NBS profiling approach, based on PCR amplification with an adapter primer for an adapter matching a restriction enzyme site and a degenerate primer targeting the conserved motifs present in the NBS domain of NBS-LRR genes, was applied. The results revealed the presence of a candidate marker that correlated to a reduction in anthracnose incidence in 72 walnut genotypes.

Fertility and precocity of Osmunda x intermedia offspring in culture

The feasibility of later-generation hybrid production in ferns has not been previously studied, although it is a significant factor in relation to reproductive isolation. Osmunda x intermedia, a hybrid between O. japonica and O. lancea, is semifertile and has moderate spore germination rates. Under the artificial conditions of this study, F2 and F3 offspring were formed. Some of the F2 offspring showed precocity, and some of the F3 offspring also showed precocity. This fertility suggests that introgressive hybridization might be ongoing in nature. This also indicates a currently unknown genetic control over the timing of fertile frond production in Osmunda.

Introgressive hybridization and the evolutionary history of the herring gull complex revealed by mitochondrial and nuclear DNA

Background

Based on extensive mitochondrial DNA (mtDNA) sequence data, we previously showed that the model of speciation among species of herring gull (Larus argentatus) complex was not that of a ring species, but most likely due more complex speciation scenario's. We also found that two species, herring gull and glaucous gull (L. hyperboreus) displayed an unexpected biphyletic distribution of their mtDNA haplotypes. It was evident that mtDNA sequence data alone were far from sufficient to obtain a more accurate and detailed insight into the demographic processes that underlie speciation of this complex, and that extensive autosomal genetic analysis was warranted.

Results

For this reason, the present study focuses on the reconstruction of the phylogeographic history of a limited number of gull species by means of a combined approach of mtDNA sequence data and 230 autosomal amplified fragment length polymorphism (AFLP) loci. At the species level, the mtDNA and AFLP genetic data were largely congruent. Not only for argentatus and hyperboreus, but also among a third species, great black-backed gull (L. marinus) we observed two distinct groups of mtDNA sequence haplotypes. Based on the AFLP data we were also able to detect distinct genetic subgroups among the various argentatus, hyperboreus, and marinus populations, supporting our initial hypothesis that complex demographic scenario's underlie speciation in the herring gull complex.

Conclusions

We present evidence that for each of these three biphyletic gull species, extensive mtDNA introgression could have taken place among the various geographically distinct subpopulations, or even among current species. Moreover, based on a large number of autosomal AFLP loci, we found evidence for distinct and complex demographic scenario's for each of the three species we studied. A more refined insight into the exact phylogeographic history within the herring gull complex is still impossible, and requires detailed autosomal sequence information, a topic of our future studies.

Population structure and linkage disequilibrium unravelled in tetraploid potato

Association mapping is considered to be an important alternative strategy for the identification of quantitative trait loci (QTL) as compared to traditional QTL mapping. A necessary prerequisite for association analysis to succeed is detailed information regarding hidden population structure and the extent of linkage disequilibrium. A collection of 430 tetraploid potato cultivars, comprising two association panels, has been analysed with 41 AFLP(®) and 53 SSR primer combinations yielding 3364 AFLP fragments and 653 microsatellite alleles, respectively. Polymorphism information content values and detected number of alleles for the SSRs studied illustrate that commercial potato germplasm seems to be equally diverse as Latin American landrace material. Genome-wide linkage disequilibrium (LD)-reported for the first time for tetraploid potato-was observed up to approximately 5 cM using r (2) higher than 0.1 as a criterion for significant LD. Within-group LD, however, stretched on average twice as far when compared to overall LD. A Bayesian approach, a distance-based hierarchical clustering approach as well as principal coordinate analysis were adopted to enquire into population structure. Groups differing in year of market release and market segment (starch, processing industry and fresh consumption) were repeatedly detected. The observation of LD up to 5 cM is promising because the required marker density is not likely to disable the possibilities for association mapping research in tetraploid potato. Population structure appeared to be weak, but strong enough to demand careful modelling of genetic relationships in subsequent marker-trait association analyses. There seems to be a good chance that linkage-based marker-trait associations can be identified at moderate marker densities.

Microsatellites reveal substantial among-population genetic differentiation and strong inbreeding in the relict fern Dryopteris aemula

BACKGROUND AND AIMS

A previous study detected no allozyme diversity in Iberian populations of the buckler-fern Dryopteris aemula. The use of a more sensitive marker, such as microsatellites, was thus needed to reveal the genetic diversity, breeding system and spatial genetic structure of this species in natural populations.

METHODS

Eight microsatellite loci for D. aemula were developed and their cross-amplification with other ferns was tested. Five polymorphic loci were used to characterize the amount and distribution of genetic diversity of D. aemula in three populations from the Iberian Peninsula and one population from the Azores.

KEY RESULTS

Most microsatellite markers developed were transferable to taxa close to D. aemula. Overall genetic variation was low (H(T) = 0.447), but was higher in the Azorean population than in the Iberian populations of this species. Among-population genetic differentiation was high (F(ST) = 0.520). All loci strongly departed from Hardy-Weinberg equilibrium. In the population where genetic structure was studied, no spatial autocorrelation was found in any distance class.

CONCLUSIONS

The higher genetic diversity observed in the Azorean population studied suggested a possible refugium in this region from which mainland Europe has been recolonized after the Pleistocene glaciations. High among-population genetic differentiation indicated restricted gene flow (i.e. lack of spore exchange) across the highly fragmented area occupied by D. aemula. The deviations from Hardy-Weinberg equilibrium reflected strong inbreeding in D. aemula, a trait rarely observed in homosporous ferns. The absence of spatial genetic structure indicated effective spore dispersal over short distances. Additionally, the cross-amplification of some D. aemula microsatellites makes them suitable for use in other Dryopteris taxa.

Rapid microsatellite isolation from a butterfly by de novo transcriptome sequencing: performance and a comparison with AFLP-derived distances

Background

The isolation of microsatellite markers remains laborious and expensive. For some taxa, such as Lepidoptera, development of microsatellite markers has been particularly difficult, as many markers appear to be located in repetitive DNA and have nearly identical flanking regions. We attempted to circumvent this problem by bioinformatic mining of microsatellite sequences from a de novo-sequenced transcriptome of a butterfly (Euphydryas editha).

Principal Findings

By searching the assembled sequence data for perfect microsatellite repeats we found 10 polymorphic loci. Although, like many expressed sequence tag-derived microsatellites, our markers show strong deviations from Hardy-Weinberg equilibrium in many populations, and, in some cases, a high incidence of null alleles, we show that they nonetheless provide measures of population differentiation consistent with those obtained by amplified fragment length polymorphism analysis. Estimates of pairwise population differentiation between 23 populations were concordant between microsatellite-derived data and AFLP analysis of the same samples (r = 0.71, p<0.00001, 425 individuals from 23 populations).

Significance

De novo transcriptional sequencing appears to be a rapid and cost-effective tool for developing microsatellite markers for difficult genomes.

Genetic variation in the allotetraploid Dryopteris corleyi (Dryopteridaceae) and its diploid parental species in the Iberian Peninsula

Studies on genetic diversity help us to unveil the evolutionary processes of species and populations and can explain several traits of diploid-polyploid complexes such as their distributions, their breeding systems, and the origin of polyploids. We examined the allozyme variation of Dryopteris aemula and D. oreades, diploid ferns with highly fragmented habitats, and the allotetraploid D. corleyi to (1) analyze the putative relationship between both diploids and the tetraploid, (2) compare the levels of genetic variation among species and determine their causes, and (3) assess the breeding system of these taxa. The allozymic pattern of D. corleyi confirms that it derived from D. aemula and D. oreades. The lack of genetic diversity in D. aemula, a species of lowland habitats, may be due to genetic drift associated with the contraction of populations in the last glaciation. By contrast, the alpine D. oreades had moderate intrapopulation genetic variation, which may derive from the expansion of populations during the last glaciation. In the latter species, low interpopulational variation suggested effective gene flow (spore exchange), and genotype frequencies in Hardy-Weinberg equilibrium indicated cross-fertilization of gametophytes. Evolutionary history appears to be an essential element in the interpretation of genetic variation of highly fragmented populations.

Effect of altitude on the genetic structure of an Alpine grass, Poa hiemata

BACKGROUND AND AIMS

The persistence of plants inhabiting restricted alpine areas under climate change will depend upon many factors including levels of genetic variation in adaptive traits, population structure, and breeding system.

METHODS

Using microsatellite markers, the genetic structure of populations of a relatively common alpine grass, Poa hiemata, is examined across three altitudinal gradients within the restricted Australian alpine zone where this species has previously been shown to exhibit local adaptation across a narrow altitudinal gradient.

KEY RESULTS

Genetic variation across six microsatellite markers revealed genetic structuring along altitudinal transects, and a reduction in genetic variation at high and low altitude extremes relative to sites central within transects. There was less genetic variation among transect sites compared with altitudinal gradients within transects, even though distances among transects were relatively larger. Central sites within transects were less differentiated than those at extremes.

CONCLUSIONS

These patterns suggest higher rates of gene flow among sites at similar altitudes than along transects, a process that could assist altitudinal adaptation. Patterns of spatial autocorrelation and isolation by distance changed with altitude and may reflect altered patterns of dispersal via pollen and/or seed. There was evidence for selfing and clonality in neighbouring plants. Levels of gene flow along transects were insufficient to prevent adaptive changes in morphological traits, given previously measured levels of selection.