Genome-wide SNPs reveal low effective population size within confined management units of the highly vagile Galapagos shark (Carcharhinus galapagensis)

Estimation of effective number of breeders and effective population size in an abundant and heavily exploited marine teleost

Obtaining reliable estimates of the effective number of breeders (N b) and generational effective population size (N e) for fishery-important species is challenging because they are often iteroparous and highly abundant, which can lead to bias and imprecision. However, recent advances in understanding of these parameters, as well as the development of bias correction methods, have improved the capacity to generate reliable estimates. We utilized samples of both single-cohort young of the year and mixed-age adults from two geographically and genetically isolated stocks of the Australasian snapper (Chrysophrys auratus) to investigate the feasibility of generating reliable N b and N e estimates for a fishery species. Snapper is an abundant, iteroparous broadcast spawning teleost that is heavily exploited by recreational and commercial fisheries. Employing neutral genome-wide SNPs and the linkage-disequilibrium method, we determined that the most reliable N b and N e estimates could be derived by genotyping at least 200 individuals from a single cohort. Although our estimates made from the mixed-age adult samples were generally lower and less precise than those based on a single cohort, they still proved useful for understanding relative differences in genetic effective size between stocks. The correction formulas applied to adjust for biases due to physical linkage of loci and age structure resulted in substantial upward modifications of our estimates, demonstrating the importance of applying these bias corrections. Our findings provide important guidelines for estimating N b and N e for iteroparous species with large populations. This work also highlights the utility of samples originally collected for stock structure and stock assessment work for investigating genetic effective size in fishery-important species.

Translating genomic advances into biodiversity conservation

A key action of the new Global Biodiversity Framework is the maintenance of genetic diversity in all species to safeguard their adaptive potential. To achieve this goal, a translational mindset, which aims to convert results of basic research into direct practical benefits, needs to be applied to biodiversity conservation. Despite much discussion on the value of genomics to conservation, a disconnect between those generating genomic resources and those applying it to biodiversity management remains. As global efforts to generate reference genomes for non-model species increase, investment into practical biodiversity applications is critically important. Applications such as understanding population and multispecies diversity and longitudinal monitoring need support alongside education for policymakers on integrating the data into evidence-based decisions. Without such investment, the opportunity to revolutionize global biodiversity conservation using genomics will not be fully realized.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

Genomic evidence indicates small island-resident populations and sex-biased behaviors of Hawaiian reef Manta Rays

BACKGROUND

Reef manta rays (Mobula alfredi) are globally distributed in tropical and subtropical seas. Their life history traits (slow growth, late maturity, low reproductive output) make them vulnerable to perturbations and therefore require informed management strategies. Previous studies have reported wide-spread genetic connectivity along continental shelves suggesting high gene flow along continuous habitats spanning hundreds of kilometers. However, in the Hawaiian Islands, tagging and photo-identification evidence suggest island populations are isolated despite proximity, a hypothesis that has not yet been evaluated with genetic data.

RESULTS

This island-resident hypothesis was tested by analyzing whole mitogenome haplotypes and 2048 nuclear single nucleotide polymorphisms (SNPs) between M. alfredi (n = 38) on Hawai'i Island and Maui Nui (the 4-island complex of Maui, Moloka'i, Lāna'i and Kaho'olawe). Strong divergence in the mitogenome (ΦST = 0.488) relative to nuclear genome-wide SNPs (neutral FST = 0.003; outlier FST = 0.186), and clustering of mitochondrial haplotypes among islands provides robust evidence that female reef manta rays are strongly philopatric and do not migrate between these two island groups. Combined with restricted male-mediated migration, equivalent to a single male moving between islands every 2.2 generations (~ 64 years), we provide evidence these populations are significantly demographically isolated. Estimates of contemporary effective population size (Ne) are 104 (95% CI: 99-110) in Hawai'i Island and 129 (95% CI: 122-136) in Maui Nui.

CONCLUSIONS

Concordant with evidence from photo identification and tagging studies, these genetic results indicate reef manta rays in Hawai'i have small, genetically-isolated resident island populations. We hypothesize that due to the Island Mass Effect, large islands provide sufficient resources to support resident populations, thereby making crossing deep channels separating island groups unnecessary. Small effective population size, low genetic diversity, and k-selected life history traits make these isolated populations vulnerable to region-specific anthropogenic threats, which include entanglement, boat strikes, and habitat degradation. The long-term persistence of reef manta rays in the Hawaiian Islands will require island-specific management strategies.

The genetics-morphology-behavior trifecta: Unraveling the single greatest limitation affecting our understanding of chondrichthyan evolution

Sharks, rays, and chimaera form the clade Chondrichthyes, an ancient group of morphologically and ecologically diverse vertebrates that has played an important role in our understanding of gnathostome evolution. Increasingly, studies seek to investigate evolutionary processes operating within the chondrichthyan crown group, with the broad aim of understanding the driving forces behind the vast phenotypic diversity observed among its constituent taxa. Genetic, morphological, and behavioral studies have all contributed to our understanding of phenotypic evolution yet are typically considered in isolation in the context of Chondrichthyes. In this viewpoint, I discuss why such isolation is prevalent in the literature, how it constrains our understanding of evolution, and how it might be overcome. I argue that integrating these core fields of organismal biology is vital if we are to understand the evolutionary processes operating in contemporary chondrichthyan taxa and how such processes have contributed to past phenotypic evolution. Despite this, the necessary tools to overcome this major limitation already exist and have been applied to other taxa.

Stepping up to genome scan allows stock differentiation in the worldwide distributed blue shark Prionace glauca

The blue shark Prionace glauca is a top predator with one of the widest geographical distributions of any shark species. It is classified as Critically Endangered in the Mediterranean Sea, and Near Threatened globally. Previous genetic studies did not reject the null hypothesis of a single global population. The blue shark was proposed as a possible archetype of the "grey zone of population differentiation," coined to designate cases where population structure may be too recent or too faint to be detected using a limited set of markers. Here, blue shark samples collected throughout its global range were sequenced using a specific RAD method (DArTseq), which recovered 37,655 genome-wide single nucleotide polymorphisms (SNPs). Two main groups emerged, with Mediterranean Sea and northern Atlantic samples (Northern population) differentiated significantly from the Indo-west Pacific samples (Southern population). Significant pairwise F_ST values indicated further genetic differentiation within the Atlantic Ocean, and between the Atlantic Ocean and the Mediterranean Sea. Reconstruction of recent demographic history suggested divergence between Northern and Southern populations occurred about 500 generations ago and revealed a drastic reduction in effective population size from a large ancestral population. Our results illustrate the power of genome scans to detect population structure and reconstruct demographic history in highly migratory marine species. Given that the management plans of the blue shark (targeted or bycatch) fisheries currently assume panmictic regional stocks, we strongly recommend that the results presented here be considered in future stock assessments and conservation strategies.

Genetic structure and relatedness of juvenile sicklefin lemon shark (Negaprion acutidens) at Dongsha Island

Negaprion acutidens (sicklefin lemon shark) is distributed in the Indo-Pacific and in close association with coral reefs. Under the protection of the Dongsha Atoll National Park, a small but well-established juvenile population of N. acutidens inhabiting coastal areas of Dongsha Island was recently observed to display site fidelity by using acoustic telemetry. This study was designed to reveal the fine scale genetic structure and relatedness within and among 5 juvenile shark cohorts inhabiting 3 sampling sites at Dongsha Island. A total 188 juveniles were caught and sampled between 2016 and 2017, and genotyped with twelve loci. They were assigned to 5 year cohorts (2013-2017) based on the body length and date they were caught, also assigned to 3 sites based on where they were caught. Among five cohorts, the percentage of unrelated pairs within a cohort is more than 62% in average, suggesting a potential high mortality during their early life stage. The results of Fst and assignment testing showed that there was no significant genetic structure between sites and cohorts indicating that there was no fine scale genetic structure, even though the juveniles possessed strong site fidelity. A small effective population size (Ne) was detected (Ne = 86.7) which indicates the presence of a potentially isolated and vulnerable population at Dongsha. These results provide the genetic diversity as a baseline for future management and conservation of N. acutidens in the South China Sea.

Assessing population genetic structure of three New Zealand stream insects using mitochondrial and nuclear DNA markers

Assessing genetic differentiation among natural populations can aid understanding of dispersal patterns and connectivity among habitats. Several molecular markers have become increasingly popular in determining population genetic structure for this purpose. Here, we compared the resolution of mitochondrial cytochrome c oxidase subunit I (COI) and nuclear single nucleotide polymorphism (SNP) markers for detecting population structure among stream insects at small spatial scales. Individuals of three endemic taxa - Coloburiscus humeralis (Ephemeroptera), Zelandobius confusus (Plecoptera), and Hydropsyche fimbriata (Trichoptera) - were collected from forested streams that flow across open pasture in the North Island of New Zealand. Both COI and SNP data indicated limited population structure across the study area, and small differences observed among these species were likely related to their putative dispersal abilities. For example, fine-scale genetic differentiation between and among neighbouring stream populations for H. fimbriata suggests that gene flow, and hence dispersal, may be more limited for this species relative to the others. Based on the generally similar results provided by both types of markers, we suggest that either COI or SNP markers can provide suitable initial estimates of fine-scale population genetic differentiation in stream insects.

Genomic time-series data show that gene flow maintains high genetic diversity despite substantial genetic drift in a butterfly species

Effective population size affects the efficacy of selection, rate of evolution by drift, and neutral diversity levels. When species are subdivided into multiple populations connected by gene flow, evolutionary processes can depend on global or local effective population sizes. Theory predicts that high levels of diversity might be maintained by gene flow, even very low levels of gene flow, consistent with species long-term effective population size, but tests of this idea are mostly lacking. Here, we show that Lycaeides buttery populations maintain low contemporary (variance) effective population sizes (e.g., ~200 individuals) and thus evolve rapidly by genetic drift. In contrast, populations harbored high levels of genetic diversity consistent with an effective population size several orders of magnitude larger. We hypothesized that the differences in the magnitude and variability of contemporary versus long-term effective population sizes were caused by gene flow of sufficient magnitude to maintain diversity but only subtly affect evolution on generational time scales. Consistent with this hypothesis, we detected low but non-trivial gene flow among populations. Furthermore, using short-term population-genomic time-series data, we documented patterns consistent with predictions from this hypothesis, including a weak but detectable excess of evolutionary change in the direction of the mean (migrant gene pool) allele frequencies across populations, and consistency in the direction of allele frequency change over time. The documented decoupling of diversity levels and short-term change by drift in Lycaeides has implications for our understanding of contemporary evolution and the maintenance of genetic variation in the wild.

Green, yellow or black? Genetic differentiation and adaptation signatures in a highly migratory marine turtle

Marine species may exhibit genetic structure accompanied by phenotypic differentiation related to adaptation despite their high mobility. Two shape-based morphotypes have been identified for the green turtle (Chelonia mydas) in the Pacific Ocean: the south-central/western or yellow turtle and north-central/eastern or black turtle. The genetic differentiation between these morphotypes and the adaptation of the black turtle to environmentally contrasting conditions of the eastern Pacific region has remained a mystery for decades. Here we addressed both questions using a reduced-representation genome approach (Dartseq; 9473 neutral SNPs) and identifying candidate outlier loci (67 outlier SNPs) of biological relevance between shape-based morphotypes from eight Pacific foraging grounds (n = 158). Our results support genetic divergence between morphotypes, probably arising from strong natal homing behaviour. Genes and enriched biological functions linked to thermoregulation, hypoxia, melanism, morphogenesis, osmoregulation, diet and reproduction were found to be outliers for differentiation, providing evidence for adaptation of C. mydas to the eastern Pacific region and suggesting independent evolutionary trajectories of the shape-based morphotypes. Our findings support the evolutionary distinctness of the enigmatic black turtle and contribute to the adaptive research and conservation genomics of a long-lived and highly mobile vertebrate.

A globally threatened shark, Carcharias taurus, shows no population decline in South Africa

Knowledge about the demographic histories of natural populations helps to evaluate their conservation status, and potential impacts of natural and anthropogenic pressures. In particular, estimates of effective population size obtained through molecular data can provide useful information to guide management decisions for vulnerable populations. The spotted ragged-tooth shark, Carcharias taurus (also known as the sandtiger or grey nurse shark), is widely distributed in warm-temperate and subtropical waters, but has suffered severe population declines across much of its range as a result of overexploitation. Here, we used multilocus genotype data to investigate the demographic history of the South African C. taurus population. Using approximate Bayesian computation and likelihood-based importance sampling, we found that the population underwent a historical range expansion that may have been linked to climatic changes during the late Pleistocene. There was no evidence for a recent anthropogenic decline. Together with census data suggesting a stable population, these results support the idea that fishing pressure and other threats have so far not been detrimental to the local C. taurus population. The results reported here indicate that South Africa could possibly harbour the last remaining, relatively pristine population of this widespread but vulnerable top predator.

Extreme seascape drives local recruitment and genetic divergence in brooding and spawning corals in remote north-west Australia

Management strategies designed to conserve coral reefs threatened by climate change need to incorporate knowledge of the spatial distribution of inter- and intra-specific genetic diversity. We characterized patterns of genetic diversity and connectivity using single nucleotide polymorphisms (SNPs) in two reef-building corals to explore the eco-evolutionary processes that sustain populations in north-west Australia. Our sampling focused on the unique reefs of the Kimberley; we collected the broadcast spawning coral Acropora aspera (n = 534) and the brooding coral Isopora brueggemanni (n = 612) across inter-archipelago (tens to hundreds of kilometres), inter-reef (kilometres to tens of kilometres) and within-reef (tens of metres to a few kilometres) scales. Initial analysis of A. aspera identified four highly divergent lineages that were co-occurring but morphologically similar. Subsequent population analyses focused on the most abundant and widespread lineage, Acropora asp-c. Although the overall level of geographic subdivision was greater in the brooder than in the spawner, fundamental similarities in patterns of genetic structure were evident. Most notably, limits to gene flow were observed at scales <35 kilometres. Further, we observed four discrete clusters and a semi-permeable barrier to dispersal that were geographically consistent between species. Finally, sites experiencing bigger tides were more connected to the metapopulation and had greater gene diversity than those experiencing smaller tides. Our data indicate that the inshore reefs of the Kimberley are genetically isolated from neighbouring oceanic bioregions, but occasional dispersal between inshore archipelagos is important for the redistribution of evolutionarily important genetic diversity. Additionally, these results suggest that networks of marine reserves that effectively protect reefs from local pressures should be spaced within a few tens of kilometres to conserve the existing patterns of demographic and genetic connectivity.

Estimating effective population size using RADseq: Effects of SNP selection and sample size

Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.

Fine-scale species distribution modelling and genotyping by sequencing to examine hybridisation between two narrow endemic plant species

Hybridization has an important and often positive role in plant evolution. However, it can also have negative consequences for species. Two closely related species of Ornduffia are endemic to the Porongurup Range in the South West Australian Global Biodiversity Hotspot. The rare Ornduffia calthifolia is found exclusively on the summits, while O. marchantii is more widely dispersed across a greater range of elevation and is not considered threatened. Hybridisation in suitable overlapping habitat has been suspected between them for decades. Here we combine genotyping by sequencing to verify hybridisation genetically, and fine scale (2 m resolution) species distribution modelling (SDM) to test if hybrids occur in suitable intersecting habitat. From a study area of c. 4700 ha, SDM identified c. 275 ha and c. 322 ha of suitable habitat for O. calthifolia and O. marchantii, respectively. We identified range overlap between species of c. 59 ha), which enveloped 32 individuals confirmed to be hybrids. While the hybrids were at the margin of suitable habitat for O. marchantii, their preference for elevated habitat was closer to the more narrowly distributed O. calthifolia. The combination of genetic data and fine scale spatial modelling approaches enabled a better understanding of hybridisation among taxa of conservation significance. However, the level to which hybrid proliferation and competition for habitat presents as a threat to O. calthifolia is currently unknown and requires priority in conservation management given the threats from global warming and disturbance by tourism.

Genomics detects population structure within and between ocean basins in a circumpolar seabird: The white-chinned petrel

The Southern Ocean represents a continuous stretch of circumpolar marine habitat, but the potential physical and ecological drivers of evolutionary genetic differentiation across this vast ecosystem remain unclear. We tested for genetic structure across the full circumpolar range of the white-chinned petrel (Procellaria aequinoctialis) to unravel the potential drivers of population differentiation and test alternative population differentiation hypotheses. Following range-wide comprehensive sampling, we applied genomic (genotyping-by-sequencing or GBS; 60,709 loci) and standard mitochondrial-marker approaches (cytochrome b and first domain of control region) to quantify genetic diversity within and among island populations, test for isolation by distance, and quantify the number of genetic clusters using neutral and outlier (non-neutral) loci. Our results supported the multi-region hypothesis, with a range of analyses showing clear three-region genetic population structure, split by ocean basin, within two evolutionary units. The most significant differentiation between these regions confirmed previous work distinguishing New Zealand and nominate subspecies. Although there was little evidence of structure within the island groups of the Indian or Atlantic oceans, a small set of highly-discriminatory outlier loci could assign petrels to ocean basin and potentially to island group, though the latter needs further verification. Genomic data hold the key to revealing substantial regional genetic structure within wide-ranging circumpolar species previously assumed to be panmictic.

A matter of taste: Spatial and ontogenetic variations on the trophic ecology of the tiger shark at the Galapagos Marine Reserve

Sharks are top predators across ocean food webs and have a major ecological role in marine ecosystems. Investigating the trophic ecology of this group of species is thus essential to understand ecosystem functioning and inform specific management actions aimed at shark conservation. The Galapagos Islands represent one of the last ocean wildernesses, where populations of sharks and other top marine predators come close to a pristine status. Here we provide the first study on the trophic ecology of the tiger shark (Galeocerdo cuvier) within the Galapagos Marine Reserve (GMR), using a combination of stable isotope analysis, satellite tracking, and passive acoustic telemetry to investigate ontogenetic and spatial variations at two regions. The mean estimated δ¹³C and δ¹⁵N at Isabela island (western region) were -13.9 ± 0.5‰ and 13.7 ± 0.7‰; and for Santa Cruz island (central region) were -13.8 ± 0.3‰ and 13.4 ± 0.7‰, respectively. Green sea turtles (Chelonia mydas) were the main prey item for large tiger sharks (>280 cm TL), while smaller sharks mainly fed on squid and pelagic fish. Tiger sharks exhibited a high degree of philopatry around green sea-turtle nesting areas, with the majority of sharks detected around green sea-turtle nesting areas for at least 10 months after their capture date, and some individuals were even present during the entire three-year study period. Although we did not report statistically significant differences between the two regions, isotopic and electronic tagging data suggest that tiger sharks in the Galapagos could be segregated into specific populations separated by geographical scales of <100 km. The high productivity of the archipelago, along with the protection from industrial fishing granted by the GMR, result in abundant and predictable sources of prey. This high food abundance, combined with the presence of suitable habitats throughout the tiger shark life cycle, might result in a reduction of migratory behaviours when compared to movement patterns of tiger sharks in other ocean basins. Additional studies using genetic tools could provide further evidence on the presence of separate management units, as it has been recently revealed for other shark species inhabiting the GMR.

Museum specimens provide reliable SNP data for population genomic analysis of a widely distributed but threatened cockatoo species

Natural history museums harbour a plethora of biological specimens which are of potential use in population and conservation genetic studies. Although technical advancements in museum genomics have enabled genome-wide markers to be generated from aged museum specimens, the suitability of these data for robust biological inference is not well characterized. The aim of this study was to test the utility of museum specimens in population and conservation genomics by assessing the biological and technical validity of single nucleotide polymorphism (SNP) data derived from such samples. To achieve this, we generated thousands of SNPs from 47 red-tailed black cockatoo (Calyptorhychus banksii) traditional museum samples (i.e. samples that were not collected with the primary intent of DNA analysis) and 113 fresh tissue samples (cryopreserved liver/muscle) using a restriction site-associated DNA marker approach (DArTseq^™ ). Thousands of SNPs were successfully generated from most of the traditional museum samples (with a mean age of 44 years, ranging from 5 to 123 years), although 38% did not provide useful data. These SNPs exhibited higher error rates and contained significantly more missing data compared with SNPs from fresh tissue samples, likely due to considerable DNA fragmentation. However, based on simulation results, the level of genotyping error had a negligible effect on inference of population structure in this species. We did identify a bias towards low diversity SNPs in older samples that appears to compromise temporal inferences of genetic diversity. This study demonstrates the utility of a RADseq-based method to produce reliable genome-wide SNP data from traditional museum specimens.

Conservation genetics of elasmobranchs of the Mexican Pacific Coast, trends and perspectives

One of the most critical threats to biodiversity is the high extinction rate driven by human activities. Reducing extinction rates requires the implementation of conservation programmes based on robust scientific data. Elasmobranchs are important ecological components of the ocean, and several species sustain substantial economic activities. Unfortunately, elasmobranchs are one of the most threatened and understudied animal taxa. The Mexican Pacific Coast (MPC) is a region with high elasmobranch diversity and is the seat of major elasmobranch fisheries. But it is also a developing region with several conservation and management challenges which require national and international attention. Here, we review the conservation genetics literature of elasmobranchs from the MPC. We present a synthesis of the works using samples from the region and emphasize the main gaps and biases in these data. In addition, we discuss the benefits and challenges of generating genomic information to improve the management and conservation of an elasmobranch biodiversity hotspot in a developing country. We found 47 elasmobranch genetic articles that cover <30% of the elasmobranch diversity in the region. These studies mainly used mitochondrial DNA sequences to analyse the genetic structure of commercially important and abundant species of the order Carcharhiniformes. Some of these papers also assessed mating systems, demographic parameters, and taxonomic uncertainties, all of which are important topics for efficient management decisions. In terms of conservation genetics, elasmobranchs from the MPC remain understudied. However, high-throughput sequencing technologies have increased the power and accessibility of genomic tools, even in developing countries such as Mexico. The tools described here provide information relevant for biodiversity conservation. Therefore, we strongly suggest that investment in genomic research will assist implementation of efficient management strategies. In time, this will reduce the extinction risk of the unique elasmobranch biodiversity from the MPC.

Taking Advantage of the Genomics Revolution for Monitoring and Conservation of Chondrichthyan Populations

Chondrichthyes (sharks, rays, skates and chimaeras) are among the oldest extant predators and are vital to top-down regulation of oceanic ecosystems. They are an ecologically diverse group occupying a wide range of habitats and are thus, exploited by coastal, pelagic and deep-water fishing industries. Chondrichthyes are among the most data deficient vertebrate species groups making design and implementation of regulatory and conservation measures challenging. High-throughput sequencing technologies have significantly propelled ecological investigations and understanding of marine and terrestrial species’ populations, but there remains a paucity of NGS based research on chondrichthyan populations. We present a brief review of current methods to access genomic and metagenomic data from Chondrichthyes and discuss applications of these datasets to increase our understanding of chondrichthyan taxonomy, evolution, ecology and population structures. Last, we consider opportunities and challenges offered by genomic studies for conservation and management of chondrichthyan populations.

'Genome skimming' with the MinION hand-held sequencer identifies CITES-listed shark species in India's exports market

Chondrichthyes - sharks, rays, skates, and chimeras, are among the most threatened and data deficient vertebrate species. Global demand for shark and ray derived products, drives unregulated and exploitative fishing practices, which are in turn facilitated by the lack of ecological data required for effective conservation of these species. Here, we describe a Next Generation Sequencing method (using the MinION, a hand-held portable sequencing device from Oxford Nanopore Technologies), and analyses pipeline for molecular ecological studies in Chondrichthyes. Using this method, the complete mitochondrial genome and nuclear intergenic and protein-coding sequences were obtained by direct sequencing of genomic DNA obtained from shark fin tissue. Recovered loci include mitochondrial barcode sequences- Cytochrome oxidase I, NADH2, 16S rRNA and 12S rRNA- and nuclear genetic loci such as 5.8S rRNA, Internal Transcribed Spacer 2, and 28S rRNA regions, which are commonly used for taxonomic identification. Other loci recovered were the nuclear protein-coding genes for antithrombin or SerpinC, Immunoglobulin lambda light chain, Preprogehrelin, selenium binding protein 1(SBP1), Interleukin-1 beta (IL-1β) and Recombination-Activating Gene 1 (RAG1). The median coverage across all genetic loci was 20x and sequence accuracy was ≥99.8% compared to reference sequences. Analyses of the nuclear ITS2 region and the mitochondrial protein-encoding loci allowed accurate taxonomic identification of the shark specimen as Carcharhinus falciformis, a CITES Appendix II species. MinION sequencing provided 1,152,211 bp of new shark genome, increasing the number of sequenced shark genomes to five. Phylogenetic analyses using both mitochondrial and nuclear loci provided evidence that Prionace glauca is nested within Carcharhinus, suggesting the need for taxonomic reassignment of P. glauca. We increased genomic information about a shark species for ecological and population genetic studies, enabled accurate identification of the shark tissue for biodiversity indexing and resolved phylogenetic relationships among multiple taxa. The method was independent of amplification bias, and adaptable for field assessments of other Chondrichthyes and wildlife species in the future.

Cost-effective assembly of the African wild dog (Lycaon pictus) genome using linked reads

Background

A high-quality reference genome assembly is a valuable tool for the study of non-model organisms. Genomic techniques can provide important insights about past population sizes and local adaptation and can aid in the development of breeding management plans. This information is important for fields such as conservation genetics, where endangered species require critical and immediate attention. However, funding for genomic-based methods can be sparse for conservation projects, as costs for general species management can consume budgets.

Findings

Here, we report the generation of high-quality reference genomes for the African wild dog (Lycaon pictus) at a low cost (<$3000), thereby facilitating future studies of this endangered canid. We generated assemblies for three individuals using the linked-read 10x Genomics Chromium system. The most continuous assembly had a scaffold and contig N50 of 21 Mb and 83 Kb, respectively, and completely reconstructed 95% of a set of conserved mammalian genes. Additionally, we estimate the heterozygosity and demographic history of African wild dogs, revealing that although they have historically low effective population sizes, heterozygosity remains high.

Conclusions

We show that 10x Genomics Chromium data can be used to effectively generate high-quality genomes from Illumina short-read data of intermediate coverage (∼25x–50x). Interestingly, the wild dog shows higher heterozygosity than other species of conservation concern, possibly due to its behavioral ecology. The availability of reference genomes for non-model organisms will facilitate better genetic monitoring of threatened species such as the African wild dog and help conservationists to better understand the ecology and adaptability of those species in a changing environment.

A genome-wide approach for uncovering evolutionary relationships of Australian Bactrocera species complexes (Diptera: Tephritidae)

Australia and Southeast Asia are hotspots of global diversity in the fruit-fly genus Bactrocera. Although a great diversity of species has been long recognised, evolutionary relationships are poorly understood, largely because previous sequencing techniques have provided insufficient historical signal for phylogenetic reconstruction. Poorly understood biogeographic history in Bactrocera has prevented a deeper understanding of migratory patterns in this economically important pest group. Using representatives from Australia and Malaysia, we tested the utility of a genome-reduction approach that generates thousands of single-nucleotide polymorphisms for phylogenetic reconstructions. This approach has high utility for species identification because of the ease of sample addition over time, and the species-level specificity able to be achieved with the markers. These data have provided a strongly supported phylogenetic tree congruent with topologies generated using more intensive sequencing approaches. In addition, our results do not support taxonomic assignments to species complex for a number of species, such as B. endiandrae in the dorsalis complex, yet find a close relationship between B. pallida and the dorsalis species. Our data have further validated non-monophyletic evolution of male response to primary attractants. We also showed at least two diversification events between Australia and Southeast Asia, indicating trans-regional dispersal in important pest species.

Next-generation metrics for monitoring genetic erosion within populations of conservation concern

Genetic erosion is a major threat to biodiversity because it can reduce fitness and ultimately contribute to the extinction of populations. Here, we explore the use of quantitative metrics to detect and monitor genetic erosion. Monitoring systems should not only characterize the mechanisms and drivers of genetic erosion (inbreeding, genetic drift, demographic instability, population fragmentation, introgressive hybridization, selection) but also its consequences (inbreeding and outbreeding depression, emergence of large-effect detrimental alleles, maladaptation and loss of adaptability). Technological advances in genomics now allow the production of data the can be measured by new metrics with improved precision, increased efficiency and the potential to discriminate between neutral diversity (shaped mainly by population size and gene flow) and functional/adaptive diversity (shaped mainly by selection), allowing the assessment of management-relevant genetic markers. The requirements of such studies in terms of sample size and marker density largely depend on the kind of population monitored, the questions to be answered and the metrics employed. We discuss prospects for the integration of this new information and metrics into conservation monitoring programmes.

Strong trans-Pacific break and local conservation units in the Galapagos shark (Carcharhinus galapagensis) revealed by genome-wide cytonuclear markers

The application of genome-wide cytonuclear molecular data to identify management and adaptive units at various spatio-temporal levels is particularly important for overharvested large predatory organisms, often characterized by smaller, localized populations. Despite being "near threatened", current understanding of habitat use and population structure of Carcharhinus galapagensis is limited to specific areas within its distribution. We evaluated population structure and connectivity across the Pacific Ocean using genome-wide single-nucleotide polymorphisms (~7200 SNPs) and mitochondrial control region sequences (945 bp) for 229 individuals. Neutral SNPs defined at least two genetically discrete geographic groups: an East Tropical Pacific (Mexico, east and west Galapagos Islands), and another central-west Pacific (Lord Howe Island, Middleton Reef, Norfolk Island, Elizabeth Reef, Kermadec, Hawaii and Southern Africa). More fine-grade population structure was suggested using outlier SNPs: west Pacific, Hawaii, Mexico, and Galapagos. Consistently, mtDNA pairwise ΦST defined three regional stocks: east, central and west Pacific. Compared to neutral SNPs (FST = 0.023-0.035), mtDNA exhibited more divergence (ΦST = 0.258-0.539) and high overall genetic diversity (h = 0.794 ± 0.014; π = 0.004 ± 0.000), consistent with the longstanding eastern Pacific barrier between the east and central-west Pacific. Hawaiian and Southern African populations group within the west Pacific cluster. Effective population sizes were moderate/high for east/west populations (738 and 3421, respectively). Insights into the biology, connectivity, genetic diversity, and population demographics informs for improved conservation of this species, by delineating three to four conservation units across their Pacific distribution. Implementing such conservation management may be challenging, but is necessary to achieve long-term population resilience at basin and regional scales.