Comparing inferences derived from microsatellite and RADseq datasets: a case study involving threatened bull trout

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Tissue-specific transcriptomes reveal potential mechanisms of microbiome heterogeneity in an ancient fish

The lake sturgeon (Acipenser fulvescens) is an ancient, octoploid fish faced with conservation challenges across its range in North America, but a lack of genomic resources has hindered molecular research in the species. To support such research, we created a transcriptomic database from 13 tissues: brain, esophagus, gill, head kidney, heart, white muscle, liver, glandular stomach, muscular stomach, anterior intestine, pyloric cecum, spiral valve and rectum. The transcriptomes for each tissue were sequenced and assembled individually from a mean of 98.3 million (±38.9 million SD) reads each. In addition, an overall transcriptome was assembled and annotated with all data used for each tissue-specific transcriptome. All assembled transcriptomes and their annotations were made publicly available as a scientific resource. The non-gut transcriptomes provide important resources for many research avenues. However, we focused our analysis on messenger ribonucleic acid (mRNA) observations in the gut because the gut represents a compartmentalized organ system with compartmentalized functions, and seven of the sequenced tissues were from each of these portions. These gut-specific analyses were used to probe evidence of microbiome regulation by studying heterogeneity in microbial genes and genera identified from mRNA annotations. Gene set enrichment analyses were used to reveal the presence of photoperiod and circadian-related transcripts in the pyloric cecum, which may support periodicity in lake sturgeon digestion. Similar analyses were used to identify different types of innate immune regulation across the gut, while analyses of unique transcripts annotated to microbes revealed heterogeneous genera and genes among different gut tissues. The present results provide a scientific resource and information about the mechanisms of compartmentalized function across gut tissues in a phylogenetically ancient vertebrate. Database URL: https://figshare.com/projects/Lake_Sturgeon_Transcriptomes/133143.

Scheme of Effective Epidemiological Investigations in Trichinella Outbreaks on Pig Farms

Trichinellosis is a parasitic, zoonotic disease caused by larvae of the genus Trichinella. Infection occurs via the consumption of raw or undercooked meat containing this parasite. Symptoms of the disease manifest as intestinal disorders, followed by facial swelling, fever, muscle pain and other symptoms, eventually leading to neurological and cardiac complications and even death. In Europe, trichinellosis is most often associated with the consumption of meat from wild boars, pigs and horses. In recent years, wild boars that are hunted illegally and not tested for Trichinella spp. have been the most common cause of trichinellosis in humans; however, there have also been cases where infected pigs have been the source of infection. When trichinellosis is suspected in humans, epidemiological measures are taken to identify the source. Similarly, an epidemiological investigation should be initiated whenever Trichinella spp. has been detected in pigs. However, commonly used actions do not provide sufficient data to determine the source of infection for pigs and to prevent further transmission. Therefore, in this article, we propose a scheme for effective epidemiological investigations into Trichinella outbreaks on pig farms that can help trace the transmission mechanisms of the parasite and that takes into account currently available testing tools. The proposed pathway can be easily adopted for epidemiological investigations in routine veterinary inspection work.

As the Goose Flies: Migration Routes and Timing Influence Patterns of Genetic Diversity in a Circumpolar Migratory Herbivore

Migration schedules and the timing of other annual events (e.g., pair formation and molt) can affect the distribution of genetic diversity as much as where these events occur. The greater white-fronted goose (Anser albifrons) is a circumpolar goose species, exhibiting temporal and spatial variation of events among populations during the annual cycle. Previous range-wide genetic assessments of the nuclear genome based on eight microsatellite loci suggest a single, largely panmictic population despite up to five subspecies currently recognized based on phenotypic differences. We used double digest restriction-site associated DNA (ddRAD-seq) and mitochondrial DNA (mtDNA) sequence data to re-evaluate estimates of spatial genomic structure and to characterize how past and present processes have shaped the patterns of genetic diversity and connectivity across the Arctic and subarctic. We uncovered previously undetected inter-population differentiation with genetic clusters corresponding to sampling locales associated with current management groups. We further observed subtle genetic clustering within each management unit that can be at least partially explained by the timing and directionality of migration events along with other behaviors during the annual cycle. The Tule Goose (A. a. elgasi) and Greenland subspecies (A. a. flavirostris) showed the highest level of divergence among all sampling locales investigated. The recovery of previously undetected broad and fine-scale spatial structure suggests that the strong cultural transmission of migratory behavior restricts gene flow across portions of the species’ range. Our data further highlight the importance of re-evaluating previous assessments conducted based on a small number of highly variable genetic markers in phenotypically diverse species.

RADseq data reveal a lack of admixture in a mouse lemur contact zone contrary to previous microsatellite results

Microsatellites have been a workhorse of evolutionary genetic studies for decades and are still commonly in use for estimating signatures of genetic diversity at the population and species level across a multitude of taxa. Yet, the very high mutation rate of these loci is a double-edged sword, conferring great sensitivity at shallow levels of analysis (e.g. paternity analysis) but yielding considerable uncertainty for deeper evolutionary comparisons. For the present study, we used reduced representation genome-wide data (restriction site-associated DNA sequencing (RADseq)) to test for patterns of interspecific hybridization previously characterized using microsatellite data in a contact zone between two closely related mouse lemur species in Madagascar (Microcebus murinus and Microcebus griseorufus). We revisit this system by examining populations in, near, and far from the contact zone, including many of the same individuals that had previously been identified as hybrids with microsatellite data. Surprisingly, we find no evidence for admixed nuclear ancestry. Instead, re-analyses of microsatellite data and simulations suggest that previously inferred hybrids were false positives and that the program NewHybrids can be particularly sensitive to erroneously inferring hybrid ancestry. Combined with results from coalescent-based analyses and evidence for local syntopic co-occurrence, we conclude that the two mouse lemur species are in fact completely reproductively isolated, thus providing a new understanding of the evolutionary rate whereby reproductive isolation can be achieved in a primate.

Molecular and morphological analyses clarify species delimitation in section Costatae and reveal Betula buggsii sp. nov. (sect. Costatae, Betulaceae) in China

BACKGROUND AND AIMS

Delineating closely related and morphologically similar species is difficult. Here, we integrate morphology, genetics, ploidy and geography to resolve species and subspecies boundaries in four trees of section Costatae (genus Betula): Betula ashburneri, B. costata, B. ermanii and B. utilis, as well as multiple subspecies and polyploid races.

METHODS

We genotyped 371 individuals (20-133 per species) from 51 populations at 15 microsatellite markers, as well as a subset of individuals, using restriction-site associated DNA sequencing and nuclear internal transcribed spacers. We determined the ploidy level of eight individuals using flow cytometry and characterized leaf variation for a subset of 109 individuals by morphometric analysis.

KEY RESULTS

Integration of multiple lines of evidence suggested a series of revisions to the taxonomy of section Costatae. Betula costata and B. ermanii were found to be valid. Molecular and leaf morphology analyses revealed little differentiation between diploid B. albosinensis and some samples of B. utilis ssp. utilis. By contrast, other B. utilis ssp. utilis samples and ssp. albosinensis formed a morphological continuum but differed based on genetics. Specifically, B. utilis ssp. albosinensis was divided into two groups with group I genetically similar to B. utilis ssp. utilis and group II, a distinct cluster, proposed as the new diploid species Betula buggsii sp. nov. Phylogenomic analysis based on 2285 620 single nucleotide polymorphisms identified a well-supported monophyletic clade of B. buggsii. Morphologically, B. buggsii is characterized by elongated lenticels and a distinct pattern of bark peeling and may be geographically restricted to the Qinling-Daba Mountains.

CONCLUSIONS

Our integrated approach identifies six taxa within section Costatae: B. ashburneri, B. buggsii, B. costata, B. utilis ssp. utilis, B. utilis ssp. albosinensis and B. ermanii. Our research demonstrates the value of an integrative approach using morphological, geographical, genetic and ploidy-level data for species delineation.

Revealing the Genetic Structure and Differentiation in Endangered Pinus bungeana by Genome-Wide SNP Markers

Understanding the genetic variation and differentiation of natural populations is essential for their protection, specifically if the species status is endangered as with Pinus bungeana. Here, we used 346,840 high density and strong specificity SNP loci to carry out genetic analyses (i.e., genetic diversity, genetic structure, phylogeny, and geographical differentiation) on 52 P. bungeana individuals from 5 populations (4 natural and one artificial) representing the main regions of the species distribution in China. Genetic diversity assessment indicated a trend of genetic diversity gradual decrease from west to east across the species distribution areas. Population genetic structure, PCA and phylogenetic analyses consistently indicated that populations in the central and eastern regions were clustered together, while those from the western regions were separated. Mantel test values indicated the presence of geographic isolation among populations, an important factor contributing to the observed genetic differentiation. The maximum likelihood tree and potential migration events inferred from TreeMix analysis indicated the presence of historical genetic exchanges between the west of Qinling Mountains and the Lvliang Mountains populations. Based on the generated genetic information, in situ and ex situ conservation strategies for P.bungeana germplasm resources are proposed, these strategies could be valuable for the conservation, protection and genetic improvement of this endangered species.

Genome-wide markers reveal differentiation between and within the cryptic sister species, sunset and vermilion rockfish

The vermilion rockfish complex, which consists of the cryptic sister species vermilion and sunset rockfish, is one of the most valuable recreational fisheries on the U.S. West Coast. These species are currently managed as a single complex, and because of uncertainty surrounding the relative contribution of each species within existing data sources, the stock status of each species is not fully known. A reliable and cost-effective method is needed to disentangle these species that will allow for the development of abundance indices, life history profiles, and catch histories that may potentially support species-specific stock assessments. Using restriction-site associated DNA sequence (RADseq) markers we generated 10,003 polymorphic loci to characterize the vermilion rockfish complex. PCA and Bayesian clustering approaches based on these loci clearly distinguished between sunset and vermilion rockfishes and identified hybrid individuals. These loci included 203 highly differentiated (FST ≥ 0.99) single nucleotide polymorphisms, which we consider candidates in the planned development of a diagnostic assay capable of distinguishing between these cryptic species. In addition to clearly delineating to species, subsets of the interspecific markers allowed for insight into intraspecific differentiation in both species. Population genetic analyses for sunset rockfish identified two weakly divergent genetic groups with similar levels of genetic diversity. Vermilion rockfish, however, were characterized by three distinct genetic groups with much stronger signals of differentiation and significantly different genetic diversities. Collectively, these data will contribute to well-informed, species-specific management strategies to protect this valuable species complex.

Trich-tracker - a practical tool to trace Trichinella spiralis transmission based on rapid, cost-effective sampling of genome-wide genetic variation

Molecular epidemiology using traditional sequencing has been notoriously difficult in inbred parasites due to a lack of genetic variation available for discriminating among parasites. Next generation sequencing techniques offer a solution to this problem by increasing the number of loci that can be sequenced. Here, we introduce Trich-tracker, a tool that makes efficient use of diagnostic variation distributed throughout the genome of Trichinella spiralis to more rapidly, and conclusively, resolve connections and distinctions among focal outbreaks of T. spiralis. In particular, we rapidly characterised genetic variation among a sample of parasites from Polish farms and wildlife, sampling genomic variation using double digest restriction site-associated DNA sequencing (ddRADseq). Approximately 400,000 bases of sequence were generated from each sample and shown to be distributed across the genome with single nucleotide polymorphisms occurring at a frequency of approximately one base in 10,000. Both phylogenetic and Bayesian clustering analyses indicated that ddRADseq genotypes formed distinct clusters for specific outbreaks and were quite distinct from wild boar samples. Two of the investigated outbreaks were more similar to each other than to other outbreak samples, suggesting a link between these outbreaks. Hence, the Trich-tracker procedure identified informative genomic variation which afforded unprecedented epidemiological resolution. Trich-tracker is very flexible tool, quickly and inexpensively mining genomes of even highly inbred populations of T. spiralis to support outbreak investigations. The simplicity of the entire procedure, and time and cost effectiveness of Trich-tracker support its practical application in ongoing Trichinella outbreaks. The discriminating power of this tool is tunable and scalable, allowing application in a variety of epidemiological contexts, and is easily adapted to other parasite systems.

Infection, genetics, and evolution of Trichinella: Historical insights and applications to molecular epidemiology

Genetic variation in pathogen populations provides the means to answer questions in disease ecology and transmission, illuminating interactions between genetic traits, environmental exposures, and disease. Such studies elucidate the phylogeny, evolution, transmission and pathogenesis of viruses, bacteria and parasites. Here, we review how such studies have fostered understanding of the biology and epidemiology of zoonotic nematode parasites in the genus Trichinella spp., which impose considerable economic and health burdens by infecting wildlife, livestock, and people. To use such data to define ongoing chains of local transmission and source traceback, researchers first must understand the extent and distribution of genetic variation resident in regional parasite populations. Thus, genetic variability illuminates a population's past as well as its present. Here we review how such data have helped define population dynamics of Trichinella spp. in wild and domesticated hosts, creating opportunities to harness genetic variation in the quest to prevent, track, and contain future outbreaks.

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.

Population genetic structure of the great star coral, Montastraea cavernosa, across the Cuban archipelago with comparisons between microsatellite and SNP markers

Coral reef habitats surrounding Cuba include relatively healthy, well-developed shallow and mesophotic (30–150 m) scleractinian communities at the cross-currents of the Tropical Western Atlantic (TWA). However, Cuba’s coral communities are not immune to the declines observed throughout the TWA, and there is limited information available regarding genetic connectivity, diversity, and structure among these populations. This represents an immense gap in our understanding of coral ecology and population dynamics at both local and regional scales. To address this gap, we evaluated the population genetic structure of the coral Montastraea cavernosa across eight reef sites surrounding Cuba. Colonies were genotyped using nine microsatellite markers and > 9,000 single nucleotide polymorphism (SNP) markers generated using the 2bRAD approach to assess fine-scale genetic structure across these sites. Both the microsatellite and SNP analyses identified patterns of genetic differentiation among sample populations. While the microsatellite analyses did not identify significant genetic structure across the seven shallow M. cavernosa sampling sites, the SNP analyses revealed significant pairwise population differentiation, suggesting that differentiation is greater between eastern and western sites. This study provides insight into methodological differences between microsatellite and SNP markers including potential trade-offs between marker-specific biases, sample size, sequencing costs, and the ability to resolve subtle patterns of population genetic structure. Furthermore, this study suggests that locations in western Cuba may play important roles in this species’ regional metapopulation dynamics and therefore may merit incorporation into developing international management efforts in addition to the local management the sites receive.

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.

Ready-to-use workflows for the implementation of genetic tools in conservation management

We present a conservation genetics tool kit, which offers two ready-to-use workflows for the routine application of genetic methods in conservation management. The workflows were optimized for work load and costs and are accompanied by an easy-to-read and richly illustrated manual with guidelines regarding sampling design, sampling of genetic material, necessary permits, laboratory methods, statistical analyses and documentation of results in a practice-oriented way. The manual also provides a detailed interpretation help for the implementation of the results in conservation management. One workflow deals with the identification of pond-breeding amphibians based on metabarcoding and environmental DNA (eDNA) from water samples. This workflow also discriminates the morphologically similar water frogs (Pelophylax sp.) and other closely related species (e.g. Triturus cristatus and T. carnifex). The second workflow studies connectivity among populations using microsatellite markers. Its statistical analyses encompass the detection of genetic groups and historical, recent and current dispersal and gene flow. Using the two workflows does not involve academic research institutes; they can be applied by environmental consultancies, laboratories from the private sector, governmental agencies or non-governmental organisations. These and additional conservation genetic workflows will hopefully foster the routine use of genetic methods in conservation management.

Evaluating the effect of reference genome divergence on the analysis of empirical RADseq datasets

The advent of high-throughput sequencing (HTS) has made genomic-level analyses feasible for nonmodel organisms. A critical step of many HTS pipelines involves aligning reads to a reference genome to identify variants. Despite recent initiatives, only a fraction of species has publically available reference genomes. Therefore, a common practice is to align reads to the genome of an organism related to the target species; however, this could affect read alignment and bias genotyping. In this study, I conducted an experiment using empirical RADseq datasets generated for two species of salmonids (Actinopterygii; Teleostei; Salmonidae) to address these questions. There are currently reference genomes for six salmonids of varying phylogenetic distance. I aligned the RADseq data to all six genomes and identified variants with several different genotypers, which were then fed into population genetic analyses. Increasing phylogenetic distance between target species and reference genome reduced the proportion of reads that successfully aligned and mapping quality. Reference genome also influenced the number of SNPs that were generated and depth at those SNPs, although the affect varied by genotyper. Inferences of population structure were mixed: increasing reference genome divergence reduced estimates of differentiation but similar patterns of population relationships were found across scenarios. These findings reveal how the choice of reference genome can influence the output of bioinformatic pipelines. It also emphasizes the need to identify best practices and guidelines for the burgeoning field of biodiversity genomics.

Comparing the Performance of Microsatellites and RADseq in Population Genetic Studies: Analysis of Data for Pike (Esox lucius) and a Synthesis of Previous Studies

Population genetic studies reveal biodiversity patterns and inform about drivers of evolutionary differentiation and adaptation, including gene flow, drift and selection. This can advance our understanding and aid decision making regarding management and conservation efforts. Microsatellites have long been used in population genetic studies. Thanks to the development of newer techniques, sequencing approaches such as restriction site associated DNA sequencing (RADseq) are on their way to replace microsatellites for some applications. However, the performance of these two marker types in population genetics have rarely been systematically compared. We utilized three neutrally and adaptively differentiated populations of anadromous pike (Esox lucius) to assess the relative performance of microsatellites and RADseq with respect to resolution and conclusiveness of estimates of population differentiation and genetic structure. To this end, the same set of individuals (N = 64) were genotyped with both RADseq and microsatellite markers. To assess effects of sample size, the same subset of 10 randomly chosen individuals from each population (N = 30 in total) were also genotyped with both methods. Comparisons of estimated genetic diversity and structure showed that both markers were able to uncover genetic structuring. The full RADseq dataset provided the clearest detection of the finer scaled genetic structuring, and the other three datasets (full and subset microsatellite, and subset RADseq) provided comparable results. A search for outlier loci performed on the full SNP dataset pointed to signs of selection potentially associated with salinity and temperature, exemplifying the utility of RADseq to inform about the importance of different environmental factors. To evaluate whether performance differences between the markers are general or context specific, the results of previous studies that have investigated population structure using both marker types were synthesized. The synthesis revealed that RADseq performed as well as, or better than microsatellites in detecting genetic structuring in the included studies. The differences in the ability to detect population structure, both in the present and the previous studies, are likely explained by the higher number of loci typically utilized in RADseq compared to microsatellite analysis, as increasing the number of markers will (regardless of the marker type) increase power and allow for clearer detection and higher resolution of genetic structure.