Genetic sex assignment in wild populations using genotyping‐by‐sequencing data: A statistical threshold approach

Establishing the sex of individuals in wild systems can be challenging and often requires genetic testing. Genotyping‐by‐sequencing (GBS) and other reduced‐representation DNA sequencing (RRS) protocols (e.g., RADseq, ddRAD) have enabled the analysis of genetic data on an unprecedented scale. Here, we present a novel approach for the discovery and statistical validation of sex‐specific loci in GBS data sets. We used GBS to genotype 166 New Zealand fur seals (NZFS, Arctocephalus forsteri) of known sex. We retained monomorphic loci as potential sex‐specific markers in the locus discovery phase. We then used (i) a sex‐specific locus threshold (SSLT) to identify significantly male‐specific loci within our data set; and (ii) a significant sex‐assignment threshold (SSAT) to confidently assign sex in silico the presence or absence of significantly male‐specific loci to individuals in our data set treated as unknowns (98.9% accuracy for females; 95.8% for males, estimated via cross‐validation). Furthermore, we assigned sex to 86 individuals of true unknown sex using our SSAT and assessed the effect of SSLT adjustments on these assignments. From 90 verified sex‐specific loci, we developed a panel of three sex‐specific PCR primers that we used to ascertain sex independently of our GBS data, which we show amplify reliably in at least two other pinniped species. Using monomorphic loci normally discarded from large SNP data sets is an effective way to identify robust sex‐linked markers for nonmodel species. Our novel pipeline can be used to identify and statistically validate monomorphic and polymorphic sex‐specific markers across a range of species and RRS data sets.     

Isolation and characterization of microsatellites in Leuciscus cephalus (Cypriniformes,Cyprinidae) and cross‐species amplification within the family Cyprinidae

Thirteen polymorphic microsatellites were isolated from Leuciscus cephalus, a widespread cyprinid species with great ecological tolerance. Together with the cross‐species amplification of six additional loci originally published for three cyprinid fish species, we optimized a multiplex panel for L. cephalus allowing the genotyping of 19 polymorphic loci. Number of alleles and heterozygosity per locus in a sample of 20 fish individuals ranged from two to 16 and from 0.05 to 0.90, respectively. These primers will be useful in determining the population structure of L. cephalus. In addition, successful cross‐amplification was obtained for several species of Cyprinidae.     

Genotype‐by‐sequencing of the plant‐pathogenic fungi Pyrenophora teres and Sphaerulina musiva utilizing Ion Torrent sequence technology

Genetic and genomics tools to characterize host–pathogen interactions are disproportionately directed to the host because of the focus on resistance. However, understanding the genetics of pathogen virulence is equally important and has been limited by the high cost of de novo genotyping of species with limited marker data. Non‐resource‐prohibitive methods that overcome the limitation of genotyping are now available through genotype‐by‐sequencing (GBS). The use of a two‐enzyme restriction‐associated DNA (RAD)‐GBS method adapted for Ion Torrent sequencing technology provided robust and reproducible high‐density genotyping of several fungal species. A total of 5783 and 2373 unique loci, ‘sequence tags’, containing 16 441 and 9992 single nucleotide polymorphisms (SNPs) were identified and characterized from natural populations of Pyrenophora teres f. maculata and Sphaerulina musiva, respectively. The data generated from the P. teres f. maculata natural population were used in association mapping analysis to map the mating‐type gene to high resolution. To further validate the methodology, a biparental population of P. teres f. teres, previously used to develop a genetic map utilizing simple sequence repeat (SSR) and amplified fragment length polymorphism (AFLP) markers, was re‐analysed using the SNP markers generated from this protocol. A robust genetic map containing 1393 SNPs on 997 sequence tags spread across 15 linkage groups with anchored reference markers was generated from the P. teres f. teres biparental population. The robust high‐density markers generated using this protocol will allow positional cloning in biparental fungal populations, association mapping of natural fungal populations and population genetics studies.     

Paralogs are revealed by proportion of heterozygotes and deviations in read ratios in genotyping‐by‐sequencing data from natural populations

Whole‐genome duplications have occurred in the recent ancestors of many plants, fish, and amphibians, resulting in a pervasiveness of paralogous loci and the potential for both disomic and tetrasomic inheritance in the same genome. Paralogs can be difficult to reliably genotype and are often excluded from genotyping‐by‐sequencing (GBS) analyses; however, removal requires paralogs to be identified which is difficult without a reference genome. We present a method for identifying paralogs in natural populations by combining two properties of duplicated loci: (i) the expected frequency of heterozygotes exceeds that for singleton loci, and (ii) within heterozygotes, observed read ratios for each allele in GBS data will deviate from the 1:1 expected for singleton (diploid) loci. These deviations are often not apparent within individuals, particularly when sequence coverage is low; but, we postulated that summing allele reads for each locus over all heterozygous individuals in a population would provide sufficient power to detect deviations at those loci. We identified paralogous loci in three species: Chinook salmon (Oncorhynchus tshawytscha) which retains regions with ongoing residual tetrasomy on eight chromosome arms following a recent whole‐genome duplication, mountain barberry (Berberis alpina) which has a large proportion of paralogs that arose through an unknown mechanism, and dusky parrotfish (Scarus niger) which has largely rediploidized following an ancient whole‐genome duplication. Importantly, this approach only requires the genotype and allele‐specific read counts for each individual, information which is readily obtained from most GBS analysis pipelines.     

GIbPSs: a toolkit for fast and accurate analyses of genotyping‐by‐sequencing data without a reference genome

Genotyping‐by‐sequencing (GBS) and related methods are increasingly used for studies of non‐model organisms from population genetic to phylogenetic scales. We present GIbPSs, a new genotyping toolkit for the analysis of data from various protocols such as RAD, double‐digest RAD, GBS, and two‐enzyme GBS without a reference genome. GIbPSs can handle paired‐end GBS data and is able to assign reads from both strands of a restriction fragment to the same locus. GIbPSs is most suitable for population genetic and phylogeographic analyses. It avoids genotyping errors due to indel variation by identifying and discarding affected loci. GIbPSs creates a genotype database that offers rich functionality for data filtering and export in numerous formats. We performed comparative analyses of simulated and real GBS data with GIbPSs and another program, pyRAD. This program accounts for indel variation by aligning homologous sequences. GIbPSs performed better than pyRAD in several aspects. It required much less computation time and displayed higher genotyping accuracy. GIbPSs retained smaller numbers of loci overall in analyses of real GBS data. It nevertheless delivered more complete genotype matrices with greater locus overlap between individuals and greater numbers of loci sampled in all individuals.     

Cross‐species comparison of microsatellite loci in the Culex pipiens complex and beyond

In the past, we have developed microsatellite loci from the two most common members of the Culex pipiens complex, Culex quinquefasciatus and Culex pipiens. Here we describe seven additional loci and present an extensive survey of a panel of 20 loci across most of the species and subspecies in the complex as well as in morphologically related species. Because we observed a high degree of polymorphism in the flanking regions, we designed new primers and surveyed multiple populations. We present alternate primers and discuss the cross‐species usefulness of these Culex microsatellite loci in a phylogenetic context.     

Application of a SSR‐GBS marker system on investigation of European Hedgehog species and their hybrid zone dynamics

By applying second‐generation sequencing technologies to microsatellite genotyping, sequence information is produced which can result in high‐resolution population genetics analysis populations and increased replicability between runs and laboratories. In the present study, we establish an approach to study the genetic structure patterns of two European hedgehog species Erinaceaus europaeus and E. roumanicus. These species are usually associated with human settlements and are good models to study anthropogenic impacts on the genetic diversity of wild populations. The short sequence repeats genotyping by sequence (SSR‐GBS) method presented uses amplicon sequences to determine genotypes for which allelic variants can be defined according to both length and single nucleotide polymorphisms (SNPs). To evaluate whether complete sequence information improved genetic structure definition, we compared this information with datasets based solely on length information. We identified a total of 42 markers which were successfully amplified in both species. Overall, genotyping based on complete sequence information resulted in a higher number of alleles, as well as greater genetic diversity and differentiation between species. Additionally, the structure patterns were slightly clearer with a division between both species and some potential hybrids. There was some degree of genetic structure within species, although only in E. roumanicus was this related to geographical distance. The statistically significant results obtained by SSR‐GBS demonstrate that it is superior to electrophoresis‐based methods for SSR genotyping. Moreover, the greater reproducibility and throughput with lower effort which can be obtained with SSR‐GBS and the possibility to include degraded DNA into the analysis, allow for continued relevance of SSR markers during the genomic era.     

Isolation of microsatellite loci and cross‐species amplifications in three gobiid fish of the genus Pomatoschistus

Gobiids of the genus Pomatoschistus are increasingly being investigated as models for adaptation to coastal environments and for mating system studies. Among the dozen currently analysed species, microsatellite primers have been characterized only for Pomatoschistus minutus. This paper describes seven new polymorphic loci isolated from Pomatoschistus marmoratus and Pomatoschistus microps, two species that hybridize. Cross‐species amplification was tested for these new loci, together with seven already published P. minutus loci. Systematic amplification of samples of each of the three species provided a first indication of their polymorphism.     

New microsatellite resources for groupers (Serranidae)

We characterized nine microsatellite loci and identified an additional 60 genomic regions containing microsatellites in the red hind grouper, Epinephelus guttatus. The nine loci were highly polymorphic, and primers designed from red hind genomic DNA produced a strong amplification product in a test panel of 16 other groupers in the genera Epinephelus and Mycteroperca collected from across the world. Most of the amplified regions were homologous to the red hind locus and a well‐defined microsatellite repeat was generally evident. The nine loci, together with the 60 uncharacterized microsatellite‐containing regions, provide a powerful tool for ecological and evolutionary studies in groupers.     

Isolation,characterization and cross‐species amplification of microsatellite DNA loci in the tropical American yam Dioscorea trifida

We developed microsatellite markers in American yam (Dioscorea trifida). A microsatellite sequence‐enriched genomic library was screened, and after sequencing, primers were designed for 20 microsatellites. Among these, eight primer pairs yielded amplification products that were both interpretable and polymorphic in 24 yam cultivars. The number of alleles per locus ranged from two to 13 and the overall expected heterozygosity was around 0.5. Six of the eight Dioscorea trifida microsatellite loci gave amplification products in other Dioscorea species.     

Identification of polymorphic microsatellite markers from RAPD product in turbot (Scophthalmus maximus) and a test of cross‐species amplification

Five polymorphic microsatellite loci have been isolated and characterized from random amplified polymorphic DNA product in turbot, Scophthalmus maximus. Twelve microsatellites were selected for designing microsatellite primers, of which five gave working primer pairs. They had between four and nine alleles. Observed and expected heterozygosities varied from 0.76 to 0.90 and from 0.63 to 0.83, respectively. Five additional fish species assessed for cross‐species amplification revealed between one and three positive amplifications and between zero and two polymorphic loci per species.     

Cross-species amplification of 44 microsatellite loci developed for Chaetodon trifascialis, C. lunulatus and C. vagabundus in 22 related butterflyfish species

Forty‐four microsatellite primers developed for three species of butterflyfish were cross‐tested against 22 related confamilial species. Amplification success and cross‐species transferability of these markers were moderately high. Between 24 and 37 loci were amplified successfully in each species, with a mean success rate per species of 71.7% (± 1.8 SE). Rates of amplification success were comparable among primers designed for the three source species, ranging from a mean success rate per species of 16.9 loci (± 0.8 SE) for Chaetodon trifascialis source loci to 13.7 loci (± 1.5 SE) for C. vagabundus source loci. Polymorphism rates were high (76.1%± 3.1 SE of all successfully amplified loci), and 10 loci were polymorphic in all successfully amplified species (Tri14, B11, C5, D3, D113, D6, D117, D120, D111, D118). The number of alleles per polymorphic locus ranged from 2 to 8, and the average number of alleles across all polymorphic loci and all species was 3.6 (± 0.07 SE). Polymorphism rates were higher overall in primers designed for C. vagabundus (89.9%± 3.9 SE). Overall cross‐testing success was lowest for Heniochus chrysostomus, the most phylogenetically divergent species. The significant cross‐testing reported here provides a valuable resource that will enable population genetics studies to be undertaken on a range of butterflyfishes without the need for expensive and time‐consuming de novo microsatellite development.     

Comparison of random and SSR-enriched shotgun pyrosequencing for microsatellite discovery and single multiplex PCR optimization in Acacia harpophylla F. Muell. Ex Benth

Streamlining the development and genotyping of microsatellites in species for which no genetic information is available represents an important technical challenge to overcome in order to enable mainstream application of state-of-the-art population genetic analysis techniques in nonmodel organisms. Using the example of Acacia harpophylla, an acacia tree endemic of north-eastern Australia, we show that high-throughput shotgun pyrosequencing technology, so-called second-generation sequencing, reduces time and cost of microsatellite marker discovery in nonmodel organisms and of their large-scale typing in natural populations. We found that 0.5% of short sequence reads generated on 454 Genome Sequencer FLX Titanium from random genome sampling and 2.2% of reads generated with prior microsatellite enrichment yielded microsatellite markers with designed polymerase chain reaction (PCR) primers, suggesting that enrichment increases efficiency of pyrosequencing when microsatellite discovery is the primary goal. Using stringent selection criteria to facilitate downstream PCR multiplex design, we identified 1435 microsatellite loci with designed primers from a total of 200,908 short sequence reads. From a subset of 96 loci tested for amplification, 38 were validated for population genetics applications, leading to the optimization of a cost-effective multiplex PCR protocol for the simultaneous typing of nine microsatellites in natural populations of A. harpophylla.     

SNP discovery in wild and domesticated populations of blue catfish,Ictalurus furcatus,using genotyping‐by‐sequencing and subsequent SNP validation

Blue catfish, Ictalurus furcatus, are valued in the United States as a trophy fishery for their capacity to reach large sizes, sometimes exceeding 45 kg. Additionally, blue catfish × channel catfish (I. punctatus) hybrid food fish production has recently increased the demand for blue catfish broodstock. However, there has been little study of the genetic impacts and interaction of farmed, introduced and stocked populations of blue catfish. We utilized genotyping‐by‐sequencing (GBS) to capture and genotype SNP markers on 190 individuals from five wild and domesticated populations (Mississippi River, Missouri, D&B, Rio Grande and Texas). Stringent filtering of SNP‐calling parameters resulted in 4275 SNP loci represented across all five populations. Population genetics and structure analyses revealed potential shared ancestry and admixture between populations. We utilized the Sequenom MassARRAY to validate two multiplex panels of SNPs selected from the GBS data. Selection criteria included SNPs shared between populations, SNPs specific to populations, number of reads per individual and number of individuals genotyped by GBS. Putative SNPs were validated in the discovery population and in two additional populations not used in the GBS analysis. A total of 64 SNPs were genotyped successfully in 191 individuals from nine populations. Our results should guide the development of highly informative, flexible genotyping multiplexes for blue catfish from the larger GBS SNP set as well as provide an example of a rapid, low‐cost approach to generate and genotype informative marker loci in aquatic species with minimal previous genetic information.     

Optimization of the genotyping‐by‐sequencing strategy for population genomic analysis in conifers

Flexibility and low cost make genotyping‐by‐sequencing (GBS) an ideal tool for population genomic studies of nonmodel species. However, to utilize the potential of the method fully, many parameters affecting library quality and single nucleotide polymorphism (SNP) discovery require optimization, especially for conifer genomes with a high repetitive DNA content. In this study, we explored strategies for effective GBS analysis in pine species. We constructed GBS libraries using HpaII, PstI and EcoRI‐MseI digestions with different multiplexing levels and examined the effect of restriction enzymes on library complexity and the impact of sequencing depth and size selection of restriction fragments on sequence coverage bias. We tested and compared UNEAK, Stacks and GATK pipelines for the GBS data, and then developed a reference‐free SNP calling strategy for haploid pine genomes. Our GBS procedure proved to be effective in SNP discovery, producing 7000–11 000 and 14 751 SNPs within and among three pine species, respectively, from a PstI library. This investigation provides guidance for the design and analysis of GBS experiments, particularly for organisms for which genomic information is lacking.     

Different filtering strategies of genotyping‐by‐sequencing data provide complementary resolutions of species boundaries and relationships in a clade of sexually deceptive orchids

Ongoing hybridization and retained ancestral polymorphism in rapidly radiating lineages could mask recent cladogenetic events. This presents a challenge for the application of molecular phylogenetic methods to resolve differences between closely related taxa. We reanalyzed published genotyping‐by‐sequencing (GBS) data to infer the phylogeny of four species within the Ophrys sphegodes complex, a recently radiated clade of orchids. We used different data filtering approaches to detect different signals contained in the dataset generated by GBS and estimated their effects on maximum likelihood trees, global F_ST and bootstrap support values. We obtained a maximum likelihood tree with high bootstrap support, separating the species by using a large dataset based on loci shared by at least 30% of accessions. Bootstrap and F_ST values progressively decreased when filtering for loci shared by a higher number of accessions. However, when filtering more stringently to retain homozygous and organellar loci, we identified two main clades. These clades group individuals independently from their a priori species assignment, but were associated with two organellar haplotype clusters. We infer that a less stringent filtering preferentially selects for rapidly evolving lineage‐specific loci, which might better delimit lineages. In contrast, when using homozygous/organellar DNA loci the signature of a putative hybridization event in the lineage prevails over the most recent phylogenetic signal. These results show that using differing filtering strategies on GBS data could dissect the organellar and nuclear DNA phylogenetic signal and yield novel insights into relationships between closely related species.     

Primers from the orders Osteoglossiform and Siluriform detect polymorphic microsatellite loci in sun-catfish, Horabagrus brachysoma (Teleostei: Bagridae)

Horabagrus brachysoma (sun‐catfish, Bagridae, Siluriformes) is a valuable ornamental and food fish. The stock structure of H. brachysoma, necessary to conserve its declining natural populations, is not known. Twenty‐five primers developed for four fish species belonging to the orders Siluriform (3) and Osteoglossiform (1) were tested and eight primers amplified microsatellite loci in H. brachysoma. The results demonstrate that cross‐priming between fish species belonging to different families and even to different orders can yield microsatellite loci. Five of eight primers each amplified two loci. However, the loci that had repeat motifs after sequencing were considered only for genotyping. Finally, eight loci were polymorphic with hree to seven alleles. Individual fish genotype data (n = 42; 21 each in two rivers) at each locus was analysed. Significant genetic heterogeneity was detected at six loci. The identified loci exhibited potential for use in population genetics application in H. brachysoma.