Whole-genome sequencing approaches for conservation biology: Advantages, limitations and practical recommendations

De Novo Genome Assemblies for Three North American Bumble Bee Species: Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii

Bumble bees are ecologically and economically important insect pollinators. Three abundant and widespread species in western North America, Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii, have been the focus of substantial research relating to diverse aspects of bumble bee ecology and evolutionary biology. We present de novo genome assemblies for each of the three species using hybrid assembly of Illumina and Oxford Nanopore Technologies sequences. All three assemblies are of high quality with large N50s (> 2.2 Mb), BUSCO scores indicating > 98% complete genes, and annotations producing 13,325 - 13,687 genes, comparing favorably with other bee genomes. Analysis of synteny against the most complete bumble bee genome, Bombus terrestris, reveals a high degree of collinearity. These genomes should provide a valuable resource for addressing questions relating to functional genomics and evolutionary biology in these species.

Low-coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp

Population divergence and gene flow are key processes in evolution and ecology. Model-based analysis of genome-wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome-level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.

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.

Genome Sequencing of the Endangered Kingdonia uniflora (Circaeasteraceae, Ranunculales) Reveals Potential Mechanisms of Evolutionary Specialization

Kingdonia uniflora, an alpine herb, has an extremely narrow distribution and represents a model for studying evolutionary mechanisms of species that have adapted to undisturbed environments for evolutionarily long periods of time. We assembled a 1,004.7-Mb draft genome (encoding 43,301 genes) of K. uniflora and found significant overrepresentation in gene families associated with DNA repair, underrepresentation in gene families associated with stress response, and loss of most plastid ndh genes. During the evolutionary process, the overrepresentation of gene families involved in DNA repair could help asexual K. uniflora reduce the accumulation of deleterious mutations, while reducing genetic diversity, which is important in responding to environment fluctuations. The underrepresentation of gene families related to stress response and functional loss of ndh genes could be due to lack or loss of ability to respond to environmental changes caused by long-term adaptation to a relatively stable ecological environment.

Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by+ ) from XX genetic females (amhr2by- ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.

Contemporary Demographic Reconstruction Methods Are Robust to Genome Assembly Quality: A Case Study in Tasmanian Devils

Reconstructing species’ demographic histories is a central focus of molecular ecology and evolution. Recently, an expanding suite of methods leveraging either the sequentially Markovian coalescent (SMC) or the site-frequency spectrum has been developed to reconstruct population size histories from genomic sequence data. However, few studies have investigated the robustness of these methods to genome assemblies of varying quality. In this study, we first present an improved genome assembly for the Tasmanian devil using the Chicago library method. Compared with the original reference genome, our new assembly reduces the number of scaffolds (from 35,975 to 10,010) and increases the scaffold N90 (from 0.101 to 2.164 Mb). Second, we assess the performance of four contemporary genomic methods for inferring population size history (PSMC, MSMC, SMC++, Stairway Plot), using the two devil genome assemblies as well as simulated, artificially fragmented genomes that approximate the hypothesized demographic history of Tasmanian devils. We demonstrate that each method is robust to assembly quality, producing similar estimates of Ne when simulated genomes were fragmented into up to 5,000 scaffolds. Overall, methods reliant on the SMC are most reliable between ∼300 generations before present (gbp) and 100 kgbp, whereas methods exclusively reliant on the site-frequency spectrum are most reliable between the present and 30 gbp. Our results suggest that when used in concert, genomic methods for reconstructing species’ effective population size histories 1) can be applied to nonmodel organisms without highly contiguous reference genomes, and 2) are capable of detecting independently documented effects of historical geological events.

Mycobiome diversity: high-throughput sequencing and identification of fungi

Fungi are major ecological players in both terrestrial and aquatic environments by cycling organic matter and channelling nutrients across trophic levels. High-throughput sequencing (HTS) studies of fungal communities are redrawing the map of the fungal kingdom by hinting at its enormous - and largely uncharted - taxonomic and functional diversity. However, HTS approaches come with a range of pitfalls and potential biases, cautioning against unwary application and interpretation of HTS technologies and results. In this Review, we provide an overview and practical recommendations for aspects of HTS studies ranging from sampling and laboratory practices to data processing and analysis. We also discuss upcoming trends and techniques in the field and summarize recent and noteworthy results from HTS studies targeting fungal communities and guilds. Our Review highlights the need for reproducibility and public data availability in the study of fungal communities. If the associated challenges and conceptual barriers are overcome, HTS offers immense possibilities in mycology and elsewhere.

A draft genome of a field-collected Steinernema feltiae strain NW

Advances in sequencing technologies have accelerated our understanding of the complex genetic network of organisms and genomic divergences that are linked to evolutionary processes. While many model organisms and laboratory strains have been sequenced, wild populations are underrepresented in the growing list of sequenced genomes. Here, we present a de novo assembly of Steinernema feltiae, strain NW, collected from a working agricultural field in south central Washington, USA. Leveraging Pacific Biosciences (PacBio) long reads, we sequenced strain NW to a high depth (99×). The resulting de novo assembly is significantly larger than the previous assembly generated from the laboratory strain SN, with a noticeable improvement in continuity and completeness. Comparative analysis of two assemblies revealed numerous single nucleotide polymorphisms (SNPs), breakpoints, and indels present between the two genomes. This alternative genome resource and annotation could benefit the research community to examine the genetic foundation of evolutionary processes as well as genomic variation among conspecific populations.

Applications of Next-Generation Sequencing Technologies and Computational Tools in Molecular Evolution and Aquatic Animals Conservation Studies: A Short Review

Aquatic ecosystems that form major biodiversity hotspots are critically threatened due to environmental and anthropogenic stressors. We believe that, in this genomic era, computational methods can be applied to promote aquatic biodiversity conservation by addressing questions related to the evolutionary history of aquatic organisms at the molecular level. However, huge amounts of genomics data generated can only be discerned through the use of bioinformatics. Here, we examine the applications of next-generation sequencing technologies and bioinformatics tools to study the molecular evolution of aquatic animals and discuss the current challenges and future perspectives of using bioinformatics toward aquatic animal conservation efforts.

Demographic Histories and Genome-Wide Patterns of Divergence in Incipient Species of Shorebirds

Understanding how incipient species are maintained with gene flow is a fundamental question in evolutionary biology. Whole genome sequencing of multiple individuals holds great potential to illustrate patterns of genomic differentiation as well as the associated evolutionary histories. Kentish (Charadrius alexandrinus) and the white-faced (C. dealbatus) plovers, which differ in their phenotype, ecology and behavior, are two incipient species and parapatrically distributed in East Asia. Previous studies show evidence of genetic diversification with gene flow between the two plovers. Under this scenario, it is of great importance to explore the patterns of divergence at the genomic level and to determine whether specific regions are involved in reproductive isolation and local adaptation. Here we present the first population genomic analysis of the two incipient species based on the de novo Kentish plover reference genome and resequenced populations. We show that the two plover lineages are distinct in both nuclear and mitochondrial genomes. Using model-based coalescence analysis, we found that population sizes of Kentish plover increased whereas white-faced plovers declined during the Last Glaciation Period. Moreover, the two plovers diverged allopatrically, with gene flow occurring after secondary contact. This has resulted in low levels of genome-wide differentiation, although we found evidence of a few highly differentiated genomic regions in both the autosomes and the Z-chromosome. This study illustrates that incipient shorebird species with gene flow after secondary contact can exhibit discrete divergence at specific genomic regions and provides basis to further exploration on the genetic basis of relevant phenotypic traits.

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.

The Value of Reference Genomes in the Conservation of Threatened Species

Conservation initiatives are now more crucial than ever-over a million plant and animal species are at risk of extinction over the coming decades. The genetic management of threatened species held in insurance programs is recommended; however, few are taking advantage of the full range of genomic technologies available today. Less than 1% of the 13505 species currently listed as threated by the International Union for Conservation of Nature (IUCN) have a published genome. While there has been much discussion in the literature about the importance of genomics for conservation, there are limited examples of how having a reference genome has changed conservation management practice. The Tasmanian devil (Sarcophilus harrisii), is an endangered Australian marsupial, threatened by an infectious clonal cancer devil facial tumor disease (DFTD). Populations have declined by 80% since the disease was first recorded in 1996. A reference genome for this species was published in 2012 and has been crucial for understanding DFTD and the management of the species in the wild. Here we use the Tasmanian devil as an example of how a reference genome has influenced management actions in the conservation of a species.

Breeding Centers, Private Ranches, and Genomics for Creating Sustainable Wildlife Populations

Human-induced changes to environments are causing species declines. Beyond preserving habitat (in situ), insurance (ex situ) populations are essential to prevent species extinctions. The Conservation Centers for Species Survival (C2S2) is leveraging space of breeding centers and private ranches to produce “source populations”—genetically diverse reservoirs that also support research and reintroductions. The initial focus is on four African antelopes. C2S2 has developed a program, the Source Population Alliance, that emphasizes animals living in spacious, naturalistic conditions in greater numbers than can be accommodated by urban zoos. Simulation modeling demonstrates how herds can rapidly increase population abundance and retain genetic diversity. Advances in genomics and resulting DNA data allow monitoring of genetic diversity and parentage as well as refined decision-making. This approach, neither pure in situ nor ex situ, but rather “sorta situ”, is an innovative way of linking public and private sector resources to ensure that endangered species survive.

Exploring a Pool-seq-only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool-seq data to generate a de novo genome assembly for mining exons, upon which Pool-seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual-based single nucleotide polymorphisms [SNPs]) and from mapping the Pool-seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F ST) between the two introduced populations exceeds that of the naturally sympatric populations (F ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( π ¯  ≈ 0.002 and π ¯  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high-quality reference assembly from a divergent species. We conclude that the Pool-seq-only approach can be suitable for detecting and quantifying genome-wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.

Optimal Management of Genetic Diversity in Subdivided Populations

One of the main objectives of conservation programs is the maintenance of genetic diversity because this provides the adaptive potential of populations to face new environmental challenges. Genetic diversity is generally assessed by means of neutral molecular markers, and it is usually quantified by the expected heterozygosity under Hardy-Weinberg equilibrium and the number of alleles per locus or allelic diversity. These two measures of genetic diversity are complementary because whereas the former is directly related to genetic variance for quantitative traits and, therefore, to the short-term response to selection and adaptation, the latter is more sensitive to population bottlenecks and relates more strongly to the long-term capacity of adaptation of populations. In the context of structured populations undergoing conservation programs, it is important to decide the optimum management strategy to preserve as much genetic diversity as possible while avoiding inbreeding. Here we examine, through computer simulations, the consequences of choosing a conservation strategy based on maximizing either heterozygosity or allelic diversity of single-nucleotide polymorphism haplotypes in a subdivided population. Our results suggest that maximization of allelic diversity can be more efficient in maintaining the genetic diversity of subdivided populations than maximization of expected heterozygosity because the former maintains a larger number of alleles while making a better control of inbreeding. Thus, maximization of allelic diversity should be a recommended strategy in conservation programs for structured populations.

Entomology in the 21st Century: Tackling Insect Invasions, Promoting Advancements in Technology, and Using Effective Science Communication-2018 Student Debates

The 2018 student debates of the Entomological Society of America were held at the Joint Annual Meeting for the Entomological Societies of America, Canada, and British Columbia in Vancouver, BC. Three unbiased introductory speakers and six debate teams discussed and debated topics under the theme 'Entomology in the 21st Century: Tackling Insect Invasions, Promoting Advancements in Technology, and Using Effective Science Communication'. This year's debate topics included: 1) What is the most harmful invasive insect species in the world? 2) How can scientists diffuse the stigma or scare factor surrounding issues that become controversial such as genetically modified organisms, agricultural biotechnological developments, or pesticide chemicals? 3) What new/emerging technologies have the potential to revolutionize entomology (other than Clustered Regularly Interspaced Short Palindromic Repeats)? Introductory speakers and debate teams spent approximately 9 mo preparing their statements and arguments and had the opportunity to share this at the Joint Annual Meeting with an engaged audience.

metapop2: Re-implementation of software for the analysis and management of subdivided populations using gene and allelic diversity

Management programmes often have to make decisions based on the analysis of the genetic properties and diversity of populations. Expected heterozygosity (or gene diversity) and population structure parameters are often used to make recommendations for conservation, such as avoidance of inbreeding or migration across subpopulations. Allelic diversity, however, can also provide complementary and useful information for conservation programmes, as it is highly sensitive to population bottlenecks, and is more related to long-term selection response than heterozygosity. Here we present a completely revised and updated re-implementation of the software metapop for the analysis of diversity in subdivided populations, as well as a tool for the management and dynamic estimation of optimal contributions in conservation programmes. This new update includes computation of allelic diversity for population analysis and management, as well as a simulation mode to forecast the consequences of taking different management strategies over time. Furthermore, the new implementation in C++ includes code optimization and improved memory usage, allowing for fast analysis of large data sets including single nucleotide polymorphism markers, as well as enhanced cross-software and cross-platform compatibility.

Identification of Arabis alpina genomic regions associated with climatic variables along an elevation gradient through whole genome scan

We performed a pooled whole-genome sequencing on samples of the alpine plant Arabis alpina, harvested in ten populations along an elevation gradient in the French Alps. A large dataset of genetic variations was produced as single nucleotide polymorphisms (SNPs). A combined genome scan approach enabled detecting genomic regions associated with a synthetic environmental variable characterizing the climate at each sampling location. Positive loci detected by two methods were retained and belong to 19 regions in the Arabis alpina genome. The most significant region harbors an ortholog of the AtNAC062 gene, encoding a membrane-bound transcription factor described as linking the cold response and pathogen resistance that may confer protection to plants under extended snow coverage at high elevations. Other genes involved in the stress response or in flowering regulation were also detected. Altogether, our results indicated that Arabis alpina represent a suitable model for studying genomic adaptation in alpine perennial plants.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

The future is now: Amplicon sequencing and sequence capture usher in the conservation genomics era

The genomics revolution has initiated a new era of population genetics where genome-wide data are frequently used to understand complex patterns of population structure and selection. However, the application of genomic tools to inform management and conservation has been somewhat rare outside a few well studied species. Fortunately, two recently developed approaches, amplicon sequencing and sequence capture, have the potential to significantly advance the field of conservation genomics. Here, amplicon sequencing refers to highly multiplexed PCR followed by high-throughput sequencing (e.g., GTseq), and sequence capture refers to using capture probes to isolate loci from reduced-representation libraries (e.g., Rapture). Both approaches allow sequencing of thousands of individuals at relatively low costs, do not require any specialized equipment for library preparation, and generate data that can be analyzed without sophisticated computational infrastructure. Here, we discuss the advantages and disadvantages of each method and provide a decision framework for geneticists who are looking to integrate these methods into their research programme. While it will always be important to consider the specifics of the biological question and system, we believe that amplicon sequencing is best suited for projects aiming to genotype <500 loci on many individuals (>1,500) or for species where continued monitoring is anticipated (e.g., long-term pedigrees). Sequence capture, on the other hand, is best applied to projects including fewer individuals or where >500 loci are required. Both of these techniques should smooth the transition from traditional genetic techniques to genomics, helping to usher in the conservation genomics era.

PCR Amplification Strategies Towards Full-length HIV-1 Genome Sequencing

The advent of next-generation sequencing has enabled greater resolution of viral diversity and improved feasibility of full viral genome sequencing allowing routine HIV-1 full genome sequencing in both research and diagnostic settings. Regardless of the sequencing platform selected, successful PCR amplification of the HIV-1 genome is essential for sequencing template preparation. As such, full HIV-1 genome amplification is a crucial step in dictating the successful and reliable sequencing downstream. Here we reviewed existing PCR protocols leading to HIV-1 full genome sequencing. In addition to the discussion on basic considerations on relevant PCR design, the advantages as well as the pitfalls of the published protocols were reviewed.

Seascape genomics of eastern oyster (Crassostrea virginica) along the Atlantic coast of Canada

Interactions between environmental factors and complex life-history characteristics of marine organisms produce the genetic diversity and structure observed within species. Our main goal was to test for genetic differentiation among eastern oyster populations from the coastal region of Canadian Maritimes against expected genetic homogeneity caused by historical events, taking into account spatial and environmental (temperature, salinity, turbidity) variation. This was achieved by genotyping 486 individuals originating from 13 locations using RADSeq. A total of 11,321 filtered SNPs were used in a combination of population genomics and environmental association analyses. We revealed significant neutral genetic differentiation (mean F ST = 0.009) between sampling locations, and the occurrence of six major genetic clusters within the studied system. Redundancy analyses (RDAs) revealed that spatial and environmental variables explained 3.1% and 4.9% of the neutral genetic variation and 38.6% and 12.2% of the putatively adaptive genetic variation, respectively. These results indicate that these environmental factors play a role in the distribution of both neutral and putatively adaptive genetic diversity in the system. Moreover, polygenic selection was suggested by genotype-environment association analysis and significant correlations between additive polygenic scores and temperature and salinity. We discuss our results in the context of their conservation and management implications for the eastern oyster.

Beyond Biodiversity: Can Environmental DNA (eDNA) Cut It as a Population Genetics Tool?

Population genetic data underpin many studies of behavioral, ecological, and evolutionary processes in wild populations and contribute to effective conservation management. However, collecting genetic samples can be challenging when working with endangered, invasive, or cryptic species. Environmental DNA (eDNA) offers a way to sample genetic material non-invasively without requiring visual observation. While eDNA has been trialed extensively as a biodiversity and biosecurity monitoring tool with a strong taxonomic focus, it has yet to be fully explored as a means for obtaining population genetic information. Here, we review current research that employs eDNA approaches for the study of populations. We outline challenges facing eDNA-based population genetic methodologies, and suggest avenues of research for future developments. We advocate that with further optimizations, this emergent field holds great potential as part of the population genetics toolkit.

MobiSeq: De novo SNP discovery in model and non-model species through sequencing the flanking region of transposable elements

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf (Canis lupus), red deer complex (Cervus sp.) and brown rat (Rattus norvegicus). MobiSeq enables the sequencing of hundreds of thousands loci across the genome and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms-enabling the exploration of diverse evolutionary and conservation questions.

Comparing genomic signatures of domestication in two Atlantic salmon (Salmo salar L.) populations with different geographical origins

Selective breeding and genetic improvement have left detectable signatures on the genomes of domestic species. The elucidation of such signatures is fundamental for detecting genomic regions of biological relevance to domestication and improving management practices. In aquaculture, domestication was carried out independently in different locations worldwide, which provides opportunities to study the parallel effects of domestication on the genome of individuals that have been selected for similar traits. In this study, we aimed to detect potential genomic signatures of domestication in two independent pairs of wild/domesticated Atlantic salmon populations of Canadian and Scottish origins, respectively. Putative genomic regions under divergent selection were investigated using a 200K SNP array by combining three different statistical methods based either on allele frequencies (LFMM, Bayescan) or haplotype differentiation (Rsb). We identified 337 and 270 SNPs potentially under divergent selection in wild and hatchery populations of Canadian and Scottish origins, respectively. We observed little overlap between results obtained from different statistical methods, highlighting the need to test complementary approaches for detecting a broad range of genomic footprints of selection. The vast majority of the outliers detected were population-specific but we found four candidate genes that were shared between the populations. We propose that these candidate genes may play a role in the parallel process of domestication. Overall, our results suggest that genetic drift may have override the effect of artificial selection and/or point toward a different genetic basis underlying the expression of similar traits in different domesticated strains. Finally, it is likely that domestication may predominantly target polygenic traits (e.g., growth) such that its genomic impact might be more difficult to detect with methods assuming selective sweeps.

High-throughput identification and diagnostics of pathogens and pests: Overview and practical recommendations

High-throughput identification technologies provide efficient tools for understanding the ecology and functioning of microorganisms. Yet, these methods have been only rarely used for monitoring and testing ecological hypotheses in plant pathogens and pests in spite of their immense importance in agriculture, forestry and plant community dynamics. The main objectives of this manuscript are the following: (a) to provide a comprehensive overview about the state-of-the-art high-throughput quantification and molecular identification methods used to address population dynamics, community ecology and host associations of microorganisms, with a specific focus on antagonists such as pathogens, viruses and pests; (b) to compile available information and provide recommendations about specific protocols and workable primers for bacteria, fungi, oomycetes and insect pests; and (c) to provide examples of novel methods used in other microbiological disciplines that are of great potential use for testing specific biological hypotheses related to pathology. Finally, we evaluate the overall perspectives of the state-of-the-art and still evolving methods for diagnostics and population- and community-level ecological research of pathogens and pests.

Does the Genomic Landscape of Species Divergence in Phaseolus Beans Coerce Parallel Signatures of Adaptation and Domestication?

Exploring the genomic architecture of species and populations divergence aids understanding how lineages evolve and adapt, and ultimately can show the repeatability of evolutionary processes. Yet, the genomic signatures associated with divergence are still relatively unexplored, leading to a knowledge gap on whether species divergence ultimately differs in its genetic architecture from divergence at other spatial scales (i.e., populations, ecotypes). Our goal in this research was to determine whether genomic islands of speciation are more prone to harbor within-species differentiation due to genomic features, suppressed recombination, smaller effective population size or increased drift, across repeated hierarchically nested levels of divergence. We used two species of Phaseolus beans with strong genepool and population sub-structure produced by multiple independent domestications each especially in Andean and Mesoamerican / Middle American geographies. We genotyped 22,531 GBS-derived SNP markers in 209 individuals of wild and cultivated Phaseolus vulgaris and Phaseolus lunatus. We identified six regions for species-associated divergence. Out of these divergence peaks, 21% were recovered in the four within-species between-genepool comparisons and in the five within-genepool wild-cultivated comparisons (some of the latter did retrieve genuine signatures of the well described multiple domestication syndromes). However, genomic regions with overall high relative differentiation (measured by FST) coincided with regions of low SNP density and regions of elevated delta divergence between-genepools (ΔDiv), independent of the scale of divergence. The divergence in chromosome Pv10 further coincided with a between-species pericentric inversion. These convergences suggest that shared variants are being recurrently fixed at replicated regions of the genome, and in a similar manner across different hierarchically nested levels of divergence, likely as result of genomic features that make certain regions more prone to accumulate islands of speciation and within-species divergence. In summary, neighboring signatures of speciation, adaptation and domestication in Phaseolus beans are influenced by ubiquitous genomic constrains, which may continue to fortuitously shape genomic differentiation at various others scales of divergence.

Along-shelf connectivity and circumpolar gene flow in Antarctic silverfish (Pleuragramma antarctica)

The Antarctic silverfish (Pleuragramma antarctica) is a critically important forage species with a circumpolar distribution and is unique among other notothenioid species for its wholly pelagic life cycle. Previous studies have provided mixed evidence of population structure over regional and circumpolar scales. The aim of the present study was to test the recent population hypothesis for Antarctic silverfish, which emphasizes the interplay between life history and hydrography in shaping connectivity. A total of 1067 individuals were collected over 25 years from different locations on a circumpolar scale. Samples were genotyped at fifteen microsatellites to assess population differentiation and genetic structuring using clustering methods, F-statistics, and hierarchical analysis of variance. A lack of differentiation was found between locations connected by the Antarctic Slope Front Current (ASF), indicative of high levels of gene flow. However, gene flow was significantly reduced at the South Orkney Islands and the western Antarctic Peninsula where the ASF is absent. This pattern of gene flow emphasized the relevance of large-scale circulation as a mechanism for circumpolar connectivity. Chaotic genetic patchiness characterized population structure over time, with varying patterns of differentiation observed between years, accompanied by heterogeneous standard length distributions. The present study supports a more nuanced version of the genetic panmixia hypothesis that reflects physical-biological interactions over the life history.

Phylogenomic Analysis of the Explosive Adaptive Radiation of the Espeletia Complex (Asteraceae) in the Tropical Andes

The subtribe Espeletiinae (Asteraceae), endemic to the high-elevations in the Northern Andes, exhibits an exceptional diversity of species, growth-forms, and reproductive strategies. This complex of 140 species includes large trees, dichotomous trees, shrubs and the extraordinary giant caulescent rosettes, considered as a classic example of adaptation in tropical high-elevation ecosystems. The subtribe has also long been recognized as a prominent case of adaptive radiation, but the understanding of its evolution has been hampered by a lack of phylogenetic resolution. Herein, we produce the first fully resolved phylogeny of all morphological groups of Espeletiinae, using whole plastomes and about a million nuclear nucleotides obtained with an original de novo assembly procedure without reference genome, and analyzed with traditional and coalescent-based approaches that consider the possible impact of incomplete lineage sorting and hybridization on phylogenetic inference. We show that the diversification of Espeletiinae started from a rosette ancestor about 2.3 Ma, after the final uplift of the Northern Andes. This was followed by two independent radiations in the Colombian and Venezuelan Andes, with a few trans-cordilleran dispersal events among low-elevation tree lineages but none among high-elevation rosettes. We demonstrate complex scenarios of morphological change in Espeletiinae, usually implying the convergent evolution of growth-forms with frequent loss/gains of various traits. For instance, caulescent rosettes evolved independently in both countries, likely as convergent adaptations to life in tropical high-elevation habitats. Tree growth-forms evolved independently three times from the repeated colonization of lower elevations by high-elevation rosette ancestors. The rate of morphological diversification increased during the early phase of the radiation, after which it decreased steadily towards the present. On the other hand, the rate of species diversification in the best-sampled Venezuelan radiation was on average very high (3.1 spp/My), with significant rate variation among growth-forms (much higher in polycarpic caulescent rosettes). Our results point out a scenario where both adaptive morphological evolution and geographical isolation due to Pleistocene climatic oscillations triggered an exceptionally rapid radiation for a continental plant group.

Range-wide genomic data synthesis reveals transatlantic vicariance and secondary contact in Atlantic cod

Recent advances in genetic and genomic analysis have greatly improved our understanding of spatial population structure in marine species. However, studies addressing phylogeographic patterns at oceanic spatial scales remain rare. In Atlantic cod (Gadus morhua), existing range-wide examinations suggest significant transatlantic divergence, although the fine-scale contemporary distribution of populations and potential for secondary contact are largely unresolved. Here, we explore transatlantic phylogeography in Atlantic cod using a data-synthesis approach, integrating multiple genome-wide single-nucleotide polymorphism (SNP) datasets representative of different regions to create a single range-wide dataset containing 1,494 individuals from 54 locations and genotyped at 796 common loci. Our analysis highlights significant transatlantic divergence and supports the hypothesis of westward post-glacial colonization of Greenland from the East Atlantic. Accordingly, our analysis suggests the presence of transatlantic secondary contact off eastern North America and supports existing perspectives on the phylogeographic history of Atlantic cod with an unprecedented combination of genetic and geographic resolution. Moreover, we demonstrate the utility of integrating distinct SNP databases of high comparability.

Emerging and Enduring Issues in Primate Conservation Genetics

Conservation genetics is a branch of conservation biology that uses molecular data to assist in the conservation and management of imperiled populations, subspecies, and species. In this review, I examine conservation action plans (CAPs)—instrumental documents designed to influence conservation policy—for selected primate species. I use the information contained in CAPs as a means to guide this review. The primary genetics-based topics that are mentioned in CAPs are genetic connectivity, inbreeding, and subspecies/species delimitation. I discuss these topics as well as historical demographic inference and hybridization using examples from wild primate species to illustrate the myriad ways in which genetics can assist in conservation efforts. I also discuss some recent technological advances such as genomic capture techniques and the potential to do molecular work in remote locations.

Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome-wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome-wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome-wide associations identified known (Twin Sister of FT, FRIGIDA-like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage-Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.

Next-generation conservation genetics and biodiversity monitoring

This special issue of Evolutionary Applications consists of 10 publications investigating the use of next-generation tools and techniques in population genetic analyses and biodiversity assessment. The special issue stems from a 2016 Next Generation Genetic Monitoring Workshop, hosted by the National Institute for Mathematical and Biological Synthesis (NIMBioS) in Tennessee, USA. The improved accessibility of next-generation sequencing platforms has allowed molecular ecologists to rapidly produce large amounts of data. However, with the increased availability of new genomic markers and mathematical techniques, care is needed in selecting appropriate study designs, interpreting results in light of conservation concerns, and determining appropriate management actions. This special issue identifies key attributes of successful genetic data analyses in biodiversity evaluation and suggests ways to improve analyses and their application in current population and conservation genetics research.

Measuring Genetic Differentiation from Pool-seq Data

The advent of high throughput sequencing and genotyping technologies enables the comparison of patterns of polymorphisms at a very large number of markers. While the characterization of genetic structure from individual sequencing data remains expensive for many nonmodel species, it has been shown that sequencing pools of individual DNAs (Pool-seq) represents an attractive and cost-effective alternative. However, analyzing sequence read counts from a DNA pool instead of individual genotypes raises statistical challenges in deriving correct estimates of genetic differentiation. In this article, we provide a method-of-moments estimator of [Formula: see text] for Pool-seq data, based on an analysis-of-variance framework. We show, by means of simulations, that this new estimator is unbiased and outperforms previously proposed estimators. We evaluate the robustness of our estimator to model misspecification, such as sequencing errors and uneven contributions of individual DNAs to the pools. Finally, by reanalyzing published Pool-seq data of different ecotypes of the prickly sculpin Cottus asper, we show how the use of an unbiased [Formula: see text] estimator may question the interpretation of population structure inferred from previous analyses.

Valuing and understanding fish populations in the Anthropocene: key questions to address

Research on the values of fish populations and fisheries has primarily focused on bio-economic aspects; a more nuanced and multidimensional perspective is mostly neglected. Although a range of social aspects is increasingly being considered in fisheries research, there is still no clear understanding as to how to include these additional values within management policies nor is there a cogent appreciation of the major knowledge gaps that should be tackled by future research. This paper results from a workshop held during the 50th anniversary symposium of the Fisheries Society of the British Isles at the University of Exeter, UK, in July 2017. Here, we aim to highlight the current knowledge gaps on the values of fish populations and fisheries thus directing future research. To this end, we present eight questions that are deeply relevant to understanding the values of fish populations and fisheries. These can be applied to all habitats and fisheries, including freshwater, estuarine and marine.