Challenges in analysis and interpretation of microsatellite data for population genetic studies

Population genetic structure and connectivity of deep-sea stony corals (Order Scleractinia) in the New Zealand region: Implications for the conservation and management of vulnerable marine ecosystems

Deep-sea stony corals, which can be fragile, long-lived, late to mature and habitat-forming, are defined as vulnerable marine ecosystem indicator taxa. Under United Nations resolutions, these corals require protection from human disturbance such as fishing. To better understand the vulnerability of stony corals (Goniocorella dumosa, Madrepora oculata, Solenosmilia variabilis) to disturbance within the New Zealand region and to guide marine protected area design, genetic structure and connectivity were determined using microsatellite loci and DNA sequencing. Analyses compared population genetic differentiation between two biogeographic provinces, amongst three subregions (north-central-south) and amongst geomorphic features. Extensive population genetic differentiation was revealed by microsatellite variation, whilst DNA sequencing revealed very little differentiation. For G. dumosa, genetic differentiation existed amongst regions and geomorphic features, but not between provinces. For M. oculata, only a north-central-south regional structure was observed. For S. variabilis, genetic differentiation was observed between provinces, amongst regions and amongst geomorphic features. Populations on the Kermadec Ridge were genetically different from Chatham Rise populations for all three species. A significant isolation-by-depth pattern was observed for both marker types in G. dumosa and also in ITS of M. oculata. An isolation-by-distance pattern was revealed for microsatellite variation in S. variabilis. Medium to high levels of self-recruitment were detected in all geomorphic populations, and rates and routes of genetic connectivity were species-specific. These patterns of population genetic structure and connectivity at a range of spatial scales indicate that flexible spatial management approaches are required for the conservation of deep-sea corals around New Zealand.

Genetic Patterns of Myrceugenia correifolia, a Rare Species of Fog-Dependent Forests of Mediterranean Chile: Is It a Climatic Relict?

Rare species frequently occur in areas with microclimatic conditions that are atypical for their regions, but that were more common in the past, and that probably have operated as climatic refugia for a long time. Myrceugenia correifolia is a rare arboreal species that grows in deep canyons and hilltops of the Coast Range of north-central Chile between 30° and 35°S. In the northern edge of its distribution M. correifolia grows in small patches of fog-dependent forest surrounding by xeric vegetation. These forest formations are thought to be remnants of an ancient and continuous rainforest that according to some authors became fragmented during aridization of the Neogene (Neogene relict) and to others during warm-dry cycles of the Pleistocene (glacial relicts). Here we asked whether the northernmost populations of M. correifolia are Neogene relicts, glacial relicts, or the result of a recent northward colonization. To answer this question we examined genetic diversity and population divergence of M. correifolia using microsatellite markers, tested various competing population history scenarios with an approximate Bayesian computation (ABC) method, and complemented these data with ecological niche modeling (ENM). We detected three genetic clusters with a distinctive latitudinal pattern (north, center, and south) and high levels of differentiation (F_ST = 0.36). Demographic inference supported an admixture event 31 kya between two populations that diverged from an ancient population 139 kya. The admixture time coincides with the beginning of a period of wet conditions in north-central Chile that extended from 33 to 19 kya and was preceded by dry and cold conditions. These results suggest that increased precipitation during glacial periods triggered northward expansion of the range of M. correifolia, with subsequent admixture between populations that remained separated during interglacial periods. Accordingly, ENM models showed that suitable habitats for M. correifolia in north-central Chile were larger and less fragmented during the Last Glacial Maximum than at present, suggesting that northernmost populations of this species are glacial relicts.

Population structure of elephant foot yams (Amorphophallus paeoniifolius (Dennst.) Nicolson) in Asia

The corms and leaves of elephant foot yams (Amorphophallus paeoniifolius (Dennst.) Nicolson) are important foods in the local diet in many Asian regions. The crop has high productivity and wide agroecological adaptation and exhibits suitability for the agroforestry system. Although the plant is assumed to reproduce via panmixia, a comprehensive study on the genetic background across regions to enhance wider consumer palatability is still lacking. Here, ten informative microsatellites were analyzed in 29 populations across regions in India, Indonesia and Thailand to understand the genetic diversity, population structure and distribution to improve breeding and conservation programs. The genetic diversity was high among and within regions. Some populations exhibited excess heterozygosity and bottlenecking. Pairwise FST indicated very high genetic differentiation across regions (FST = 0.274), and the Asian population was unlikely to be panmictic. Phylogenetic tree construction grouped the populations according to country of origin with the exception of the Medan population from Indonesia. The current gene flow was apparent within the regions but was restricted among the regions. The present study revealed that Indonesia and Thailand populations could be alternative centers of the gene pool, together with India. Consequently, regional action should be incorporated in genetic conservation and breeding efforts to develop new varieties with global acceptance.

New insight into hybridization and unidirectional introgression between Ammodytes japonicus and Ammodytes heian (Trachiniformes, Ammodytidae)

Based on northern (NOL) and southern (SOL) mitochondrial lineages, recently, it proposed the new species Ammodytes heian and revived the species name Ammodytes japonicus to describe sand lances from the northwestern Pacific Ocean. This study used molecular methods to investigate genetic relationships between the two sand lance species in Korea and Japan. In total, 154 specimens were collected from four locations in Korea (Baengnyeongdo in the Yellow Sea, Tongyeong in the Korean Strait, and Jumunjin and Gijang in the East Sea), and 50 specimens were collected from a single location in Japan (Wakkanai in the Okhotsk Sea). Mitochondrial DNA analysis demonstrated that the individuals from Baengnyeongdo and Tongyeong all belonged to the SOL, whereas those from Gijang, Jumunjin, and Wakkanai included individuals from both the NOL and SOL (over 75% NOL). Population structure analyses were performed on the same individuals using seven microsatellite DNA markers. The population structure analysis based on 201 specimens identified two clusters (named as northern group and southern group), with the admixture proportion (q) of < 0.1 for the northern group in the Backyeongdo and Tongyeong sand lances and < 0.1 for the southern group in the Wakkanai sand lances. The high heterogeneity indicated that the former was probably A. japonicus and the latter probably A. heian. However, the admixture proportion in the Jumunjin and Gijang sand lances was 0.71–0.75 for the southern group, indicating that hybridization and unidirectional introgression from SOL to NOL occurs in southwestern margin of the East Sea. Our findings illustrate the speciation process based on different patterns of gene flow between Korean and Japanese sand lance, which is strongly influenced by both the paleo-climatic change and the contemporary local oceanic current pattern.

Genetic structure of Cercospora beticola populations on Beta vulgaris in New York and Hawaii

Cercospora leaf spot (CLS), caused by Cercospora beticola, is a major disease of Beta vulgaris worldwide. No sexual stage is known for C. beticola but in its asexual form it overwinters on infected plant debris as pseudostromata, and travels short distances by rain splash-dispersed conidiospores. Cercospora beticola infects a broad range of host species and may be seedborne. The relative contribution of these inoculum sources to CLS epidemics on table beet is not well understood. Pathogen isolates collected from table beet, Swiss chard and common lambsquarters in mixed-cropping farms and monoculture fields in New York and Hawaii, USA, were genotyped (n = 600) using 12 microsatellite markers. All isolates from CLS symptoms on lambsquarters were identified as C. chenopodii. Sympatric populations of C. beticola derived from Swiss chard and table beet were not genetically differentiated. Results suggested that local (within field) inoculum sources may be responsible for the initiation of CLS epidemics in mixed-cropping farms, whereas external sources of inoculum may be contributing to CLS epidemics in the monoculture fields in New York. New multiplex PCR assays were developed for mating-type determination for C. beticola. Implications of these findings for disease management are discussed.

New gSSR and EST-SSR markers reveal high genetic diversity in the invasive plant Ambrosia artemisiifolia L. and can be transferred to other invasive Ambrosia species

Ambrosia artemisiifolia L., (common ragweed), is an annual invasive and highly troublesome plant species originating from North America that has become widespread across Europe. New sets of genomic and expressed sequence tag (EST) based simple sequence repeats (SSRs) markers were developed in this species using three approaches. After validation, 13 genomic SSRs and 13 EST-SSRs were retained and used to characterize the genetic diversity and population genetic structure of Ambrosia artemisiifolia populations from the native (North America) and invasive (Europe) ranges of the species. Analysing the mating system based on maternal families did not reveal any departure from complete allogamy and excess homozygosity was mostly due the presence of null alleles. High genetic diversity and patterns of genetic structure in Europe suggest two main introduction events followed by secondary colonization events. Cross-species transferability of the newly developed markers to other invasive species of the Ambrosia genus was assessed. Sixty-five percent and 75% of markers, respectively, were transferable from A. artemisiifolia to Ambrosia psilostachya and Ambrosia tenuifolia. 40% were transferable to Ambrosia trifida, this latter species being seemingly more phylogenetically distantly related to A. artemisiifolia than the former two.

APMicroDB: A microsatellite database of Acyrthosiphon pisum

Pea aphids represent a complex genetic system that could be used for QTL analysis, genetic diversity and population genetics studies. Here, we described the development of first microsatellite repeat database of the pea aphid (APMicroDB), accessible at “http://deepaklab.com/aphidmicrodb”. We identified 3,40,233 SSRs using MIcroSAtellite (MISA) tool that was distributed in 14,067 (out of 23,924) scaffold of the pea aphid. We observed 89.53% simple repeats of which 73.41% were mono-nucleotide, followed by di-nucleotide repeats. This database stored information about the repeats kind, GC content, motif type (mono - hexa), genomic location etc. We have also incorporated the primer information derived from Primer3 software of the 250bp flanking region of the identified marker. Blast tool is also provided for searching the user query sequence for identified marker and their primers. This work has an immense use for scientific community working in the field of agricultural pest management, QTL mapping, and host-pathogen interaction analysis.

Tiger on the prowl: Invasion history and spatio-temporal genetic structure of the Asian tiger mosquito Aedes albopictus (Skuse 1894) in the Indo-Pacific

Background

Within the last century, increases in human movement and globalization of trade have facilitated the establishment of several highly invasive mosquito species in new geographic locations with concurrent major environmental, economic and health consequences. The Asian tiger mosquito, Aedes albopictus, is an extremely invasive and aggressive daytime-biting mosquito that is a major public health threat throughout its expanding range.

Methodology/Principal findings

We used 13 nuclear microsatellite loci (on 911 individuals) and mitochondrial COI sequences to gain a better understanding of the historical and contemporary movements of Ae. albopictus in the Indo-Pacific region and to characterize its population structure. Approximate Bayesian computation (ABC) was employed to test competing historical routes of invasion of Ae. albopictus within the Southeast (SE) Asian/Australasian region. Our ABC results show that Ae. albopictus was most likely introduced to New Guinea via mainland Southeast Asia, before colonizing the Solomon Islands via either Papua New Guinea or SE Asia. The analysis also supported that the recent incursion into northern Australia’s Torres Strait Islands was seeded chiefly from Indonesia. For the first time documented in this invasive species, we provide evidence of a recently colonized population (the Torres Strait Islands) that has undergone rapid temporal changes in its genetic makeup, which could be the result of genetic drift or represent a secondary invasion from an unknown source.

Conclusions/Significance

There appears to be high spatial genetic structure and high gene flow between some geographically distant populations. The species' genetic structure in the region tends to favour a dispersal pattern driven mostly by human movements. Importantly, this study provides a more widespread sampling distribution of the species’ native range, revealing more spatial population structure than previously shown. Additionally, we present the most probable invasion history of this species in the Australasian region using ABC analysis.

The Asian tiger mosquito, Aedes albopictus, is an aggressive mosquito that has expanded globally in the last century, chiefly due to the increase of human movements. It is a major public health concern due to its role in transmitting dengue, chikungunya and Zika viruses. Understanding how populations of Ae. albopictus are genetically related and how they have been introduced into new regions is important for controlling them and assessing their disease risk; few studies have explored this in the Indo-Pacific. In our study, we sampled a broader range of populations of Ae. albopictus in the Indo-Pacific to explore genetic patterns and to investigate the likely route of invasion through Australasia. We uncovered clear genetic groups throughout this region, but also found that some geographically distant populations are closely related, likely due to human-associated movements. We also found, that Ae. albopictus likely colonized New Guinea from mainland Southeast (SE) Asia before spreading to the Solomon Islands via either PNG or SE Asia. In contrast, the populations in Australia’s Torres Strait Islands were introduced from Indonesia. Interestingly, we found major genetic changes over time in some Torres Strait populations, less than a decade after its introduction, potentially reflective of a random reduction in genetic diversity (genetic drift) or a secondary invasion.

Genetic diversity of Atlantic Bluefin tuna in the Mediterranean Sea: insights from genome-wide SNPs and microsatellites

Background

Elucidating the patterns of the Atlantic Bluefin tuna [ABFT, Thunnus thynnus (Linnaeus, 1758)] population structure constitutes a challenging task of great importance. Most of the unique challenges stem from its biology, as well as the attributes of the marine realm in which it disperses. Accurate information is urgently needed for stock assessment, and the identification of critical features to the persistence and adaptation of populations in order to formulate and adopt effective strategies for ABFT conservation and management. Conclusions of a great number of ABFT genetic studies on the Mediterranean Sea stock structure are rather controversial and not yet conclusive. In this study, ABFT genomic diversity was investigated in the Mediterranean Sea, which is the most important area for the species’ reproduction.

Results

Analyzing genome-wide SNPs and microsatellites from ABFT samples collected throughout the Mediterranean Sea did not provide strong evidence of genetic structure, pointing towards the existence of a single panmictic unit. An alternative view would recognize a failure to reject the null hypothesis of a panmictic unit as an effect of the study’s sampling design, the type of markers used, and the effectiveness/suitability of analysis methods in respect to the species biological characteristics or any combination of the above.

Conclusions

Unravelling the drivers of ABFT population diversity would require the consideration of important aspects of the species spawning behavior for the determination of the appropriate sampling design. Novel approaches and methods of analysis that will bring together experts in genetics/-omics, ecology and oceanography are deemed necessary. Analyzing ABFT genetic data under the discipline of seascape genetics could provide the analysis framework under which major abiotic and biotic forces controlling ABFT recruitment could be identified, elucidating the complicated population dynamics of the species, while multiple and continuous fisheries monitoring should in all cases be considered as a prerequisite in order to achieve efficient and long-term ABFT conservation.

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.

Estimating genomic diversity and population differentiation - an empirical comparison of microsatellite and SNP variation in Arabidopsis halleri

BACKGROUND

Microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. However, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. Here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms (SNPs) to assess and quantify potential marker-specific biases and derive recommendations for future studies. Overall, we genotyped 180 Arabidopsis halleri individuals from nine populations using 20 microsatellite markers. Twelve of these markers were originally developed for Arabidopsis thaliana (cross-species markers) and eight for A. halleri (species-specific markers). We further characterized 2 million SNPs across the genome with a pooled whole-genome re-sequencing approach (Pool-Seq).

RESULTS

Our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity (SSR-H e) was not significantly correlated with genome-wide SNP diversity estimates (SNP-H e and θ Watterson) in A. halleri. Instead, microsatellite allelic richness (A r) was a better proxy for genome-wide SNP diversity. Estimates of genetic differentiation among populations (F ST) based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from SNPs. Possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. In contrast, genome-wide SNP data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide SNPs were used. Further, we inferred that a few thousand random SNPs are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation.

CONCLUSIONS

We recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based SNP data to draw conclusions on genome-wide diversity patterns. In species comparable to A. halleri, a few thousand SNPs are sufficient to achieve this goal.

High genetic diversity and connectivity in Colossoma macropomum in the Amazon basin revealed by microsatellite markers

Colossoma macropomum is the second largest scaled fish of the Amazon. It is economically important for commercial fisheries and for aquaculture, but few studies have examined the diversity and genetic structure of natural populations of this species. The aim of this study was to investigate the levels of genetic variability and connectivity that exist between three natural populations of C. macropomum from the Amazon basin. In total, 247 samples were collected from the municipalities of Tefé, Manaus, and Santarém. The populations were genotyped using a panel of 12 multiplex microsatellite markers. The genetic diversity found in these populations was high and similar to other populations described in the literature. These populations showed a pattern of high gene flow associated with the lack of a genetic structure pattern, indicating that the number of migrants per generation and recent migration rates are high. The values of the FST, RST, and exact test of differentiation were not significant for pairwise comparisons between populations. The Bayesian population clustering analysis indicated a single population. Thus, the data provide evidence for high genetic diversity and high gene flow among C. macropomum populations in the investigated region of the Amazon basin. This information is important for programs aiming at the conservation of natural populations.

Next Gen Pop Gen: implementing a high-throughput approach to population genetics in boarfish (Capros aper)

The recently developed approach for microsatellite genotyping by sequencing (GBS) using individual combinatorial barcoding was further improved and used to assess the genetic population structure of boarfish (Capros aper) across the species' range. Microsatellite loci were developed de novo and genotyped by next-generation sequencing. Genetic analyses of the samples indicated that boarfish can be subdivided into at least seven biological units (populations) across the species' range. Furthermore, the recent apparent increase in abundance in the northeast Atlantic is better explained by demographic changes within this area than by influx from southern or insular populations. This study clearly shows that the microsatellite GBS approach is a generic, cost-effective, rapid and powerful method suitable for full-scale population genetic studies-a crucial element for assessment, sustainable management and conservation of valuable biological resources.

Patterns of population structure at microsatellite and mitochondrial DNA markers in the franciscana dolphin (Pontoporia blainvillei)

The franciscana dolphin, Pontorporia blainvillei, is an endemic cetacean of the Atlantic coast of South America. Its coastal distribution and restricted movement patterns make this species vulnerable to anthropogenic factors, particularly to incidental bycatch. We used mitochondrial DNA control region sequences, 10 microsatellites, and sex data to investigate the population structure of the franciscana dolphin from a previously established management area, which includes the southern edge of its geographic range. F‐statistics and Bayesian cluster analyses revealed the existence of three genetically distinct populations. Based on the microsatellite loci, similar levels of genetic variability were found in the area; 13 private alleles were found in Monte Hermoso, but none in Claromecó. When considering the mitochondrial DNA control region sequences, lower levels of genetic diversity were found in Monte Hermoso, when compared to the other localities. Low levels of gene flow were found between most localities. Additionally, no evidence of isolation by distance nor sex‐biased dispersal was detected in the study area. In view of these results showing that populations from Necochea/Claromecó, Monte Hermoso, and Río Negro were found to be genetically distinct and the available genetic information for the species previously published, Argentina would comprise five distinct populations: Samborombón West/Samborombón South, Cabo San Antonio/Buenos Aires East, Necochea/Claromecó/Buenos Aires Southwest, Monte Hermoso, and Río Negro. In order to ensure the long‐term survival of the franciscana dolphin, management and conservation strategies should be developed considering each of these populations as different management units.

The importance of correctly identifying the process responsible for spatial genetic structure in Leopard: A response to McManus and Smuts (2016)

Microsatellite analyses suggest that spatial genetic structure among six leopard-sampling sites in southern Africa is the result of isolation by distance.

Assessing performance of single-sample molecular genetic methods to estimate effective population size: empirical evidence from the endangered Gochu Asturcelta pig breed

Estimating effective population size (N_e) using linkage disequilibrium (LD) information (N_e(LD)) has the operational advantage of using a single sample. However, N_e(LD) estimates assume discrete generations and its performance are constrained by demographic issues. However, such concerns have received little empirical attention so far. The pedigree of the endangered Gochu Asturcelta pig breed includes individuals classified into discrete filial generations and individuals with generations overlap. Up to 780 individuals were typed with a set of 17 microsatellites. Performance of N_e(LD) was compared with N_e estimates obtained using genealogical information, molecular coancestry (N_e(M)) and a temporal (two‐sample) method (N_e(JR)). Molecular‐based estimates of N_e exceeded those obtained using pedigree data. Estimates of N_e(LD) for filial generations F3 and F4 (17.0 and 17.3, respectively) were lower and steadier than those obtained using yearly or biannual samplings. N_e(LD) estimated for samples including generations overlap could only be compared with those obtained for the discrete filial generations when sampling span approached a generation interval and demographic correction for bias was applied. Single‐sample N_e(M) estimates were lower than their N_e(LD) counterparts. N_e(M) estimates are likely to partially reflect the number of founders rather than population size. In any case, estimates of LD and molecular coancestry tend to covary and, therefore, N_e(M) and N_e(LD) can hardly be considered independent. Demographically adjusted estimates of N_e(JR) and N_e(LD) took comparable values when: (1) the two samples used for the former were separated by one equivalent to discrete generations in the pedigree and (2) sampling span used for the latter approached a generation interval. Overall, the empirical evidence given in this study suggested that the advantage of using single‐sample methods to obtain molecular‐based estimates of N_e is not clear in operational terms. Estimates of N_e obtained using methods based in molecular information should be interpreted with caution.

A Galaxy-based bioinformatics pipeline for optimised, streamlined microsatellite development from Illumina next-generation sequencing data

Microsatellites are useful tools for ecologists and conservationist biologists, but are taxa-specific and traditionally expensive and time-consuming to develop. New methods using next-generation sequencing (NGS) have reduced these problems, but the plethora of software available for processing NGS data may cause confusion and difficulty for researchers new to the field of bioinformatics. We developed a bioinformatics pipeline for microsatellite development from Illumina paired-end sequences, which is packaged in the open-source bioinformatics tool Galaxy. This optimises and streamlines the design of a microsatellite panel and provides a user-friendly graphical user interface. The pipeline utilises existing programs along with our own novel program and wrappers to: quality-filter and trim reads (Trimmomatic); generate sequence quality reports (FastQC); identify potentially-amplifiable microsatellite loci (Pal_finder); design primers (Primer3); assemble pairs of reads to enhance marker amplification success rates (PANDAseq); and filter optimal loci (Pal_filter). The complete pipeline is freely available for use via a pre-configured Galaxy instance, accessible at https://palfinder.ls.manchester.ac.uk.

The repeatability of genome-wide mutation rate and spectrum estimates

Over the last decade, mutation studies have grown in popularity due to the affordability and accessibility of whole genome sequencing. As the number of species in which spontaneous mutation has been directly estimated approaches 20 across two domains of life, questions arise over the repeatability of results in such experiments. Five species were identified in which duplicate mutation studies have been performed. Across these studies the difference in estimated spontaneous mutation rate is at most, weakly significant (p < 0.01). However, a highly significant (p < 10(-5)), threefold difference in the rate of insertions/deletions (indels) exists between two recent studies in Schizosaccharomyces pombe. Upon investigation of the ancestral genome sequence for both studies, a possible anti-mutator allele was identified. The observed variation in indel rate may imply that the use of indel markers, such as microsatellites, for the investigation of genetic diversity within and among populations may be inappropriate because of the assumption of uniform mutation rate within a species.

Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level

BACKGROUND

The amount and structure of genetic diversity in dessert apple germplasm conserved at a European level is mostly unknown, since all diversity studies conducted in Europe until now have been performed on regional or national collections. Here, we applied a common set of 16 SSR markers to genotype more than 2,400 accessions across 14 collections representing three broad European geographic regions (North + East, West and South) with the aim to analyze the extent, distribution and structure of variation in the apple genetic resources in Europe.

RESULTS

A Bayesian model-based clustering approach showed that diversity was organized in three groups, although these were only moderately differentiated (FST = 0.031). A nested Bayesian clustering approach allowed identification of subgroups which revealed internal patterns of substructure within the groups, allowing a finer delineation of the variation into eight subgroups (FST = 0.044). The first level of stratification revealed an asymmetric division of the germplasm among the three groups, and a clear association was found with the geographical regions of origin of the cultivars. The substructure revealed clear partitioning of genetic groups among countries, but also interesting associations between subgroups and breeding purposes of recent cultivars or particular usage such as cider production. Additional parentage analyses allowed us to identify both putative parents of more than 40 old and/or local cultivars giving interesting insights in the pedigree of some emblematic cultivars.

CONCLUSIONS

The variation found at group and subgroup levels may reflect a combination of historical processes of migration/selection and adaptive factors to diverse agricultural environments that, together with genetic drift, have resulted in extensive genetic variation but limited population structure. The European dessert apple germplasm represents an important source of genetic diversity with a strong historical and patrimonial value. The present work thus constitutes a decisive step in the field of conservation genetics. Moreover, the obtained data can be used for defining a European apple core collection useful for further identification of genomic regions associated with commercially important horticultural traits in apple through genome-wide association studies.

Population Structure and Evolution after Speciation of the Hokkaido Salamander (Hynobius retardatus)

The Hokkaido salamander (Hynobius retardatus) is endemic to Hokkaido Island, Japan, and shows intriguing flexible phenotypic plasticity and regional morphological diversity. However, to date, allozymes and partial mitochondria DNA sequences have provided only an outline of its demographic histories and the pattern of its genetic diversification. To understand the finer details of the population structure of this species and its evolution since speciation, we genotyped five regional populations by using 12 recently developed microsatellite polymorphic markers. We found a clear population structure with low gene flow among the five populations, but a close genetic relationship between the Teshio and Kitami populations. Our demographic analysis suggested that Teshio and Erimo had the largest effective population sizes among the five populations. These findings regarding the population structure and demography of H. retardatus improve our understanding of the faunal phylogeography on Hokkaido Island and also provide fundamental genetic information that will be useful for future studies.

Genetic diversity estimates point to immediate efforts for conserving the endangered Tibetan sheep of India

Tibetan is a valuable Himalayan sheep breed classified as endangered. Knowledge of the level and distribution of genetic diversity in Tibetan sheep is important for designing conservation strategies for their sustainable survival and to preserve their evolutionary potential. Thus, for the first time, genetic variability in the Tibetan population was accessed with twenty five inter-simple sequence repeat markers. All the microsatellites were polymorphic and a total of 148 alleles were detected across these loci. The observed number of alleles across all the loci was more than the effective number of alleles and ranged from 3 (BM6506) to 11 (BM6526) with 5.920 ± 0.387 mean number of alleles per locus. The average observed heterozygosity was less than the expected heterozygosity. The observed and expected heterozygosity values ranged from 0.150 (BM1314) to 0.9 (OarCP20) with an overall mean of 0.473 ± 0.044 and from 0.329 (BM8125) to 0.885 (BM6526) with an overall mean 0.672 ± 0.030, respectively. The lower heterozygosity pointed towards diminished genetic diversity in the population. Thirteen microsatellite loci exhibited significant (P < 0.05) departures from the Hardy-Weinberg proportions in the population. The estimate of heterozygote deficiency varied from - 0.443 (OarCP20) to 0.668 (OarFCB128) with a mean positive value of 0.302 ± 0.057. A normal 'L' shaped distribution of mode-shift test and non-significant heterozygote excess on the basis of different models suggested absence of recent bottleneck in the existing Tibetan population. In view of the declining population of Tibetan sheep (less than 250) in the breeding tract, need of the hour is immediate scientific management of the population so as to increase the population hand in hand with retaining the founder alleles to the maximum possible extent.

Bulk development and stringent selection of microsatellite markers in the western flower thrips Frankliniella occidentalis

Recent improvements in next-generation sequencing technologies have enabled investigation of microsatellites on a genome-wide scale. Faced with a huge amount of candidates, the use of appropriate marker selection criteria is crucial. Here, we used the western flower thrips Frankliniella occidentalis for an empirical microsatellite survey and validation; 132,251 candidate microsatellites were identified, 92,102 of which were perfect. Dinucleotides were the most abundant category, while (AG)n was the most abundant motif. Sixty primer pairs were designed and validated in two natural populations, of which 30 loci were polymorphic, stable, and repeatable, but not all in Hardy-Weinberg equilibrium (HWE) and linkage equilibrium. Four marker panels were constructed to understand effect of marker selection on population genetic analyses: (i) only accept loci with single nucleotide insertions (SNI); (ii) only accept the most polymorphic loci (MP); (iii) only accept loci that did not deviate from HWE, did not show SNIs, and had unambiguous peaks (SS) and (iv) all developed markers (ALL). Although the MP panel resulted in microsatellites of highest genetic diversity followed by the SNI, the SS performed best in individual assignment. Our study proposes stringent criteria for selection of microsatellites from a large-scale number of genomic candidates for population genetic studies.

Diversity and population-genetic properties of copy number variations and multicopy genes in cattle

The diversity and population genetics of copy number variation (CNV) in domesticated animals are not well understood. In this study, we analysed 75 genomes of major taurine and indicine cattle breeds (including Angus, Brahman, Gir, Holstein, Jersey, Limousin, Nelore, and Romagnola), sequenced to 11-fold coverage to identify 1,853 non-redundant CNV regions. Supported by high validation rates in array comparative genomic hybridization (CGH) and qPCR experiments, these CNV regions accounted for 3.1% (87.5 Mb) of the cattle reference genome, representing a significant increase over previous estimates of the area of the genome that is copy number variable (∼2%). Further population genetics and evolutionary genomics analyses based on these CNVs revealed the population structures of the cattle taurine and indicine breeds and uncovered potential diversely selected CNVs near important functional genes, including AOX1, ASZ1, GAT, GLYAT, and KRTAP9-1. Additionally, 121 CNV gene regions were found to be either breed specific or differentially variable across breeds, such as RICTOR in dairy breeds and PNPLA3 in beef breeds. In contrast, clusters of the PRP and PAG genes were found to be duplicated in all sequenced animals, suggesting that subfunctionalization, neofunctionalization, or overdominance play roles in diversifying those fertility-related genes. These CNV results provide a new glimpse into the diverse selection histories of cattle breeds and a basis for correlating structural variation with complex traits in the future.

Circumpolar Genetic Structure and Recent Gene Flow of Polar Bears: A Reanalysis

Recently, an extensive study of 2,748 polar bears (Ursus maritimus) from across their circumpolar range was published in PLOS ONE, which used microsatellites and mitochondrial haplotypes to apparently show altered population structure and a dramatic change in directional gene flow towards the Canadian Archipelago-an area believed to be a future refugium for polar bears as their southernmost habitats decline under climate change. Although this study represents a major international collaborative effort and promised to be a baseline for future genetics work, methodological shortcomings and errors of interpretation undermine some of the study's main conclusions. Here, we present a reanalysis of this data in which we address some of these issues, including: (1) highly unbalanced sample sizes and large amounts of systematically missing data; (2) incorrect calculation of FST and of significance levels; (3) misleading estimates of recent gene flow resulting from non-convergence of the program BayesAss. In contrast to the original findings, in our reanalysis we find six genetic clusters of polar bears worldwide: the Hudson Bay Complex, the Western and Eastern Canadian Arctic Archipelago, the Western and Eastern Polar Basin, and-importantly-we reconfirm the presence of a unique and possibly endangered cluster of bears in Norwegian Bay near Canada's expected last sea-ice refugium. Although polar bears' abundance, distribution, and population structure will certainly be negatively affected by ongoing-and increasingly rapid-loss of Arctic sea ice, these genetic data provide no evidence of strong directional gene flow in response to recent climate change.

Are microsatellite fragment lengths useful for population-level studies? The case of Polygala lewtonii (Polygalaceae)

PREMISE OF THE STUDY

Microsatellites, despite being commonly used population-level markers, contain biases because scoring relies solely on fragment length. Their complexity can lead to homoplasy, the effects of which are poorly understood. We investigated the impact of using fragment lengths, repeats, or flanking region sequences on common population-level analyses.

METHODS

Five polymorphic microsatellite markers amplified across the central Florida scrub endemic Polygala lewtonii (Polygalaceae) and its close, widespread congener P. polygama. We genotyped 147 individuals of P. lewtonii and 156 of P. polygama, and sequenced the amplicons of four markers across all individuals. We ran basic statistics, spatial clustering analysis, historical demographics, and migration tests.

RESULTS

One population of intermediate morphology was genetically clearly identified as P. polygama, making it the southernmost population of this species. Statistics were comparable between the fragment length and repeat numbers, with some notable differences. Flanking regions exhibited surprisingly high polymorphism between species, and between geographically distant conspecific populations.

DISCUSSION

The increasing use of markers developed in one species and amplified in another is only a good practice if precautions are taken, notably the sequencing of the fragments between species and populations. Flanking region sequences are a useful marker at the interspecific level.

Historical and current introgression in a Mesoamerican hummingbird species complex: a biogeographic perspective

The influence of geologic and Pleistocene glacial cycles might result in morphological and genetic complex scenarios in the biota of the Mesoamerican region. We tested whether berylline, blue-tailed and steely-blue hummingbirds, Amazilia beryllina, Amazilia cyanura and Amazilia saucerottei, show evidence of historical or current introgression as their plumage colour variation might suggest. We also analysed the role of past and present climatic events in promoting genetic introgression and species diversification. We collected mitochondrial DNA (mtDNA) sequence data and microsatellite loci scores for populations throughout the range of the three Amazilia species, as well as morphological and ecological data. Haplotype network, Bayesian phylogenetic and divergence time inference, historical demography, palaeodistribution modelling, and niche divergence tests were used to reconstruct the evolutionary history of this Amazilia species complex. An isolation-with-migration coalescent model and Bayesian assignment analysis were assessed to determine historical introgression and current genetic admixture. mtDNA haplotypes were geographically unstructured, with haplotypes from disparate areas interdispersed on a shallow tree and an unresolved haplotype network. Assignment analysis of the nuclear genome (nuDNA) supported three genetic groups with signs of genetic admixture, corresponding to: (1) A. beryllina populations located west of the Isthmus of Tehuantepec; (2) A. cyanura populations between the Isthmus of Tehuantepec and the Nicaraguan Depression (Nuclear Central America); and (3) A. saucerottei populations southeast of the Nicaraguan Depression. Gene flow and divergence time estimates, and demographic and palaeodistribution patterns suggest an evolutionary history of introgression mediated by Quaternary climatic fluctuations. High levels of gene flow were indicated by mtDNA and asymmetrical isolation-with-migration, whereas the microsatellite analyses found evidence for three genetic clusters with distributions corresponding to isolation by the Isthmus of Tehuantepec and the Nicaraguan Depression and signs of admixture. Historical levels of migration between genetically distinct groups estimated using microsatellites were higher than contemporary levels of migration. These results support the scenario of secondary contact and range contact during the glacial periods of the Pleistocene and strongly imply that the high levels of structure currently observed are a consequence of the limited dispersal of these hummingbirds across the isthmus and depression barriers.

Genetic Distinctiveness of Rye In situ Accessions from Portugal Unveils a New Hotspot of Unexplored Genetic Resources

Rye (Secale cereale L.) is a cereal crop of major importance in many parts of Europe and rye breeders are presently very concerned with the restrict pool of rye genetic resources available. Such narrowing of rye genetic diversity results from the presence of "Petkus" pool in most modern rye varieties as well as "Petkus" × "Carsten" heterotic pool in hybrid rye breeding programs. Previous studies on rye's genetic diversity revealed moreover a common genetic background on landraces (ex situ) and cultivars, regardless of breeding level or geographical origin. Thus evaluation of in situ populations is of utmost importance to unveil "on farm" diversity, which is largely undervalued. Here, we perform the first comprehensive assessment of rye's genetic diversity and population structuring using cultivars, ex situ landraces along a comprehensive sampling of in situ accessions from Portugal, through a molecular-directed analysis using SSRs markers. Rye genetic diversity and population structure analysis does not present any geographical trend but disclosed marked differences between genetic backgrounds of in situ accessions and those of cultivars/ex situ collections. Such genetic distinctiveness of in situ accessions highlights their unexplored potential as new genetic resources, which can be used to boost rye breeding strategies and the production of new varieties. Overall, our study successfully demonstrates the high prospective impact of comparing genetic diversity and structure of cultivars, ex situ, and in situ samples in ascertaining the status of plant genetic resources (PGR).

Evolution of Nine Microsatellite Loci in the Fungus Fusarium oxysporum

The evolution of nine microsatellites and one minisatellite was investigated in the fungus Fusarium oxysporum and sister taxa Fusarium redolens and Fusarium verticillioides. Compared to other organisms, fungi have been reported to contain fewer and less polymorphic microsatellites. Mutational patterns over evolutionary time were studied for these ten loci by mapping changes in core repeat numbers onto a phylogeny based on the sequence of the conserved translation elongation factor 1-α gene. The patterns of microsatellite formation, expansion, and interruption by base substitutions were followed across the phylogeny, showing that these loci are evolving in a manner similar to that of microsatellites in other eukaryotes. Most mutations could be fit to a stepwise mutation model, but a few appear to have involved multiple repeat units. No evidence of gene conversion was seen at the minisatellite locus, which may also be mutating by replication slippage. Some homoplastic numbers of repeat units were observed for these loci, and polymorphisms in the regions flanking the microsatellites may provide better genetic markers for population genetics studies of these species.

Genome-wide discovery and development of polymorphic microsatellites from Leishmania panamensis parasites circulating in central Panama

Background

The parasite Leishmania panamensis is the main cause of leishmaniasis in Panama. The disease is largely uncontrolled, with a rising incidence and no appropriate control measures. While microsatellites are considered some of the best genetic markers to study population genetics and molecular epidemiology in these and other parasites, none has been developed for L. panamensis.

Findings

Here we have developed and tested a new panel of microsatellites for this species, based on high-throughput genome-wide screening. The new set of microsatellites is composed of seventeen loci, mainly spanning trinucleotide or longer motifs. We have evaluated the sensitivity and specificity of the panel based on a sample of 27 isolates obtained from cutaneous leishmaniasis patients from central Panama and also several reference species from both L. (Leishmania) and L. (Viannia) subgenera. The genetic equilibrium was assessed both intra- and inter-loci, while the reproductive mode was evaluated using several tests. The new SSR panel shows high polymorphism and sensitivity, as well as good specificity. The preliminary data described here for L. panamensis suggest extensive departure from Hardy-Weinberg proportions, significant linkage disequilibrium and strong deficit of heterozygotes. Several recombination tests involving multilocus linkage disequilibrium and a phylogenetic approach allowed rejection of frequent recombination in our dataset.

Conclusions

The genome-wide strategy described here proved to be useful to identify and test new polymorphic SSR loci in Leishmania. The new panel of polymorphic microsatellites is a valuable contribution to the existing molecular markers for the study of genetic structure and other aspects of this important species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1153-2) contains supplementary material, which is available to authorized users.

The mind, the lab, and the field: Three kinds of populations in scientific practice

Scientists use models to understand the natural world, and it is important not to conflate model and nature. As an illustration, we distinguish three different kinds of populations in studies of ecology and evolution: theoretical, laboratory, and natural populations, exemplified by the work of R. A. Fisher, Thomas Park, and David Lack, respectively. Biologists are rightly concerned with all three types of populations. We examine the interplay between these different kinds of populations, and their pertinent models, in three examples: the notion of "effective" population size, the work of Thomas Park on Tribolium populations, and model-based clustering algorithms such as Structure. Finally, we discuss ways to move safely between three distinct population types while avoiding confusing models and reality.

Genetic diversity and population structure of Polistes nimpha based on DNA microsatellite markers

The Eurasiatic Polistes nimpha belongs to primitively eusocial wasps for which no data are available on its population's genetic structure and relatedness/relationships of individuals. The purpose of this research is to determine the amplification efficiency in P. nimpha of microsatellite primers developed for P. dominula and using these primers, to explore genetic diversity, population structure and relatedness/relationship of P. nimpha in the context of its reproductive options. Eight out of twelve microsatellite markers analyzed on 59 individuals (pupae and larvae) were polymorphic (mean PIC = 0.545) and mutated following the stepwise mutation model. The Bayesian clustering method gave the probability of >0.898 of there being 10 clusters within the pooled sample of 15 nests. In two or three nest clusters, full- and/or half-siblings and unrelated individuals occurred. A significant correlation between genetic and geographic distances was detected. There are three main possibilities that come into play to explain our genetic results and direct field observations: cooperative nest foundation, visitations, and/or usurpation events. So far there is no conclusive evidence to exclude or support any of these possibilities.

To what extent do human-altered landscapes retain population connectivity? Historical changes in gene flow of wetland fish Pungitius pungitius

Understanding how human-altered landscapes affect population connectivity is valuable for conservation planning. Natural connectivity among wetlands, which is maintained by floods, is disappearing owing to farmland expansion. Using genetic data, we assessed historical changes in the population connectivity of the ninespine stickleback within a human-altered wetland system. We predicted that: (i) the contemporary gene flow maintained by the artificial watercourse network may be restricted to a smaller spatial scale compared with the gene flow preceding alteration, and (ii) the contemporary gene flow is dominated by the downstream direction owing to the construction of low-head barriers. We evaluated the potential source population in both timescales. Seventeen studied populations were grouped into four genetically different clusters, and we estimated the migration rates among these clusters. Contemporary migration was restricted to between neighbouring clusters, although a directional change was not detected. Furthermore, we consistently found the same potential source cluster, from past to present, characterized by large amounts of remnant habitats connected by artificial watercourses. These findings highlight that: (i) artificial connectivity can sustain the short-distance connectivity of the ninespine stickleback, which contributes to maintaining the potential source populations; however, (ii) population connectivity throughout the landscape has been prevented by agricultural developments.

Genetic diversity and relationship of Indian cattle inferred from microsatellite and mitochondrial DNA markers

BACKGROUND

Indian agriculture is an economic symbiosis of crop and livestock production with cattle as the foundation. Sadly, the population of indigenous cattle (Bos indicus) is declining (8.94% in last decade) and needs immediate scientific management. Genetic characterization is the first step in the development of proper management strategies for preserving genetic diversity and preventing undesirable loss of alleles. Thus, in this study we investigated genetic diversity and relationship among eleven Indian cattle breeds using 21 microsatellite markers and mitochondrial D loop sequence.

RESULTS

The analysis of autosomal DNA was performed on 508 cattle which exhibited sufficient genetic diversity across all the breeds. Estimates of mean allele number and observed heterozygosity across all loci and population were 8.784 ± 0.25 and 0.653 ± 0.014, respectively. Differences among breeds accounted for 13.3% of total genetic variability. Despite high genetic diversity, significant inbreeding was also observed within eight populations. Genetic distances and cluster analysis showed a close relationship between breeds according to proximity in geographic distribution. The genetic distance, STRUCTURE and Principal Coordinate Analysis concluded that the Southern Indian Ongole cattle are the most distinct among the investigated cattle populations. Sequencing of hypervariable mitochondrial DNA region on a subset of 170 cattle revealed sixty haplotypes with haplotypic diversity of 0.90240, nucleotide diversity of 0.02688 and average number of nucleotide differences as 6.07407. Two major star clusters for haplotypes indicated population expansion for Indian cattle.

CONCLUSIONS

Nuclear and mitochondrial genomes show a similar pattern of genetic variability and genetic differentiation. Various analyses concluded that the Southern breed 'Ongole' was distinct from breeds of Northern/ Central India. Overall these results provide basic information about genetic diversity and structure of Indian cattle which should have implications for management and conservation of indicine cattle diversity.

Dendritic connectivity shapes spatial patterns of genetic diversity: a simulation-based study

Landscape features notoriously affect spatial patterns of biodiversity. For instance, in dendritic ecological networks (such as river basins), dendritic connectivity has been proposed to create unique spatial patterns of biodiversity. Here, we compared genetic datasets simulated under a lattice-like, a dendritic and a circular landscape to test the influence of dendritic connectivity on neutral genetic diversity. The circular landscape had a level of connectivity similar to that of the dendritic landscape, so as to isolate the influence of dendricity on genetic diversity. We found that genetic diversity and differentiation varied strikingly among the three landscapes. For instance, the dendritic landscape generated higher total number of alleles and higher global Fst than the lattice-like landscape, and these indices also varied between the dendritic and the circular landscapes, suggesting an effect of dendricity. Furthermore, in the dendritic landscape, allelic richness was higher in highly connected demes (e.g. confluences in rivers) than in low-connected demes (e.g. upstream and downstream populations), which was not the case in the circular landscape, hence confirming the major role of dendricity. This led to bell-shaped distributions of allelic richness along an upstream-downstream gradient. Conversely, genetic differentiation (Fst ) was lower in highly than in low-connected demes (which was not observed in circular landscape), and significant patterns of isolation by distance (IBD) were also observed in the dendritic landscape. We conclude that in dendritic networks, the combined influence of dendricity and connectivity generates unique spatial patterns of neutral genetic diversity, which has implications for population geneticists and conservationists.