Multi-individual microsatellite identification: A multiple genome approach to microsatellite design (MiMi)

Novel Microsatellite Loci, Cross-Species Validation of Multiplex Assays, and By-Catch Mitochondrial Genomes on Ochthebius Beetles from Supratidal Rockpools

Here we focus on designing, for the first time, microsatellite markers for evolutionary and ecological research on aquatic beetles from the genus Ochthebius (Coleoptera, Hydraenidae). Some of these non-model species, with high cryptic diversity, exclusively inhabit supratidal rockpools, extreme and highly dynamic habitats with important anthropogenic threats. We analysed 15 individuals of four species (O. lejolisii, O. subinteger, O. celatus, and O. quadricollis) across 10 localities from the Mediterranean coasts of Spain and Malta. Using next-generation sequencing technology, two libraries were constructed to interpret the species of the two subgenera present consistently (Ochthebius s. str., O. quadricollis; and Cobalius, the rest of the species). Finally, 20 markers (10 for each subgenus) were obtained and successfully tested by cross-validation in the four species under study. As a by-catch, we could retrieve the complete mitochondrial genomes of O. lejolisii, O. quadricollis, and O. subinteger. Interestingly, the mitochondrial genome of O. quadricollis exhibited high genetic variability compared to already published data. The novel SSR panels and mitochondrial genomes for Ochthebius will be valuable in future research on species identification, diversity, genetic structure, and population connectivity in highly dynamic and threatened habitats such as supratidal coastal rockpools.

Genotypic detection of barriers to rat dispersal: Rattus rattus behind a peninsula predator-proof fence

Clear delimitation of management units is essential for effective management of invasive species. Analysis of population genetic structure of target species can improve identification and interpretation of natural and artificial barriers to dispersal. In Aotearoa New Zealand where the introduced ship rat (Rattus rattus) is a major threat to native biodiversity, effective suppression of pest numbers requires removal and limitation of reinvasion from outside the managed population. We contrasted population genetic structure in rat populations over a wide scale without known barriers, with structure over a fine scale with potential barriers to dispersal. MtDNA D-loop sequences and microsatellite genotypes resolved little genetic structure in southern North Island population samples of ship rat 100 km apart. In contrast, samples from major islands differed significantly for both mtDNA and nuclear markers. We also compared ship rats collected within a small peninsula reserve bounded by sea, suburbs and, more recently, a predator fence with rats in the surrounding forest. Here, mtDNA did not differ but genotypes from 14 nuclear loci were sufficient to distinguish the fenced population. This suggests that natural (sea) and artificial barriers (town, fence) are effectively limiting gene flow among ship rat populations over the short distance (~ 500 m) between the peninsula reserve and surrounding forest. The effectiveness of the fence alone is not clear given it is a recent feature and no historical samples exist; resampling population genetic diversity over time will improve understanding. Nonetheless, the current genetic isolation of the fenced rat population suggests that rat eradication is a sensible management option given that reinvasion appears to be limited and could probably be managed with a biosecurity programme.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10530-023-03004-8.

Pipeline for developing polymorphic microsatellites in species without reference genomes

Microsatellites, also known as simple sequence repeats (SSRs), are the preferred type of marker for many genetic applications. In conjunction with the ongoing development of next-generation sequencing, several bioinformatic tools have been developed for identifying SSRs from genomic or transcriptomic sequences. Although these tools are handy for generating polymorphic SSRs, their application almost always depends on an existing reference genome or self-assembly of the reference genome. With this in mind, we propose a pipeline for developing polymorphic SSRs that may be applied to species without reference genomes. Using a species without a reference genome (black Amur bream; Megalobrama terminalis Richardson, 1846) as a model, our pipeline was able to effectively discover polymorphic SSRs. Under different R parameters of a reference-free single nucleotide polymorphisms (SNPs) caller (ebwt2InDel), a total of 258, 208, 102, and 11 polymorphic SSRs were mined. To quantify the accuracy of the polymorphic SSRs detected using our pipeline, we analyzed 25 SSRs with PCR experiments. All primers were successfully amplified, and most SSRs (23 SSRs, 92%) were polymorphic. From the 36 individual black Amur bream, we acquired an average of 3.36 alleles per locus, ranging from one to 11. This demonstrates the effectiveness of our pipeline in identifying polymorphic SSRs and designing primers for SSR genotyping. Ultimately, our pipeline can effectively mine polymorphic SSRs for species without reference genomes, complementing SSR mining approaches based on reference genomes and helping to resolve biological issues that accompany these methods.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03313-0.

Expression Profiles of Microsatellites in Fruit Tissues of Akebia trifoliata and Development of Efficient EST-SSR Markers

Akebia trifoliata, a member of the family Lardizabalaceae, has high exploitation potential for multiple economic purposes, so genetic improvements to meet requirements for commercial demand are needed. However, this progress is largely impeded by a lack of effective selection markers. In this study, we obtained 271.49 Gb of clean transcriptomic data from 12 samples (three tissues at four developmental stages) of A. trifoliata fruit. We identified 175,604, 194,370, and 207,906 SSRs from the de novo assembled 416,363, 463,756, and 491,680 unigene sequences obtained from the flesh, seed, and rind tissues, respectively. The profile and proportion of SSR motifs expressed in each fruit tissue and developmental stage were remarkably similar, but many trinucleotide repeats had differential expression levels among different tissues or at different developmental stages. In addition, we successfully designed 16,869 functional EST-SSR primers according to the annotated unigenes. Finally, 94 and 72 primer pairs out of 100 randomly selected primer pairs produced clear bands and polymorphic bands, respectively. These results were also used to elucidate the expression profiles of different tissues at various stages. Additionally, we provided a set of effective, polymorphic, and reliable EST-SSR markers sufficient for accelerating the discovery of metabolic and pathway-specific functional genes for genetic improvement and increased commercial productivity.

Characterization of Microsatellites in the Akebia trifoliata Genome and Their Transferability and Development of a Whole Set of Effective, Polymorphic, and Physically Mapped Simple Sequence Repeat Markers

Akebia trifoliata is a perennial climbing woody liana plant with a high potential for commercial exploitation and theoretical research. Similarly, microsatellites (simple sequence repeats, SSRs) also have dual roles: as critical markers and as essential elements of the eukaryotic genome. To characterize the profile of SSRs and develop molecular markers, the high-quality assembled genome of A. trifoliata was used. Additionally, to determine the potential transferability of SSR loci, the genomes of Amborella trichopoda, Oryza sativa, Vitis vinifera, Arabidopsis thaliana, Papaver somniferum, and Aquilegia coerulea were also used. We identified 434,293 SSRs with abundant short repeats, such as 290,868 (66.98%) single-nucleotide repeats (SNRs) and 113,299 (26.09%) dinucleotide repeats (DNRs) in the A. trifoliata genome. 398,728 (91.81%) SSRs on 344,283 loci were physically mapped on the chromosomes, and a positive correlation (r = 0.98) was found between the number of SSRs and chromosomal length. Additionally, 342,916 (99.60%) potential SSR markers could be designed from the 344,283 physically mapped loci, while only 36,160 could be viewed as high-polymorphism-potential (HPP) markers, findings that were validated by PCR. Finally, SSR loci exhibited broad potential transferability, particularly DNRs such as the "AT/AT" and "AG/CT" loci, among all angiosperms, a finding that was not related to the genetic divergence distance. Practically, we developed a whole set of effective, polymorphic, and physically anchored molecular markers and found that, evolutionarily, DNRs could be responsible for microsatellite origin and protecting gene function.

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.

Automating microsatellite screening and primer design from multi-individual libraries using Micro-Primers

Analysis of intra- and inter-population diversity has become important for defining the genetic status and distribution patterns of a species and a powerful tool for conservation programs, as high levels of inbreeding could lead into whole population extinction in few generations. Microsatellites (SSR) are commonly used in population studies but discovering highly variable regions across species' genomes requires demanding computation and laboratorial optimization. In this work, we combine next generation sequencing (NGS) with automatic computing to develop a genomic-oriented tool for characterizing SSRs at the population level. Herein, we describe a new Python pipeline, named Micro-Primers, designed to identify, and design PCR primers for amplification of SSR loci from a multi-individual microsatellite library. By combining commonly used programs for data cleaning and microsatellite mining, this pipeline easily generates, from a fastq file produced by high-throughput sequencing, standard information about the selected microsatellite loci, including the number of alleles in the population subset, and the melting temperature and respective PCR product of each primer set. Additionally, potential polymorphic loci can be identified based on the allele ranges observed in the population, to easily guide the selection of optimal markers for the species. Experimental results show that Micro-Primers significantly reduces processing time in comparison to manual analysis while keeping the same quality of the results. The elapsed times at each step can be longer depending on the number of sequences to analyze and, if not assisted, the selection of polymorphic loci from multiple individuals can represent a major bottleneck in population studies.

Development of the first microsatellite markers using high-throughput sequencing for a hexaploid coastal species, Tournefortia argentea L. f. (Boraginaceae)

BACKGROUD

Tournefortia argentea L. f. is a hexaploid shrub or tree species with ecological and evolutionary significances, which forms the fringe of vegetation closest to the sea on tropical coral islands. Previous studies have never addressed on genetic information, and thus genomic resources remain scarce.

METHODS AND RESULTS

We used nine individuals from different islands to identify polymorphic microsatellites of T. argentea by Illumina high-throughput sequencing. Thirty-five polymorphic microsatellite markers were developed. Characteristics of each locus were tested using 48 individuals collected from three populations of T. argentea. A total of 320 alleles were found across the 35 microsatellite loci. The number of alleles per locus ranged from 5 to 15, with an average of 9.1. Observed and expected heterozygosities in each locus per population varied from 0.000 to 1.000 and from 0.000 to 0.893, respectively.

CONCLUSIONS

In this study, we report the development of 35 polymorphic microsatellite markers based on Illumina high-throughput sequencing. These markers will facilitate the investigations of genetic diversity, population structures and evolutionary history of T. argentea.

Multi-individual microsatellite identification: A multiple genome approach to microsatellite design (MiMi)

Bespoke microsatellite marker panels are increasingly affordable and tractable to researchers and conservationists. The rate of microsatellite discovery is very high within a shotgun genomic data set, but extensive laboratory testing of markers is required for confirmation of amplification and polymorphism. By incorporating shotgun next-generation sequencing data sets from multiple individuals of the same species, we have developed a new method for the optimal design of microsatellite markers. This new tool allows us to increase the rate at which suitable candidate markers are selected by 58% in direct comparisons and facilitate an estimated 16% reduction in costs associated with producing a novel microsatellite panel. Our method enables the visualisation of each microsatellite locus in a multiple sequence alignment allowing several important quality checks to be made. Polymorphic loci can be identified and prioritised. Loci containing fragment-length-altering mutations in the flanking regions, which may invalidate assumptions regarding the model of evolution underlying variation at the microsatellite, can be avoided. Priming regions containing point mutations can be detected and avoided, helping to reduce sample-site-marker specificity arising from genetic isolation, and the likelihood of null alleles occurring. We demonstrate the utility of this new approach in two species: an echinoderm and a bird. Our method makes a valuable contribution towards minimising genotyping errors and reducing costs associated with developing a novel marker panel. The Python script to perform our method of multi-individual microsatellite identification (MiMi) is freely available from GitHub (https://github.com/graemefox/mimi).