The report of my death was an exaggeration: A review for researchers using microsatellites in the 21st century

MIPs: multi-locus intron polymorphisms in species identification and population genomics

The study of species groups in which the presence of interspecific hybridization or introgression phenomena is known or suspected involves analysing shared bi-parentally inherited molecular markers. Current methods are based on different categories of markers among which the classical microsatellites or the more recent genome wide approaches for the analyses of thousands of SNPs or hundreds of microhaplotypes through high throughput sequencing. Our approach utilizes intron-targeted amplicon sequencing to characterise multi-locus intron polymorphisms (MIPs) and assess genetic diversity. These highly variable intron regions, combined with inter-specific transferable loci, serve as powerful multiple-SNP markers potentially suitable for various applications, from species and hybrid identification to population comparisons, without prior species knowledge. We developed the first panel of MIPs highly transferable across fish genomes, effectively distinguishing between species, even those closely related, and populations with different structures. MIPs offer versatile, hypervariable nuclear markers and promise to be especially useful when multiple nuclear loci must be genotyped across different species, such as for the monitoring of interspecific hybridization. Moreover, the relatively long sequences obtained ease the development of single-locus PCR-based diagnostic markers. This method, here demonstrated in teleost fishes, can be readily applied to other taxa, unlocking a new source of genetic variation.

Development of microsatellite markers for the soft tick Ornithodoros phacochoerus

BACKGROUND

Soft ticks of the genus Ornithodoros are responsible for the maintenance and transmission of the African swine fever (ASF) virus in the sylvatic and domestic viral cycles in Southern Africa. They are also the main vectors of the Borrelia species causing relapsing fevers. Currently, no genetic markers are available for Afrotropical Ornithodoros ticks. As ASF spreads globally, such markers are needed to assess the role of ticks in the emergence of new outbreaks. The aim of this study is to design microsatellite markers that could be used for ticks of the Ornithodoros moubata complex, particularly Ornithodoros phacochoerus, to assess population structure and tick movements in ASF endemic areas.

METHODS

A total of 151 markers were designed using the O. moubata and O. porcinus genomes after elimination of repeated sequences in the genomes. All designed markers were tested on O. phacochoerus and O. porcinus DNA to select the best markers.

RESULTS

A total of 24 microsatellite markers were genotyped on two populations of O. phacochoerus and on individuals from four other Ornithodoros species. Nineteen markers were selected to be as robust as possible for population genetic studies on O. phacochoerus.

CONCLUSIONS

The microsatellite markers developed here represent the first genetic tool to study nidicolous populations of O. phacochoerus.

Ex situ conservation of two rare oak species using microsatellite and SNP markers

Plant collections held by botanic gardens and arboreta are key components of ex situ conservation. Maintaining genetic diversity in such collections allows them to be used as resources for supplementing wild populations. However, most recommended minimum sample sizes for sufficient ex situ genetic diversity are based on microsatellite markers, and it remains unknown whether these sample sizes remain valid in light of more recently developed next-generation sequencing (NGS) approaches. To address this knowledge gap, we examine how ex situ conservation status and sampling recommendations differ when derived from microsatellites and single nucleotide polymorphisms (SNPs) in garden and wild samples of two threatened oak species. For Quercus acerifolia, SNPs show lower ex situ representation of wild allelic diversity and slightly lower minimum sample size estimates than microsatellites, while results for each marker are largely similar for Q. boyntonii. The application of missing data filters tends to lead to higher ex situ representation, while the impact of different SNP calling approaches is dependent on the species being analyzed. Measures of population differentiation within species are broadly similar between markers, but larger numbers of SNP loci allow for greater resolution of population structure and clearer assignment of ex situ individuals to wild source populations. Our results offer guidance for future ex situ conservation assessments utilizing SNP data, such as the application of missing data filters and the usage of a reference genome, and illustrate that both microsatellites and SNPs remain viable options for botanic gardens and arboreta seeking to ensure the genetic diversity of their collections.

Assessment of the Genetic Diversity and Structure of the Korean Endemic Freshwater Fish Microphysogobio longidorsalis (Gobioninae) Using Microsatellite Markers: A First Glance from Population Genetics

Microphysogobio longidorsalis is endemic to South Korea and inhabits small areas of the Namhangang, Bukhangang, and Imjingang Rivers in the Hangang River water system. Endemic species usually are more vulnerable than species with a wide distribution. Notably, there is a lack of basic conservation data for M. longidorsalis. We analyzed 19 microsatellite loci in six populations of M. longidorsalis in South Korea to characterize their population structure and genetic diversity. The genetic diversity of the microsatellites was 0.741-0.779, which is lower than that of other freshwater fishes. The pairwise genetic differentiation of microsatellite (FST) values ranged from 0.007 to 0.041, suggesting low genetic differentiation between the populations. The Jojongicheon stream population (CP) had an effective population size of <100. Therefore, conservation efforts are required to prevent inbreeding depression in M. longidorsalis. Discriminant analysis of principal components showed that the Hangang River water system would be a single management unit (MU). Our findings provide fundamental genetic insights for the formulation of conservation strategies for M. longidorsalis.

Genetic Diversity of Apple Clonal Rootstocks from the Collection of the Michurinsk State Agrarian University Based on SSR Markers

The Michurinsk State Agrarian University (Michurinsk SAU) is one of the leading centers for breeding apple clonal rootstocks. A diverse collection of apple rootstocks, founded in 1930s by V.I. Budagovsky, is maintained at the Michurinsk SAU. In the present study, 87 rootstocks from this collection were analyzed using 18 SSR markers to assess their genetic diversity and relatedness. The detected polymorphism level was rather high compared to the previous estimates of apple rootstock genetic variability. A total of 199 alleles were detected with an average of 11.1 alleles per locus. Among the detected alleles, 67 (33.67%) were rare and 43 (21.61%) were unique. The average PIC value was 0.73, and the expected and observed heterozygosity averaged 0.76 and 0.69, respectively. All the studied accessions except two could be identified with the used marker set. Cluster analysis revealed several groups according to the rootstocks' pedigrees and genetic origin. Furthermore, Structure analysis revealed two main groups of the studied rootstock accessions. No significant differentiation of the studied sample according to dwarfing ability was detected, while weak differentiation was detected according to leaf color. SSR genotyping data can be used for rootstock fingerprinting and pedigree verification and will facilitate collection management. In addition, data on the genetic diversity and structure of the studied collection may be useful for further development of the Michurinsk SAU rootstock breeding program.

Population size affected by environmental variability impacts genetics, traits, and plant performance in Trifolium montanum L

Population size, genetic diversity, and performance have fundamental importance for ecology, evolution, and nature conservation of plant species. Despite well-studied relationships among environmental, genetic, and intraspecific trait variation (ITV), the influence of population size on these aspects is less understood. To assess the sources of population size variation, but also its impact on genetic, functional trait, and performance aspects, we conducted detailed population size estimations, assessed 23 abiotic and biotic environmental habitat factors, performed population genetic analyses using nine microsatellite markers, and recorded nine functional traits based on 260 Trifolium montanum individuals from 13 semi-dry grassland locations of Central Europe. Modern statistical analyses based on a multivariate framework (path analysis) with preselected linear regression models revealed that the variation of abiotic factors (in contrast to factors per se) almost completely, significantly explained fluctuations in population size (R 2 = .93). In general, abiotic habitat variation (heterogeneity) was not affected by habitat area. Population size significantly explained genetic diversity (N A: R 2 = .42, H o: R 2 = .67, H e: R 2 = .43, and I: R 2 = .59), inbreeding (F IS: R 2 = .35), and differentiation (G ST: R 2 = .20). We also found that iFDCV (ITV) was significantly explained by abiotic habitat heterogeneity, and to a lesser extent by genetic diversity H e (R 2 = .81). Nevertheless, habitat heterogeneity did not statistically affect genetic diversity. This may be due to the use of selectively neutral microsatellite markers, and possibly by insufficient abiotic selective pressures on habitats examined. Small T. montanum populations in nonoptimal habitats were characterized by reduced genetic and functional trait diversity, and elevated genetic inbreeding and differentiation. This indicates reduced adaptability to current and future environmental changes. The long-term survival of small populations with reduced genetic diversity and beginning inbreeding will be highly dependent on habitat protection and adequate land-use actions.

Characterization of 21 microsatellite loci for the precocious, grass-shrimp trematode Alloglossidium renale

We developed microsatellite markers to use in studying the population genetics of the trematode Alloglossidium renale, a fluke with a precocious life cycle where sexual maturation occurs in a grass shrimp. Among 21 tested loci in a Mississippi population sample, 14 were polymorphic, 12 of which significantly deviated from Hardy-Weinberg Equilibrium (HWE). We estimated identity disequilibrium (ID) to confirm whether the deviations from HWE were due to significant amounts of selfing or due to technical factors. The selfing rate derived from F_IS was 86.6%, whereas the selfing rate obtained by ID was 83.9%, indicating that the deviation in HWE was due to a high amount of selfing within the population. These markers will be useful for ecological and evolutionary studies of A. renale especially in relation to the interplay of hermaphroditic mating systems, inbreeding depression, and transmission dynamics.

Unravelling the genetic diversity and population structure of common walnut in the Iranian Plateau

BACKGROUND

Common walnut (Juglans regia L.) has a long cultivation history, given its highly valuable wood and rich nutritious nuts. The Iranian Plateau has been considered as one of the last glaciation refugia and a centre of origin and domestication for the common walnut. However, a prerequisite to conserve or utilize the genetic resources of J. regia in the plateau is a comprehensive evaluation of the genetic diversity that is conspicuously lacking. In this regard, we used 31 polymorphic simple sequence repeat (SSR) markers to delineate the genetic variation and population structure of 508 J. regia individuals among 27 populations from the Iranian Plateau.

RESULTS

The SSR markers expressed a high level of genetic diversity (H_O = 0.438, and H_E = 0.437). Genetic differentiation among the populations was moderate (F_ST = 0.124), and genetic variation within the populations (79%) significantly surpassed among populations (21%). The gene flow (N_m = 1.840) may have remarkably influenced the population genetic structure of J. regia, which can be attributed to anthropological activities and wind dispersal of pollen. The STRUCTURE analysis divided the 27 populations into two main clusters. Comparing the neighbor-joining and principal coordinate analysis dendrograms and the Bayesian STRUCTURE analysis revealed the general agreement between the population subdivisions and the genetic relationships among the populations. However, a few geographically close populations dispersed into different clusters. Further, the low genetic diversity of the Sulaymaniyah (SMR) population of Iraq necessitates urgent conservation by propagation and seedling management or tissue culture methods; additionally, we recommend the indispensable preservation of the Gonabad (RGR) and Arak (AKR) populations in Iran.

CONCLUSIONS

These results reflected consistent high geographical affinity of the accession across the plateau. Our findings suggest that gene flow is a driving factor influencing the genetic structure of J. regia populations, whereas ecological and geological variables did not act as strong barriers. Moreover, the data reported herein provide new insights into the population structure of J. regia germplasm, which will help conserve genetic resources for the future, hence improving walnut breeding programs' efficiency.

Plant mitochondrial introns as genetic markers - conservation and variation

Plant genomes are comprised of nuclear, plastid and mitochondrial components characterized by different patterns of inheritance and evolution. Genetic markers from the three genomes provide complementary tools for investigations of inheritance, genetic relationships and phenotypic contributions. Plant mitochondrial genomes are challenging for universal marker development because they are highly variable in terms of size, gene order and intergenic sequences and highly conserved with respect to protein-coding sequences. PCR amplification of introns with primers that anneal to conserved, flanking exons is effective for the development of polymorphic nuclear genome markers. The potential for plant mitochondrial intron polymorphisms to distinguish between congeneric species or intraspecific varieties has not been systematically investigated and is possibly constrained by requirements for intron secondary structure and interactions with co-evolved organelle intron splicing factors. To explore the potential for broadly applicable plant mitochondrial intron markers, PCR primer sets based upon conserved sequences flanking 11 introns common to seven angiosperm species were tested across a range of plant orders. PCR-amplified introns were screened for indel polymorphisms among a group of cross-compatible Citrus species and relatives; two Raphanus sativus mitotypes; representatives of the two Phaseolus vulgaris gene pools; and congeneric pairs of Cynodon, Cenchrus, Solanum, and Vaccinium species. All introns were successfully amplified from each plant entry. Length polymorphisms distinguishable by gel electrophoresis were common among genera but infrequent within genera. Sequencing of three introns amplified from 16 entries identified additional short indel polymorphisms and nucleotide substitutions that separated Citrus, Cynodon, Cenchrus and Vaccinium congeners, but failed to distinguish Solanum congeners or representatives of the Phaseolus vulgaris major gene pools. The ability of primer sets to amplify a wider range of plant species' introns and the presence of intron polymorphisms that distinguish congeners was confirmed by in silico analysis. While mitochondrial intron variation is limited in comparison to nuclear introns, these exon-based primer sets provide robust tools for the amplification of mitochondrial introns across a wide range of plant species wherein useful polymorphisms can be identified.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Phylogeographic structure of common sage (Salvia officinalis L.) reveals microrefugia throughout the Balkans and colonizations of the Apennines

Studying the population-genetic and phylogeographic structures of a representative species of a particular geographical region can not only provide us with information regarding its evolutionary history, but also improve our understanding of the evolutionary processes underlying the patterns of species diversity in that area. By analysing eight highly polymorphic microsatellite loci and two chloroplast DNA regions, we have investigated the influence of Pleistocene climate fluctuations on the evolutionary history of Salvia officinalis L. (common sage). The populations with the highest genetic diversity were located in the central parts of the Balkan distribution range. A large group of closely related haplotypes was distributed throughout the Balkans and the central Apennines, while the private lineage occupied the southern Apennines. In addition, two highly differentiated lineages were scattered only over the Balkans. The results suggest that a single refugium of the studied species from the last glacial period was located in the central part of the range in the Balkans. Numerous microrefugia, probably spanning several glaciation cycles, were scattered across the Balkans, while colonisation of the Apennines from the Balkans occurred at least on two occasions.

Tracking population genetic signatures of local extinction with herbarium specimens

BACKGROUND AND AIMS

Habitat degradation and landscape fragmentation dramatically lower population sizes of rare plant species. Decreasing population sizes may, in turn, negatively affect genetic diversity and reproductive fitness, which can ultimately lead to local extinction of populations. Although such extinction vortex dynamics have been postulated in theory and modelling for decades, empirical evidence from local extinctions of plant populations is scarce. In particular, comparisons between current vs. historical genetic diversity and differentiation are lacking despite their potential to guide conservation management.

METHODS

We studied the population genetic signatures of the local extinction of Biscutella laevigata subsp. gracilis populations in Central Germany. We used microsatellites to genotype individuals from 15 current populations, one ex situ population, and 81 herbarium samples from five extant and 22 extinct populations. In the current populations, we recorded population size and fitness proxies, collected seeds for a germination trial and conducted a vegetation survey. The latter served as a surrogate for habitat conditions to study how habitat dissimilarity affects functional connectivity among the current populations.

KEY RESULTS

Bayesian clustering revealed similar gene pool distribution in current and historical samples but also indicated that a distinct genetic cluster was significantly associated with extinction probability. Gene flow was affected by both the spatial distance and floristic composition of population sites, highlighting the potential of floristic composition as a powerful predictor of functional connectivity which may promote decision-making for reintroduction measures. For an extinct population, we found a negative relationship between sampling year and heterozygosity. Inbreeding negatively affected germination.

CONCLUSIONS

Our study illustrates the usefulness of historical DNA to study extinction vortices in threatened species. Our novel combination of classical population genetics together with data from herbarium specimens, an ex situ population and a germination trial underlines the need for genetic rescue measures to prevent extinction of B. laevigata in Central Germany.

Characterization of polymorphic microsatellite markers for the Japanese endangered land snail Mandarina

Objective

Mandarina is an endangered land snail genus of the oceanic Ogasawara archipelago. On Chichijima Island, the largest inhabited island in Ogasawara, this genus is almost extinct in the wild due to predation by invasive species. Although ex situ conservation programs started in 2010, genetic diversity and population structure remain unclear due to a lack of genetic markers with sufficient genetic variation. In this study, we designed polymorphic microsatellite markers of Mandarina to enable genetic analysis and to develop appropriate conservation plans.

Results

Twenty-three polymorphic microsatellite markers were identified from the genomic DNA of wild samples of Mandarina mandarina. We assessed the genetic diversity of each marker. In 16 markers, neither linkage disequilibrium nor deviation from Hardy–Weinberg equilibrium was detected. These 16 markers were tested for multiplex PCR using low-density DNA extracted non-lethally from captive samples of M. mandarina, M. chichijimana and M. suenoae. Of the 16 markers, 15, 12 and 9 were usable for multiplex PCR, respectively. Genetic analysis using these microsatellite loci will be an important resource for the conservation of Mandarina.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-06147-4.

A comparative study of population genetic structure reveals patterns consistent with selection at functional microsatellites in common sunflower

Microsatellites, also known as short tandem repeats (STRs), have long been considered non-functional, neutrally evolving regions of the genome. Recent findings suggest that they can function as drivers of rapid adaptive evolution. Previous work on the common sunflower identified 479 transcribed microsatellites where allele length significantly correlates with gene expression (eSTRs) in a stepwise manner. Here, a population genetic approach is used to test whether eSTR allele length variation is under selection. Genotypic variation among and within populations at 13 eSTRs was compared with that at 19 anonymous microsatellites in 672 individuals from 17 natural populations of sunflower from across a cline running from Saskatchewan to Oklahoma (distance of approximately 1600 km). Expected heterozygosity, allelic richness, and allelic diversity were significantly lower at eSTRs, a pattern consistent with higher relative rates of purifying selection. Further, an analysis of variation in microsatellite allele lengths (lnRV), and heterozygosities (lnRH), indicate recent selective sweeps at the eSTRs. Mean microsatellite allele lengths at four eSTRs within populations are significantly correlated with latitude consistent with the predictions of the tuning-knob model which predicts stepwise relationships between microsatellite allele length and phenotypes. This finding suggests that shorter or longer alleles at eSTRs may be favored in climatic extremes. Collectively, our results imply that eSTRs are likely under selection and that they may be playing a role in facilitating local adaptation across a well-defined cline in the common sunflower.

Revisiting the comparative phylogeography of unglaciated eastern North America: 15 years of patterns and progress

In a landmark comparative phylogeographic study, “Comparative phylogeography of unglaciated eastern North America,” Soltis et al. (Molecular Ecology, 2006, 15, 4261) identified geographic discontinuities in genetic variation shared across taxa occupying unglaciated eastern North America and proposed several common biogeographical discontinuities related to past climate fluctuations and geographic barriers. Since 2006, researchers have published many phylogeographical studies and achieved many advances in genotyping and analytical techniques; however, it is unknown how this work has changed our understanding of the factors shaping the phylogeography of eastern North American taxa. We analyzed 184 phylogeographical studies of eastern North American taxa published between 2007 and 2019 to evaluate: (1) the taxonomic focus of studies and whether a previously detected taxonomic bias towards studies focused on vertebrates has changed over time, (2) the extent to which studies have adopted genotyping technologies that improve the resolution of genetic groups (i.e., NGS DNA sequencing) and analytical approaches that facilitate hypothesis‐testing (i.e., divergence time estimation and niche modeling), and (3) whether new studies support the hypothesized biogeographic discontinuities proposed by Soltis et al. (Molecular Ecology, 2006, 15, 4261) or instead support new, previously undetected discontinuities. We observed little change in taxonomic focus over time, with studies still biased toward vertebrates. Although many technological and analytical advances became available during the period, uptake was slow and they were employed in only a small proportion of studies. We found variable support for previously identified discontinuities and identified one new recurrent discontinuity. However, the limited resolution and taxonomic breadth of many studies hindered our ability to clarify the most important climatological or geographical factors affecting taxa in the region. Broadening the taxonomic focus to include more non‐vertebrate taxa, employing technologies that improve genetic resolution, and using analytical approaches that improve hypothesis testing are necessary to strengthen our inference of the forces shaping the phylogeography of eastern North America.

Genetic diversity and sex‐biased dispersal in the brown spotted pitviper ( Protobothrops mucrosquamatus ): Evidence from microsatellite markers

Dispersal plays a vital role in the geographical distribution, population genetic structure, quantity dynamics, and evolution of a species. Sex‐biased dispersal is common among vertebrates and many studies have documented a tendency toward male‐biased dispersal in mammals and female‐biased dispersal in birds. However, dispersal patterns in reptiles remain poorly understood. In this study, we explored the genetic diversity and dispersal patterns of the widely distributed Asian pitviper Protobothrops mucrosquamatus. In total, 16 polymorphic microsatellite loci were screened in 150 snakes (48 males, 44 females, 58 samples without sex information) covering most of their distribution. Microsatellite analysis revealed high genetic diversity in P. mucrosquamatus. Bayesian clustering of population assignment identified two major clusters for all populations, somewhat inconsistent with the mitochondrial DNA phylogeny of P. mucrosquamatus reported in previous research. Analyses based on 92 sex‐determined and 37 samples of P. mucrosquamatus from three small sites in Sichuan, China (Mingshan, Yibin, and Zizhong) consistently suggested female‐biased dispersal in P. mucrosquamatus, which is the first example of this pattern in snakes. The female‐biased dispersal patterns in P. mucrosquamatus may be explained by local resource competition.

A pipeline for effectively developing highly polymorphic simple sequence repeats markers based on multi-sample genomic data

Simple sequence repeats (SSRs) are widely used genetic markers in ecology, evolution, and conservation even in the genomics era, while a general limitation to their application is the difficulty of developing polymorphic SSR markers. Next-generation sequencing (NGS) offers the opportunity for the rapid development of SSRs; however, previous studies developing SSRs using genomic data from only one individual need redundant experiments to test the polymorphisms of SSRs. In this study, we designed a pipeline for the rapid development of polymorphic SSR markers from multi-sample genomic data. We used bioinformatic software to genotype multiple individuals using resequencing data, detected highly polymorphic SSRs prior to experimental validation, significantly improved the efficiency and reduced the experimental effort. The pipeline was successfully applied to a globally threatened species, the brown eared-pheasant (Crossoptilon mantchuricum), which showed very low genomic diversity. The 20 newly developed SSR markers were highly polymorphic, the average number of alleles was much higher than the genomic average. We also evaluated the effect of the number of individuals and sequencing depth on the SSR mining results, and we found that 10 individuals and ~10X sequencing data were enough to obtain a sufficient number of polymorphic SSRs, even for species with low genetic diversity. Furthermore, the genome assembly of NGS data from the optimal number of individuals and sequencing depth can be used as an alternative reference genome if a high-quality genome is not available. Our pipeline provided a paradigm for the application of NGS technology to mining and developing molecular markers for ecological and evolutionary studies.

Phylogeography and population genetic structure of the European roe deer in Switzerland following recent recolonization

In the early 1800s, the European roe deer (Capreolus capreolus) was probably extirpated from Switzerland, due to overhunting and deforestation. After a federal law was enacted in 1875 to protect lactating females and young, and limiting the hunting season, the roe deer successfully recovered and recolonized Switzerland. In this study, we use mitochondrial DNA and nuclear DNA markers to investigate the recolonization and assess contemporary genetic structure in relation to broad topographic features, in order to understand underlying ecological processes, inform future roe deer management strategies, and explore the opportunity for development of forensic traceability tools. The results concerning the recolonization origin support natural, multidirectional immigration from neighboring countries. We further demonstrate that there is evidence of weak genetic differentiation within Switzerland among topographic regions. Finally, we conclude that the genetic data support the recognition of a single roe deer management unit within Switzerland, within which there is a potential for broad-scale geographic origin assignment using nuclear markers to support law enforcement.

Comparing the population history of Neotropical annual species: The role of climate change and hybridization between Tithonia tubaeformis and T. rotundifolia (Asteraceae)

Historical climate changes and interspecific gene flow have played an important role in shaping the distribution and genetic diversity of the biota in the Neotropics. In this study, we explored the role of both the Pleistocene climate changes and hybridization on the contemporary geographic structure of two Neotropical weedy species: Tithonia tubaeformis and T. rotundifolia. Distribution shifts under past and current climate conditions were explored through ecological niche modelling (ENM). We then tested these hypotheses using chloroplast microsatellite (cpSSR) data in T. tubaeformis and compared the patterns deduced with those previously reported for T. rotundifolia using the same cpSSR loci. Lastly, we searched for shared haplotypes between species. Both species exhibited significant downwards altitudinal shifts during the last interglacial (LIG) and the Last Glacial Maximum (LGM). For instance, T. rotundifolia showed large suitable habitat areas since the LIG at xeric conditions in western Mesoamerica. Tithonia tubaeformis showed a ~five-fold range contraction during the LIG compared to current climate conditions, followed by a range expansion in the LGM. Despite the large shared refugial areas predicted through ENM, we found a low number of shared haplotypes, suggesting a minor role of hybridization in shaping the geographic structure of these species. Our results provide additional patterns of the population history of the northern Neotropics during the Quaternary, and we suggest that weedy widespread species are a well-suited group for the study of the effects of historic climatic changes on the biota of this region.

Development of a multiplex microsatellite marker set for the study of the solitary red mason bee, Osmia bicornis (Megachilidae)

Background

Solitary bees, such as the red mason bee (Osmia bicornis), provide important ecosystem services including pollination. In the face of global declines of pollinator abundance, such haplodiploid Hymenopterans have a compounded extinction risk due to the potential for limited genetic diversity. In order to assess the genetic diversity of Osmia bicornis populations, we developed microsatellite markers and characterised them in two populations.

Methods and results

Microsatellite sequences were mined from the recently published Osmia bicornis genome, which was assembled from DNA extracted from a single male bee originating from the United Kingdom. Sequences were identified that contained dinucleotide, trinucleotide, and tetranucleotide repeat regions. Seventeen polymorphic microsatellite markers were designed and tested, sixteen of which were developed into four multiplex PCR sets to facilitate cheap, fast and efficient genotyping and were characterised in unrelated females from Germany (n = 19) and England (n = 14).

Conclusions

The microsatellite markers are highly informative, with a combined exclusion probability of 0.997 (first parent), which will enable studies of genetic structure and diversity to inform conservation efforts in this bee.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-021-06796-x.

A Semi-Automated SNP-Based Approach for Contaminant Identification in Biparental Polyploid Populations of Tropical Forage Grasses

Artificial hybridization plays a fundamental role in plant breeding programs since it generates new genotypic combinations that can result in desirable phenotypes. Depending on the species and mode of reproduction, controlled crosses may be challenging, and contaminating individuals can be introduced accidentally. In this context, the identification of such contaminants is important to avoid compromising further selection cycles, as well as genetic and genomic studies. The main objective of this work was to propose an automated multivariate methodology for the detection and classification of putative contaminants, including apomictic clones (ACs), self-fertilized individuals, half-siblings (HSs), and full contaminants (FCs), in biparental polyploid progenies of tropical forage grasses. We established a pipeline to identify contaminants in genotyping-by-sequencing (GBS) data encoded as allele dosages of single nucleotide polymorphism (SNP) markers by integrating principal component analysis (PCA), genotypic analysis (GA) measures based on Mendelian segregation, and clustering analysis (CA). The combination of these methods allowed for the correct identification of all contaminants in all simulated progenies and the detection of putative contaminants in three real progenies of tropical forage grasses, providing an easy and promising methodology for the identification of contaminants in biparental progenies of tetraploid and hexaploid species. The proposed pipeline was made available through the polyCID Shiny app and can be easily coupled with traditional genetic approaches, such as linkage map construction, thereby increasing the efficiency of breeding programs.

A 2b-RAD parentage analysis pipeline for complex and mixed DNA samples

Next-generation sequencing technology has revolutionized genotyping in many fields of study, yet parentage analysis often still relies on microsatellite markers that are costly to generate and are currently available only for a limited number of species. 2b-RAD sequencing (2b-RAD) is a DNA sequencing technique developed for ecological population genomics that utilizes type IIB restriction enzymes to generate consistent, uniform fragments across samples. This technology is inexpensive, effective with low DNA inputs, and robust to DNA degradation. Here, we developed a probabilistic genotyping-by-sequencing genetic testing pipeline for parentage analysis by using 2b-RAD for inferring familial relationships from mixed DNA samples and populations. Our approach to partial paternity assignment utilizes a novel weighted outlier paternity index (WOPI) adapted for next-generation sequencing data and an identity-by-state (IBS) matrix-based clustering method for pedigree reconstruction. The combination of these two parentage assignment methods overcomes two major obstacles faced by other genetic testing methods: 1) It allows detection of parentage when closely related or inbred individuals are in the alleged parent population (e.g., in laboratory strains); and 2) it resolves mixed DNA samples. We successfully demonstrate this novel approach by correctly inferring paternity for samples pooled from multiple offspring (i.e., entire clutches) in a highly inbred population of an East African cichlid fish. The unique advantages of 2b-RAD in combination with our bioinformatics pipeline enable straightforward and cost-effective parentage analysis in any species regardless of genomic resources available.

Non-linear genetic diversity and notable population differentiation caused by low gene flow of bermudagrass [Cynodon dactylon (L.) Pers.] along longitude gradients

Background

Environmental variation related to ecological habitat is the main driver of plant adaptive divergence. Longitude plays an important role in the formation of plant population structure, indicating that environmental differentiation can significantly shape population structure.

Methods

Genetic diversity and population genetic structure were estimated using 105 expressed sequence tag-derived simple sequence repeat (EST-SSR) loci. A total of 249 C. dactylon (L.) Pers. (common bermudagrass) individuals were sampled from 13 geographic sites along the longitude (105°57′34″–119°27′06″E).

Results

There was no obvious linear trend of intra-population genetic diversity along longitude and the intra-population genetic diversity was not related to climate in this study. Low gene flow (Nm = 0.7701) meant a rich genetic differentiation among populations of C. dactylon along longitude gradients. Significantly positive Mantel correlation (r = 0.438, P = 0.001) was found between genetic distance and geographical interval while no significant partial Mantel correlation after controlling the effect of mean annual precipitation, which indicated geographic distance correlated with mean annual precipitation affect genetic distance. The genetic diversity of C. dactylon with higher ploidy level was higher than that with lower ploidy level and groups of individuals with higher ploidy level were separated further away by genetic distance from the lower ploidy levels. Understanding the different genetic bases of local adaptation comparatively between latitude and longitude is one of the core findings in the adaptive evolution of plants.

Genetic Diversity and Population Structure of Metaphire vulgaris Based on the Mitochondrial COI Gene and Microsatellites

The earthworm species Metaphire vulgaris (a member of the Clitellata class) is widely distributed across China, and has important ecological functions and medicinal value. However, investigations into its genetic diversity and differentiation are scarce. Consequently, we evaluated the genetic diversity of five populations of M. vulgaris (GM, HD, NYYZ, QDDY, and QDY) in Yancheng, China via the mitochondrial COI gene and the novel microsatellites developed there. A total of nine haplotypes were obtained by sequencing the mitochondrial COI gene, among which NYYZ and QDDY populations had the greatest number of haplotypes (nh = 5). Further, the nucleotide diversity ranged from 0.00437 to 0.1243. The neighbor-joining trees and the TCS network of haplotypes indicated that earthworm populations within close geographical range were not genetically isolated at these small scale distances. Results of the identification of microsatellite molecular markers revealed that the allele number in 12 microsatellite loci ranged from 4 to 13. The observed heterozygosity ranged from 0.151 to 0.644, whereas the expected heterozygosity ranged from 0.213 to 0.847. The polymorphism data content of most sites was >0.5, which indicated that the designed sites had high polymorphism. Structural analysis results indicated that GM, HD, and NYYZ had similar genetic structures across the five populations. The Nei's genetic distance between HD and NYYZ populations was the smallest (D s = 0.0624), whereas that between HD and QDY populations was the largest (D s = 0.2364). The UPGMA tree showed that HD were initially grouped with NYYZ, followed by GM, and then with QDDY. Furthermore, cross-species amplification tests were conducted for Metaphire guillelmi, which indicated that the presented markers were usable for this species. This study comprised a preliminary study on the genetic diversity of M. vulgaris, which provides basic data for future investigations into this species.

STRs: Ancient Architectures of the Genome beyond the Sequence

Short tandem repeats (STRs) are commonly defined as short runs of repetitive nucleotides, consisting of tandemly repeating 2-6- bp motif units, which are ubiquitously distributed throughout genomes. Functional STRs are polymorphic in the population, and their variations influence gene expression, which subsequently may result in pathogenic phenotypes. To understand STR phenotypic effects and their functional roles, we describe four different mutational mechanisms including the unequal crossing-over model, gene conversion, retrotransposition mechanism and replication slippage. Due to the multi-allelic nature, small length, abundance, high variability, codominant inheritance, nearly neutral evolution, extensive genome coverage and simple assaying of STRs, these markers are widely used in various types of biological research, including population genetics studies, genome mapping, molecular epidemiology, paternity analysis and gene flow studies. In this review, we focus on the current knowledge regarding STR genomic distribution, function, mutation and applications.

Microsatellites as Agents of Adaptive Change: An RNA-Seq-Based Comparative Study of Transcriptomes from Five Helianthus Species

Mutations that provide environment-dependent selective advantages drive adaptive divergence among species. Many phenotypic differences among related species are more likely to result from gene expression divergence rather than from non-synonymous mutations. In this regard, cis-regulatory mutations play an important part in generating functionally significant variation. Some proposed mechanisms that explore the role of cis-regulatory mutations in gene expression divergence involve microsatellites. Microsatellites exhibit high mutation rates achieved through symmetric or asymmetric mutation processes and are abundant in both coding and non-coding regions in positions that could influence gene function and products. Here we tested the hypothesis that microsatellites contribute to gene expression divergence among species with 50 individuals from five closely related Helianthus species using an RNA-seq approach. Differential expression analyses of the transcriptomes revealed that genes containing microsatellites in non-coding regions (UTRs and introns) are more likely to be differentially expressed among species when compared to genes with microsatellites in the coding regions and transcripts lacking microsatellites. We detected a greater proportion of shared microsatellites in 5′UTRs and coding regions compared to 3′UTRs and non-coding transcripts among Helianthus spp. Furthermore, allele frequency differences measured by pairwise FST at single nucleotide polymorphisms (SNPs), indicate greater genetic divergence in transcripts containing microsatellites compared to those lacking microsatellites. A gene ontology (GO) analysis revealed that microsatellite-containing differentially expressed genes are significantly enriched for GO terms associated with regulation of transcription and transcription factor activity. Collectively, our study provides compelling evidence to support the role of microsatellites in gene expression divergence.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers

Background

Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM.

Results

This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (p_GSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM.

Conclusions

This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01805-w.

Contrasting effects of local environment and grazing pressure on the genetic diversity and structure of Artemisia frigida

Drylands count among the most globally extensive biomes, and while many desert and dry rangeland ecosystems are under threat, genetic structures of dryland species are still rarely studied. Artemisia frigida is one of the most widely distributed plant species in the temperate rangelands of Eurasia and North America, and it also dominates in many habitats of Mongolia due to its tolerance to low temperatures, drought and disturbance. Local environmental conditions and grazing pressure can influence species performance and affect spatial patterns of genetic diversity in contrasting ways, and our study set out to evaluate such effects on the genetic diversity and structure of A. frigida. We first developed new species-specific Simple Sequence Repeats (SSRs) markers using whole genome sequencing. We then analysed 11 populations of A. frigida that had been sampled along a large climatic gradient in Mongolia, which were sub-structured according to three levels of grazing intensity. Estimates of genetic diversity at the population level were high (HO = 0.56, HE = 0.73) and tended to increase with higher precipitation and soil nutrient availability. Grazing had no effect on genetic diversity, however, a high number of grazing-specific indicator alleles was found at grazed sites. Genetic differentiation among populations was extremely low (global GST = 0.034). Analysis of Molecular Variance revealed 5% variance between populations along the climatic gradient, with 3% of the variance being partitioned among different grazing intensity levels. We found no relationship between geographic and genetic distances, and thus no isolation by distance in this widely distributed species. The relatively low genetic structuring suggests that considerable gene flow exists among A. frigida populations across the rangelands of Mongolia, in spite of the pervasive grazing in the region.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Comparing Genetic Diversity in Three Threatened Oaks

Genetic diversity is a critical resource for species’ survival during times of environmental change. Conserving and sustainably managing genetic diversity requires understanding the distribution and amount of genetic diversity (in situ and ex situ) across multiple species. This paper focuses on three emblematic and IUCN Red List threatened oaks (Quercus, Fagaceae), a highly speciose tree genus that contains numerous rare species and poses challenges for ex situ conservation. We compare the genetic diversity of three rare oak species—Quercus georgiana, Q. oglethorpensis, and Q. boyntonii—to common oaks; investigate the correlation of range size, population size, and the abiotic environment with genetic diversity within and among populations in situ; and test how well genetic diversity preserved in botanic gardens correlates with geographic range size. Our main findings are: (1) these three rare species generally have lower genetic diversity than more abundant oaks; (2) in some cases, small population size and geographic range correlate with genetic diversity and differentiation; and (3) genetic diversity currently protected in botanic gardens is inadequately predicted by geographic range size and number of samples preserved, suggesting non-random sampling of populations for conservation collections. Our results highlight that most populations of these three rare oaks have managed to avoid severe genetic erosion, but their small size will likely necessitate genetic management going forward.

Impact of different numbers of microsatellite markers on population genetic results using SLAF-seq data for Rhododendron species

Microsatellites (simple sequence repeats, SSRs) are co-dominant nuclear markers that are widely used in population genetic studies. Population genetic parameters from different studies might be significantly influenced by differences in marker number. In our study, 265 sequences with polymorphic microsatellites were obtained from SLAF-seq data. Then, subpopulations containing different numbers (5, 6, 7,…, 15, 20, 25, 30, 35, 40) of markers were genotyped 10 times to investigate the impact of marker numbers on population genetic diversity results. Our results show that genotyping with less than 11 or 12 microsatellite markers lead to significant deviations in the population genetic diversity or genetic structure results. In order to provide markers for population genetic and conservation studies for Rhododendron, 26 SSR primers were designed and validated in three species.

Stevia rebaudiana germplasm characterization using microsatellite markers and steviol glycosides quantification by HPLC

Stevia rebaudiana Bertoni, Asteraceae, is an herbaceous perennial plant native to Paraguay. This species is considered since ancient times a medicinal plant with important bioactive compounds and pharmacologic and food properties, namely diterpenes glycosides. The high natural sweetener potential stevioside and rebaudioside A produced by S. rebaudiana plants are suitable sucrose substitutes, and their obtention is influenced by environmental, phytosociological, and genetic factors. The plants’ genetic profile and sweet potential depiction are needed for suitable plant selection for improvement and deployment. Thirty-one S. rebaudiana accessions grown in the same plot where leaves samples were collected in early 2019, were genotyped using six microsatellite markers, including two steviol glycosides biosynthesis functionally involved markers. Additionally, an aqueous extract of each sample was obtained in a water bath and purified by SPE for stevioside and rebaudioside A quantification by normal phase HPLC. Stevioside and rebaudioside A contents varied between 0.53–7.36% (w w⁻¹) and 0.37–3.60% (w w⁻¹), respectively. Two genotypes displayed interesting ratios of rebaudioside A/stevioside (number 3 and 33). The level of genetic similarity between genotypes was tested through a pairwise similarity coefficient, and two groups of individuals had the same fingerprinting. Strong relatedness was found within genotypes, possibly due to cloning, thus, influx of new germplasm ought to be made to prevent mating between relatives, and for further selection and genetic improvement.

Performance comparison of two reduced-representation based genome-wide marker-discovery strategies in a multi-taxon phylogeographic framework

Multi-locus genetic data are pivotal in phylogenetics. Today, high-throughput sequencing (HTS) allows scientists to generate an unprecedented amount of such data from any organism. However, HTS is resource intense and may not be accessible to wide parts of the scientific community. In phylogeography, the use of HTS has concentrated on a few taxonomic groups, and the amount of data used to resolve a phylogeographic pattern often seems arbitrary. We explore the performance of two genetic marker sampling strategies and the effect of marker quantity in a comparative phylogeographic framework focusing on six species (arthropods and plants). The same analyses were applied to data inferred from amplified fragment length polymorphism fingerprinting (AFLP), a cheap, non-HTS based technique that is able to straightforwardly produce several hundred markers, and from restriction site associated DNA sequencing (RADseq), a more expensive, HTS-based technique that produces thousands of single nucleotide polymorphisms. We show that in four of six study species, AFLP leads to results comparable with those of RADseq. While we do not aim to contest the advantages of HTS techniques, we also show that AFLP is a robust technique to delimit evolutionary entities in both plants and animals. The demonstrated similarity of results from the two techniques also strengthens biological conclusions that were based on AFLP data in the past, an important finding given the wide utilization of AFLP over the last decades. We emphasize that whenever the delimitation of evolutionary entities is the central goal, as it is in many fields of biodiversity research, AFLP is still an adequate technique.

A new one-tube reaction kit for the SSR genotyping of apple (Malus × domestica Borkh.)

Though apple genotyping is mainly used for scientific and breeding purposes, it can also be adopted by national authorities to control the authenticity of apple cultivars. To facilitate the introduction of routine apple genotyping into practice, a new apple simple sequence repeat (SSR) genotyping kit was developed (called the Ap17 in. SSR Genotyping Kit). The kit combines 17 SSR markers including those recommended by the Working Group of the European Cooperative Programme for Plant Genetic Resources (ECPGR), covering all apple linkage groups in a one-tube reaction format, using a fragment analysis method to simplify the genotyping procedure. The kit was successfully tested using 880 unique diploid apple germplasm accessions; the kit can also readily discriminate triploid and tetraploid samples. The total probability of identity for the kit and the sample collection used was calculated to be 1.73 × 10^-22. Tables for converting results to enable genotype comparisons between currently-used genotyping systems and the Ap17 in. kit are provided. The kit is ideally suited for validation in laboratories genotyping a large number of apple samples, saving time, costs, and labor, while minimizing technical and human errors.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

Population genetics and host specificity of Varroa destructor mites infesting eastern and western honeybees

In a globalized world, parasites are often brought in contact with new potential hosts. When parasites successfully shift host, severe diseases can emerge at a large cost to society. However, the evolutionary processes leading to successful shifts are rarely understood, hindering risk assessment, prevention, or mitigation of their effects. Here, we screened populations of Varroa destructor, an ectoparasitic mite of the honeybee genus Apis, to investigate their genetic structure and reproductive potential on new and original hosts. From the patterns identified, we deduce the factors that influenced the macro- and microevolutionary processes that led to the structure observed. Among the mite variants identified, we found two genetically similar populations that differed in their reproductive abilities and thus in their host specificity. These lineages could interbreed, which represents a threat due to the possible increased virulence of the parasite on its original host. However, interbreeding was unidirectional from the host-shifted to the nonshifted native mites and could thus lead to speciation of the former. The results improve our understanding of the processes affecting the population structure and evolution of this economically important mite genus and suggest that introgression between shifted and nonshifted lineages may endanger the original host.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10340-020-01322-7.

Comparative analysis of genomic- and EST-SSRs in European plum (Prunus domestica L.): implications for the diversity analysis of polyploids

Simple sequence repeats (SSRs) are among the most useful DNA markers in plant science. The aim of this study was to compare the features and usefulness of genomic SSR (gSSR) and EST-SSRs in European plum (Prunus domestica L.), an economically important, hexaploid stone fruit crop globally cultivated to produce fleshy fruits and derived foodstuff. The analysis of an ample set of morphologically diverse varieties indicated that gSSRs and EST-SSRs provide different estimates of some of the locus-based indicators of diversity. Moreover, the two classes of SSRs gave different, weakly correlated, estimations of distance-based parameters with gSSRs being more powerful for discriminating purposes. The two SSR classes provide complementary information in European plum, making the contribution of EST-SSRs useful not only as non-neutral markers. The differences between SSR classes are discussed considering the neutral and non-neutral evolution, and the polyploidy and asexual propagation of the cultivated tree varieties.

SSRgenotyper: A simple sequence repeat genotyping application for whole-genome resequencing and reduced representational sequencing projects

PREMISE

Many programs can identify simple sequence repeat (SSR) motifs in genomic data. SSRgenotyper extends SSR identification to en masse genotyping from resequencing data for diversity panels and linkage mapping populations.

METHODS AND RESULTS

SSRgenotyper will find and genotype SSRs from SAM files and an SSR reference FASTA. Several outputs are possible, including a simple table with the SSR marker name, position, and SSR alleles, defined by the repeat number of the repeat motif. Specific output files include a GENEPOP-formatted file for downstream genetic diversity analyses and a traditional A, H, B mapping file output that is phased to the parents of the population for biparental linkage map construction. Linkage maps produced using SSRgenotyper genotypes were highly collinear with physical maps and correctly inferred known phylogenies.

CONCLUSIONS

SSRgenotyper provides an easy-to-use, accurate, and scalable SSR genotyping platform for whole-genome resequencing data. SSRgenotyper is freely available at https://github.com/dlewis27/SSRgenotyper.

A reference allelic ladder for Western Capercaillie (Tetrao urogallus) and Black Grouse (Tetrao tetrix) enables linking grouse genetic data across studies

Rapid anthropogenic climate change and progressing habitat degradation are considered top threats to biodiversity. The employment of demanding umbrella species as indicators for ecosystem health is a popular and cost-effective strategy that facilitates continuous monitoring and evaluation within a long-term conservation management scheme. The Western Capercaillie (Tetrao urogallus) and the Black Grouse (Tetrao tetrix) are both considered viable candidates due to their extensive habitat requirements, the possibility for conservative, non-invasive sampling, and their broad popular appeal. Regional population surveys based on genetic data from Short Sequence Repeat (SSR) analysis are being conducted throughout the Palearctic. However, to ensure reliable comparability among laboratories, standardization is required. Here, we report a catalogue of fifty fully characterized reference alleles from twelve SSR loci and the construction of a customizable allelic ladder for genotyping and individualization in Western Capercaillie and Black Grouse. This methodological improvement will help to cost-efficiently generate and collate supraregional data from different grouse surveys and thereby contribute to conservation management. Reference alleles and ladders can be obtained on demand.

Population Genetics of Polyploid Complex Psidium cattleyanum Sabine (Myrtaceae): Preliminary Analyses Based on New Species-Specific Microsatellite Loci and Extension to Other Species of the Genus

Polyploidy is a phenomenon that alters the genetic diversity of populations and has been reported as one of the most important evolutionary forces for plant diversification. The Psidium cattleyanum complex comprises a group of wild populations with several ploidy levels reported in the literature. The multiple cytotypes, associated with its wide distribution area, make this species a potential key model for understanding evolutionary processes related to polyploidization. In this study, we isolated and characterized nuclear microsatellite markers of P. cattleyanum and tested their transferability to other nine species of the genus. We performed a preliminary analysis of genetic diversity and population structure in three populations of P. cattleyanum. The three populations analyzed had different chromosome numbers, being polyploid cytotypes (2n = 6x = 66, 2n = 7x = 77 and 2n = 8x = 88). We designed 46 primer pairs and successfully amplified 37 markers, from which the 10 best were selected for analysis. Considering both the PIC and DP values, most of markers were highly informative. The new SSR markers were used to assess the levels of genetic diversity of the populations and detected one population with predominance of sexual reproduction. DAPC analysis pointed the formation of three groups, which corresponded to the populations analyzed. The markers were successfully amplified in related species, with some species presenting 80% transferability. By producing this panel of polymorphic microsatellites, we contribute to the understanding evolution in groups of natural polyploids for future studies.

Comparative assessment of range-wide patterns of genetic diversity and structure with SNPs and microsatellites: A case study with Iberian amphibians

Reduced representation genome sequencing has popularized the application of single nucleotide polymorphisms (SNPs) to address evolutionary and conservation questions in nonmodel organisms. Patterns of genetic structure and diversity based on SNPs often diverge from those obtained with microsatellites to different degrees, but few studies have explicitly compared their performance under similar sampling regimes in a shared analytical framework. We compared range‐wide patterns of genetic structure and diversity in two amphibians endemic to the Iberian Peninsula: Hyla molleri and Pelobates cultripes, based on microsatellite (18 and 14 loci) and SNP (15,412 and 33,140 loci) datasets of comparable sample size and spatial extent. Model‐based clustering analyses with STRUCTURE revealed minor differences in genetic structure between marker types, but inconsistent values of the optimal number of populations (K) inferred. SNPs yielded more repeatable and less admixed ancestries with increasing K compared to microsatellites. Genetic diversity was weakly correlated between marker types, with SNPs providing a better representation of southern refugia and of gradients of genetic diversity congruent with the demographic history of both species. Our results suggest that the larger number of loci in a SNP dataset can provide more reliable inferences of patterns of genetic structure and diversity than a typical microsatellite dataset, at least at the spatial and temporal scales investigated.

Conservation Genomics of a Threatened Rhododendron: Contrasting Patterns of Population Structure Revealed From Neutral and Selected SNPs

Though it is well-acknowledged that next generation sequencing (NGS) technologies can provide further insights into plant conservation management than traditional molecular markers, studies employing NGS to address conservation genomics and subsequent conservation strategies for threatened plants are still rare. Rhododendron is the largest genus of woody plants in China, and many species are threatened, however, to date there has been no conservation genetic research using NGS in this genus. In the present study, we investigated the conservation genetics of R. cyanocarpum, a threatened species endemic to the Cangshan Mountains in Yunnan, China, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. Due to the availability of sufficient SNPs, we were able to distinguish between neutral and putatively selected SNPs and were able to further investigate the genetic diversity, population structure, and differentiation in R. cyanocarpum, as well as make an estimation of its demographic history. A total of 6,584 SNPs were obtained, of which 5,729 were neutral (detected using Tajima’s D). In terms of the 5,729 neutral SNPs, R. cyanocarpum had a higher genetic diversity (π = 0.0702 ± 0.0017, H_e = 0.0675 ± 0.0016) than other plant species assessed using Rad-seq methods, while population differentiation (F_st from 0.0314 to 0.0452) was weak. Interestingly, contrasting patterns of population structure were revealed from all neutral and selected SNPs, with distinct genetic clusters forming for all SNPs and neutral SNPs, but no distinct subgroups for selected ones. Moreover, we were able to detect changes in effective population size (N_e) of R. cyanocarpum from 150,000 years ago, including a bottleneck event ca. 60,000 years ago, followed by recovery of N_e over a short period, and a subsequent gradual decline in N_e to date. Implications for conserving R. cyanocarpum based on these main results are then discussed.

Development of novel microsatellite markers using RAD sequencing technology for diversity assessment of rambutan (Nephelium lappaceum L.) germplasm

The trend of microsatellite marker discovery and development revolved as a result of the advancement of next generation sequencing (NGS) technology as it has developed numerous microsatellites within a short period of time at a low cost. This study generated microsatellite markers using RAD sequencing technologies for the understudied Nephelium lappaceum. A total of 1403 microsatellite markers were successfully designed, which consisted of 853 di-, 525 tri-, 17 tetra-, 5 penta-, and 3 hexanucleotide microsatellite markers. Subsequently, selection of 39 microsatellites was made for the evaluation of genetic diversity of the selected 22 rambutan varieties. Twelve microsatellites, which exhibited high call rates across the samples, were used to assess the diversity of the aforementioned rambutan varieties. The analysis of 12 microsatellites revealed the presence of 72 alleles and six alleles per locus in average. Furthermore, the polymorphic information content (PIC) value ranged from 0.326 (NlaSSR20) to 0.832 (NlaSSR32), which included an average of 0.629 per locus, while the generated Neighbour Joining dendrogram showed two major clusters. The pairwise genetic distance of shared alleles exhibited a range of values from 0.046 (R134↔R170) to 0.818 (R5↔R170), which suggested highest dissimilarity detected between R5 and R170. Notably, these research findings would useful for varietal identification, proper management and conservation of the genetic resources, and exploitation and utilization in future breeding programs.

Population size estimates based on the frequency of genetically assigned parent-offspring pairs within a subsample

Estimating population density as precise as possible is a key premise for managing wild animal species. This can be a challenging task if the species in question is elusive or, due to high quantities, hard to count. We present a new, mathematically derived estimator for population size, where the estimation is based solely on the frequency of genetically assigned parent-offspring pairs within a subsample of an ungulate population. By use of molecular markers like microsatellites, the number of these parent-offspring pairs can be determined. The study's aim was to clarify whether a classical capture-mark-recapture (CMR) method can be adapted or extended by this genetic element to a genetic-based capture-mark-recapture (g-CMR). We numerically validate the presented estimator (and corresponding variance estimates) and provide the R-code for the computation of estimates of population size including confidence intervals. The presented method provides a new framework to precisely estimate population size based on the genetic analysis of a one-time subsample. This is especially of value where traditional CMR methods or other DNA-based (fecal or hair) capture-recapture methods fail or are too difficult to apply. The DNA source used is basically irrelevant, but in the present case the sampling of an annual hunting bag is to serve as data basis. In addition to the high quality of muscle tissue samples, hunting bags provide additional and essential information for wildlife management practices, such as age, weight, or sex. In cases where a g-CMR method is ecologically and hunting-wise appropriate, it enables a wide applicability, also through its species-independent use.

Development of a Large Gene-Associated SSR Marker Set and in-Depth Genetic Characterization in Scarlet Sage

Salvia splendens, scarlet or tropical sage, is a tender perennial herbaceous flowering plant popularly grown in public and private gardens all over the world. In this study, we developed a set of simple sequence repeats (SSRs) from genome-wide sequences to assess the genetic diversity and population structure among 112 cultivars. We obtained 364,379 SSRs by mining scarlet sage’s recently published whole genome sequence; 14,545 gene-associated SSR loci were identified in 2 kb gene flanking regions. Among the 768 gene-associated SSR primer sets we screened, 576 loci successfully amplified in DNA pools of 3–4 different cultivars, of which 271 remained polymorphic when tested across eight individual plants. We searched for the related gene functions attributable to these gene-associated SSRs using diverse databases, resulting in 259 Non-redundant matching sequences, 205 individual Gene Ontology (GO) terms, 236 assigned to eukaryotic orthologous groups, and 67 KEGG-annotated (Kyoto Encyclopedia of Genes and Genomes) sequences. We finally selected 41 polymorphic SSR loci to infer genetic diversity and population structure among 112 S. splendens accessions. Based on the developed gene-associated SSRs, clustering analyses consistently revealed two distinct genetic groups within the core collection of S. splendens cultivars. This work developed and characterized an exhaustive set of genome-wide gene-associated SSR markers for scarlet sage. These SSRs can provide species identification, genetic diversity and population structure information for S. splendens, and will therefore be important tools for the management and protection of S. splendens germplasm.

Transcribed microsatellite allele lengths are often correlated with gene expression in natural sunflower populations

Microsatellites are common in genomes of most eukaryotic species. Due to their high mutability, an adaptive role for microsatellites has been considered. However, little is known concerning the contribution of microsatellites towards phenotypic variation. We used populations of the common sunflower (Helianthus annuus) at two latitudes to quantify the effect of microsatellite allele length on phenotype at the level of gene expression. We conducted a common garden experiment with seed collected from sunflower populations in Kansas and Oklahoma followed by an RNA-Seq experiment on 95 individuals. The effect of microsatellite allele length on gene expression was assessed across 3,325 microsatellites that could be consistently scored. Our study revealed 479 microsatellites at which allele length significantly correlates with gene expression (eSTRs). When irregular allele sizes not conforming to the motif length were removed, the number of eSTRs rose to 2,379. The percentage of variation in gene expression explained by eSTRs ranged from 1%-86% when controlling for population and allele-by-population interaction effects at the 479 eSTRs. Of these eSTRs, 70.4% are in untranslated regions (UTRs). A gene ontology (GO) analysis revealed that eSTRs are significantly enriched for GO terms associated with cis- and trans-regulatory processes. Our findings suggest that a substantial number of transcribed microsatellites can influence gene expression.