Always look on both sides: phylogenetic information conveyed by simple sequence repeat allele sequences

The sequence of the repetitive motif influences the frequency of multistep mutations in Short Tandem Repeats

Microsatellites, or Short Tandem Repeats (STRs), are subject to frequent length mutations that involve the loss or gain of an integer number of repeats. This work aimed to investigate the correlation between STRs' specific repetitive motif composition and mutational dynamics, specifically the occurrence of single- or multistep mutations. Allelic transmission data, comprising 323,818 allele transfers and 1,297 mutations, were gathered for 35 Y-chromosomal STRs with simple structure. Six structure groups were established: ATT, CTT, TCTA/GATA, GAAA/CTTT, CTTTT, and AGAGAT, according to the repetitive motif present in the DNA leading strand of the markers. Results show that the occurrence of multistep mutations varies significantly among groups of markers defined by the repetitive motif. The group of markers with the highest frequency of multistep mutations was the one with repetitive motif CTTTT (25% of the detected mutations) and the lowest frequency corresponding to the group with repetitive motifs TCTA/GATA (0.93%). Statistically significant differences (α = 0.05) were found between groups with repetitive motifs with different lengths, as is the case of TCTA/GATA and ATT (p = 0.0168), CTT (p < 0.0001) and CTTTT (p < 0.0001), as well as between GAAA/CTTT and CTTTT (p = 0.0102). The same occurred between the two tetrameric groups GAAA/CTTT and TCTA/GATA (p < 0.0001) - the first showing 5.7 times more multistep mutations than the second. When considering the number of repeats of the mutated paternal alleles, statistically significant differences were found for alleles with 10 or 12 repeats, between GATA and ATT structure groups. These results, which demonstrate the heterogeneity of mutational dynamics across repeat motifs, have implications in the fields of population genetics, epidemiology, or phylogeography, and whenever STR mutation models are used in evolutionary studies in general.

Joint analysis of microsatellites and flanking sequences enlightens complex demographic history of interspecific gene flow and vicariance in rear-edge oak populations

Inference of recent population divergence requires fast evolving markers and necessitates to differentiate shared genetic variation caused by ancestral polymorphism and gene flow. Theoretical research shows that the use of compound marker systems integrating linked polymorphisms with different mutational dynamics, such as a microsatellite and its flanking sequences, can improve estimation of population structure and inference of demographic history, especially in the case of complex population dynamics. However, empirical application in natural populations has so far been limited by lack of suitable methods for data collection. A solution comes from the development of sequence-based microsatellite genotyping which we used to study molecular variation at 36 sequenced nuclear microsatellites in seven Quercus canariensis and four Q. faginea rear-edge populations across Algeria. We aim to decipher their taxonomic relationship, past evolutionary history and recent demographic trajectory. First, we compare the estimation of population genetics parameters and simulation-based inference of demographic history from microsatellite sequence alone, flanking sequence alone or the combination of linked microsatellite and flanking sequence variation. Second, we apply random forest approximate Bayesian computation to identify which of these sequence types is most informative. Whereas analysing microsatellite variation alone indicates recent interspecific gene flow, additional information gained by integrating nucleotide variation in flanking sequences, by reducing homoplasy, suggests ancient interspecific gene flow followed by drift in isolation instead. The weight of each polymorphism in the inference also demonstrates the value of linked variations with contrasted mutation dynamic to improve estimation of both demographic and mutational parameters.

High genetic and chemical diversity of wild hop populations from Central Italy with signals of a genetic structure influenced by both sexual and asexual reproduction

In order to investigate the intraspecific diversity of wild Humulus lupulus (hop) in Central Italy, 12 populations were evaluated for their genetic polymorphism by means of 13 SSR loci together with six commercial cultivars as a reference. High levels of polymorphism were found across the populations, being 140 the number of multilocus genotypes over 159 samples analyzed. Moreover, the observed heterozygosity was higher than expected in most of the populations. High levels of gene flow were thus envisaged to occur within and among wild populations, and our sampling strategy allowed us to gain insights on the propagation modes of this species, i.e. clonal versus sexual propagation. Nevertheless, a genetic structure of populations with at least five genetically different clusters was disclosed. Private alleles were observed in both wild and cultivated hops. Chemical analysis of bittering and aromatic quality of female flowers from a subset of 8 wild populations revealed a high variability among plants, especially for essential oil components. Overall, the high variability of wild accessions here examined represent a valid source to be exploited in future breeding programs for new or improved hop cultivars development.

Fast sequence-based microsatellite genotyping development workflow

Application of high-throughput sequencing technologies to microsatellite genotyping (SSRseq) has been shown to remove many of the limitations of electrophoresis-based methods and to refine inference of population genetic diversity and structure. We present here a streamlined SSRseq development workflow that includes microsatellite development, multiplexed marker amplification and sequencing, and automated bioinformatics data analysis. We illustrate its application to five groups of species across phyla (fungi, plant, insect and fish) with different levels of genomic resource availability. We found that relying on previously developed microsatellite assay is not optimal and leads to a resulting low number of reliable locus being genotyped. In contrast, de novo ad hoc primer designs gives highly multiplexed microsatellite assays that can be sequenced to produce high quality genotypes for 20-40 loci. We highlight critical upfront development factors to consider for effective SSRseq setup in a wide range of situations. Sequence analysis accounting for all linked polymorphisms along the sequence quickly generates a powerful multi-allelic haplotype-based genotypic dataset, calling to new theoretical and analytical frameworks to extract more information from multi-nucleotide polymorphism marker systems.

Advances in genotyping microsatellite markers through sequencing and consequences of scoring methods for Ceratonia siliqua (Leguminosae)

PREMISE OF THE STUDY

Simple sequence repeat (SSR) or microsatellite markers have been used in a broad range of studies mostly scoring alleles on the basis of amplicon size as a proxy for the number of repeat units of an SSR motif. However, additional sources of variation within the SSR or in the flanking regions have largely remained undetected.

METHODS

In this study, we implemented a next-generation sequencing-based genotyping approach in a newly characterized set of 18 nuclear SSR markers for the carob tree, Ceratonia siliqua. Our aim was to evaluate the effect of three different methods of scoring molecular variation present within microsatellite markers on the genetic diversity and structure results.

RESULTS

The analysis of the sequences of 77 multilocus genotypes from four populations revealed SSR variation and additional sources of polymorphism in 87% of the loci analyzed (42 single-nucleotide polymorphisms and five insertion/deletion polymorphisms), as well as divergent paralog copies in two loci. Ignoring sequence variation under standard amplicon size genotyping resulted in incorrect identification of 69% of the alleles, with important effects on the genetic diversity and structure estimates.

DISCUSSION

Next-generation sequencing allows the detection and scoring of SSRs, single-nucleotide polymorphisms, and insertion/deletion polymorphisms to increase the resolution of population genetic studies.

SSR-seq: Genotyping of microsatellites using next-generation sequencing reveals higher level of polymorphism as compared to traditional fragment size scoring

Microsatellites (or simple sequence repeats, SSR) are widely used markers in population genetics. Traditionally, genotyping was and still is carried out through recording fragment length. Now, next-generation sequencing (NGS) makes it easy to obtain also sequence information for the loci of interest. This avoids misinterpretations that otherwise could arise due to size homoplasy. Here, an NGS strategy is described that allows to genotype hundreds of individuals at many custom-designed SSR loci simultaneously, combining multiplex PCR, barcoding, and Illumina sequencing. We created three different datasets for which alleles were coded according to (a) length of the repetitive region, (b) total fragment length, and (c) sequence identity, in order to evaluate the eventual benefits from having sequence data at hand, not only fragment length data. For each dataset, genetic diversity statistics, as well as F ST and R ST values, were calculated. The number of alleles per locus, as well as observed and expected heterozygosity, was highest in the sequence identity dataset, because of single-nucleotide polymorphisms and insertions/deletions in the flanking regions of the SSR motif. Size homoplasy was found to be very common, amounting to 44.7%-63.5% (mean over all loci) in the three study species. Thus, the information obtained by next-generation sequencing offers a better resolution than the traditional way of SSR genotyping and allows for more accurate evolutionary interpretations.

Are species-pairs diverging lineages? A nine-locus analysis uncovers speciation among species-pairs of the Lobaria meridionalis-group (Ascomycota)

In spite of considerable effort to verify the theory of species-pairs, uncertainty still exists about the relationship between sexually or vegetatively reproducing populations of morphologically indistinguishable, sympatric lichen species. The current paper studies putative species-pairs within the Asian Lobaria meridionalis-group, using a nine-locus and time calibrated species-tree approach. Analyses demonstrate that pairs of sexually or vegetatively reproducing lineages split into highly supported monophyletic clades-confirming molecularly the species-pair concept for the L. meridionalis-group. In the broader context of evolution and speciation dynamics in lichenized fungi, this paper attempts to synthesize molecular findings from the last two decades to promote a more modern perception of the species-pair concept. Taxonomically, eight species were found to currently conform to the L. meridionalis-group, which differentiated during the Pliocene and Pleistocene. The coincidence of paleoclimatic events with estimated dates of divergence support a bioclimatic hypothesis for the evolution of species in the L. meridionalis-group, which also explains their current eco-geographic distribution patterns. Greater recognition for species with a long and independent evolutionary history, which merit high conservation priority, will be especially critical for preserving geographically restricted endemics from Southeast Asia, where habitat loss is driving rapid declines.

Development of genomic tools in a widespread tropical tree, Symphonia globulifera L.f.: a new low-coverage draft genome, SNP and SSR markers

Population genetic studies in tropical plants are often challenging because of limited information on taxonomy, phylogenetic relationships and distribution ranges, scarce genomic information and logistic challenges in sampling. We describe a strategy to develop robust and widely applicable genetic markers based on a modest development of genomic resources in the ancient tropical tree species Symphonia globulifera L.f. (Clusiaceae), a keystone species in African and Neotropical rainforests. We provide the first low-coverage (11X) fragmented draft genome sequenced on an individual from Cameroon, covering 1.027 Gbp or 67.5% of the estimated genome size. Annotation of 565 scaffolds (7.57 Mbp) resulted in the prediction of 1046 putative genes (231 of them containing a complete open reading frame) and 1523 exact simple sequence repeats (SSRs, microsatellites). Aligning a published transcriptome of a French Guiana population against this draft genome produced 923 high-quality single nucleotide polymorphisms. We also preselected genic SSRs in silico that were conserved and polymorphic across a wide geographical range, thus reducing marker development tests on rare DNA samples. Of 23 SSRs tested, 19 amplified and 18 were successfully genotyped in four S. globulifera populations from South America (Brazil and French Guiana) and Africa (Cameroon and São Tomé island, F_ST  = 0.34). Most loci showed only population-specific deviations from Hardy-Weinberg proportions, pointing to local population effects (e.g. null alleles). The described genomic resources are valuable for evolutionary studies in Symphonia and for comparative studies in plants. The methods are especially interesting for widespread tropical or endangered taxa with limited DNA availability.

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.

A review of the prevalence, utility, and caveats of using chloroplast simple sequence repeats for studies of plant biology

Microsatellites occur in all plant genomes and provide useful markers for studies of genetic diversity and structure. Chloroplast microsatellites (cpSSRs) are frequently targeted because they are more easily isolated than nuclear microsatellites. Here, we quantified the frequency and uses of cpSSRs based on a literature review of over 400 studies published 1995-2013. These markers are an important and economical tool for plant biologists and continue to be used alongside modern genomics approaches to study genetic diversity and structure, evolutionary history, and hybridization in native and agricultural species. Studies using species-specific primers reported a greater number of polymorphic loci than those employing universal primers. A major disadvantage to cpSSRs is fragment size homoplasy; therefore, we documented its occurrence at several cpSSR loci within and between species of Acmispon (Fabaceae). Based on our empirical data set, we recommend targeted sequencing of a subset of samples combined with fragment genotyping as a cost-efficient, data-rich approach to the use of cpSSRs and as a test of homoplasy. The availability of genomic resources for plants aids in the development of primers for new study systems, thereby enhancing the utility of cpSSRs across plant biology.

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

Advancing technologies have facilitated the ever-widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity-based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration-specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context-dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non-human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high-throughput technologies in population genetics.

Isolation and characterization of microsatellite markers and analysis of genetic diversity in Chinese jujube (Ziziphus jujuba Mill.)

Chinese jujube (Ziziphus jujuba Mill, 2n = 2× = 24, Rhamnaceae) is an economically important Chinese native species. It has high nutritional value, and its medicinal properties have led to extensive use in traditional oriental medicine. The characterization of genotypes using molecular markers is important for genetic studies and plant breeding. However, few simple sequence repeat (SSR) markers are available for this species. In this study, 1,488 unique SSR clones were isolated from Z. jujuba 'Dongzao' using enriched genomic libraries coupled with a three-primer colony PCR screening strategy, yielding a high enrichment rate of 73.3%. Finally, 1,188 (80.87%) primer pairs were amplified successfully in the size expected for 'Dongzao'. A total of 350 primer pairs were further selected and evaluated for their ability to detect polymorphisms across a panel of six diverse cultivars; among these, 301 primer pairs detected polymorphisms, and the polymorphism information content (PIC) value across all loci ranged from 0.15 to 0.82, with an average of 0.52. An analysis of 76 major cultivars employed in Chinese jujube production using 31 primer pairs revealed comparatively high genetic diversity among these cultivars. Within-population differences among individuals accounted for 98.2% of the observed genetic variation. Neighbor-joining clustering divided the cultivars into three main groups, none of which correspond to major geographic regions, suggesting that the genetics and geographical origin of modern Chinese jujube cultivars might not be linked. The current work firstly reports the large-scale development of Chinese jujube SSR markers. The development of these markers and their polymorphic information represent a significant improvement in the available Chinese jujube genomic resources and will facilitate both genetic and breeding applications, further accelerating the development of new cultivars.

Improved short-sequence-repeat genotyping of Mycobacterium avium subsp. paratuberculosis by using matrix-assisted laser desorption ionization-time of flight mass spectrometry

Accurate sequence analysis of mononucleotide repeat regions is difficult, complicating the use of short sequence repeats (SSRs) as a tool for bacterial strain discrimination. Although multiple SSR loci in the genome of Mycobacterium avium subsp. paratuberculosis allow genotyping of M. avium subsp. paratuberculosis isolates with high discriminatory power, further characterization of the most discriminatory loci is limited due to inherent difficulties in sequencing mononucleotide repeats. Here, a method was evaluated using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) as an alternative to Sanger sequencing to further differentiate the dominant mycobacterial interspersed repetitive-unit (MIRU)-variable-number tandem-repeat (VNTR) M. avium subsp. paratuberculosis type (n = 37) in Canadian dairy herds by targeting a highly discriminatory mononucleotide SSR locus. First, PCR-amplified DNA was digested with two restriction enzymes to yield a sufficiently small fragment containing the SSR locus. Second, MALDI-TOF MS was performed to identify the mass, and thus repeat length, of the target. Sufficiently intense, discriminating spectra were obtained to determine repeat lengths up to 15, an improvement over the limit of 11 using traditional sequencing techniques. Comparison to synthetic oligonucleotides and Sanger sequencing results confirmed a valid and reproducible assay that increased discrimination of the dominant M. avium subsp. paratuberculosis MIRU-VNTR type. Thus, MALDI-TOF MS was a reliable, fast, and automatable technique to accurately resolve M. avium subsp. paratuberculosis genotypes based on SSRs.