Next Gen Pop Gen: implementing a high-throughput approach to population genetics in boarfish (Capros aper)

The genomic basis and environmental correlates of local adaptation in the Atlantic horse mackerel (Trachurus trachurus)

Understanding how populations adapt to their environment is increasingly important to prevent biodiversity loss due to overexploitation and climate change. Here we studied the population structure and genetic basis of local adaptation of Atlantic horse mackerel, a commercially and ecologically important marine fish that has one of the widest distributions in the eastern Atlantic. We analyzed whole-genome sequencing and environmental data of samples collected from the North Sea to North Africa and the western Mediterranean Sea. Our genomic approach indicated low population structure with a major split between the Mediterranean Sea and the Atlantic Ocean and between locations north and south of mid-Portugal. Populations from the North Sea are the most genetically distinct in the Atlantic. We discovered that most population structure patterns are driven by a few highly differentiated putatively adaptive loci. Seven loci discriminate the North Sea, two the Mediterranean Sea, and a large putative inversion (9.9 Mb) on chromosome 21 underlines the north-south divide and distinguishes North Africa. A genome-environment association analysis indicates that mean seawater temperature and temperature range, or factors correlated to them, are likely the main environmental drivers of local adaptation. Our genomic data broadly support the current stock divisions, but highlight areas of potential mixing, which require further investigation. Moreover, we demonstrate that as few as 17 highly informative SNPs can genetically discriminate the North Sea and North African samples from neighboring populations. Our study highlights the importance of both, life history and climate-related selective pressures in shaping population structure patterns in marine fish. It also supports that chromosomal rearrangements play a key role in local adaptation with gene flow. This study provides the basis for more accurate delineation of the horse mackerel stocks and paves the way for improving stock assessments.

A baseline for the genetic stock identification of Atlantic herring, Clupea harengus, in ICES Divisions 6.a, 7.b-c

Atlantic herring in International Council for Exploration of the Sea (ICES) Divisions 6.a, 7.b-c comprises at least three populations, distinguished by temporal and spatial differences in spawning, which have until recently been managed as two stocks defined by geographical delineators. Outside of spawning the populations form mixed aggregations, which are the subject of acoustic surveys. The inability to distinguish the populations has prevented the development of separate survey indices and separate stock assessments. A panel of 45 single-nucleotide polymorphisms, derived from whole-genome sequencing, were used to genotype 3480 baseline spawning samples (2014-2021). A temporally stable baseline comprising 2316 herring from populations known to inhabit Division 6.a was used to develop a genetic assignment method, with a self-assignment accuracy greater than 90%. The long-term temporal stability of the assignment model was validated by assigning archive (2003-2004) baseline samples (270 individuals) with a high level of accuracy. Assignment of non-baseline samples (1514 individuals) from Divisions 6.a, 7.b-c indicated previously unrecognized levels of mixing of populations outside of the spawning season. The genetic markers and assignment models presented constitute a 'toolbox' that can be used for the assignment of herring caught in mixed survey and commercial catches in Division 6.a into their population of origin with a high level of accuracy.

Dispersal ability, habitat characteristics, and sea-surface circulation shape population structure of Cingula trifasciata (Gastropoda: Rissoidae) in the remote Azores Archipelago

Background

In the marine realm, dispersal ability is among the major factors shaping the distribution of species. In the Northeast Atlantic Ocean, the Azores Archipelago is home to a multitude of marine invertebrates which, despite their dispersal limitations, maintain gene flow among distant populations, with complex evolutionary and biogeographic implications. The mechanisms and factors underlying the population dynamics and genetic structure of non-planktotrophic gastropods within the Azores Archipelago and related mainland populations are still poorly understood. The rissoid Cingula trifasciata is herewith studied to clarify its population structure in the Northeast Atlantic Ocean and factors shaping it, with a special focus in intra-archipelagic dynamics.

Results

Coupling microsatellite genotyping by amplicon sequencing (SSR-GBAS) and mitochondrial datasets, our results suggest the differentiation between insular and continental populations of Cingula trifasciata, supporting previously raised classification issues and detecting potential cryptic diversity. The finding of connectivity between widely separated populations was startling. In unique ways, dispersal ability, habitat type, and small-scale oceanographic currents appear to be the key drivers of C. trifasciata’s population structure in the remote Azores Archipelago. Dispersal as non-planktotrophic larvae is unlikely, but its small-size adults easily engage in rafting. Although the typical habitat of C. trifasciata, with low hydrodynamics, reduces the likelihood of rafting, individuals inhabiting algal mats are more prone to dispersal. Sea-surface circulation might create dispersal pathways for rafts, even between widely separated populations/islands.

Conclusions

Our results show that gene flow of a marine non-planktotrophic gastropod within a remote archipelago can reveal unanticipated patterns, such that the understanding of life in such areas is far from well-understood. We expect this work to be the starting of the application of SSR-GBAS in other non-model marine invertebrates, providing insights on their population dynamics at distinct geographical scales and on hidden diversity. How transversal is the role played by the complex interaction between functional traits, ecological features, and sea-surface circulation in the population structure of marine invertebrates can be further addressed by expanding this approach to more taxa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01862-1.

Evidence for multiple introductions of an invasive wild bee species currently under rapid range expansion in Europe

BACKGROUND

Invasive species are increasingly driving biodiversity decline, and knowledge of colonization dynamics, including both drivers and dispersal modes, are important to prevent future invasions. The bee species Megachile sculpturalis (Hymenoptera: Megachilidae), native to East-Asia, was first recognized in Southeast-France in 2008, and has since spread throughout much of Europe. The spread is very fast, and colonization may result from multiple fronts.

RESULT

To track the history of this invasion, codominant markers were genotyped using Illumina sequencing and the invasion history and degree of connectivity between populations across the European invasion axis were investigated. Distinctive genetic clusters were detected with east-west differentiations in Middle-Europe.

CONCLUSION

We hypothesize that the observed cluster formation resulted from multiple, independent introductions of the species to the European continent. This study draws a first picture of an early invasion stage of this wild bee and forms a foundation for further investigations, including studies of the species in their native Asian range and in the invaded range in North America.

Genetic diversity and genetic structure of Acer monspessulanum L. across Zagros forests of Iran using molecular markers

Acer monspessulanum L. is an important tree species found in the temperate Zagros forests of Iran. Despite its importance, the long-term persistence of its small and fragmented populations is jeopardised by genetic erosion and hence, monitoring its genetic resource and variability is practically required for providing conservation measures of the species germplasm in Zagros woodland ecosystem. The present study aimed to provide the first data on genetic diversity and genetic differentiation pattern of 19 natural populations comprising 240 individuals of A. monspessulanum across its growing area in Zagros forests using three molecular tools including inter-simple sequence repeats (ISSR), start codon targeted (SCoT), and simple sequence repeat (SSR). In total, ISSR and SCoT primers generated a total of 141 and 121 clear and scorable bands for analysis with the polymorphism rate of 90.50 and 90.02% and a mean of 10.85 and 11 fragments per marker, respectively. In addition, 73 alleles were achieved using 10 polymorphic SSR loci from the studied accessions with 100% polymorphism, ranging between 5 and 10 alleles per locus. Average percentage of polymorphic alleles per population for ISSR, SCoT, and SSR data was 84.02%, 83%, and 100%, respectively, and generally, Nei's gene diversity (H) and Shannon's index of diversity (I) values for all populations demonstrated moderate to high levels of genetic diversity (H = 0.267-0.707; I = 0.38-1.38). The AMOVA results of the three marker systems attributed higher genetic variation to individuals within in each population than among populations. Furthermore, overall G_ST value for all populations detected the moderate to high levels of genetic differentiation, indicating a limited gene flow occurrence among the populations. STRUCTURE analysis (K = 5) clustered the populations into four to five distinct groups, in accordance with geographical distances. These results could represent an important contribution for effective germplasm characterization and could be eventually used in in situ or ex situ conservation of A. monspessulanum genetic resources.

The Potential of HTS Approaches for Accurate Genotyping in Grapevine (Vitis vinifera L.)

The main challenge associated with genotyping based on conventional length polymorphisms is the cross-laboratory standardization of allele sizes. This step requires the inclusion of standards and manual sizing to avoid false results. Capillary electrophoresis (CE) approaches limit the information to the length polymorphism and do not allow the determination of a complete marker sequence. As an alternative, high-throughput sequencing (HTS) offers complete information regarding marker sequences and their flanking regions. In this work, we investigated the suitability of a semi-quantitative sequencing approach for microsatellite genotyping using Illumina paired-end technology. Twelve microsatellite loci that are well established for grapevine CE typing were analysed on 96 grapevine samples from six different countries. We redesigned primers to the length of the amplicon for short sequencing (~100 bp). The primer pair was flanked with a 10 bp overhang for the introduction of barcodes on both sides of the amplicon to enable high multiplexing. The highest data peaks were determined as simple sequence repeat (SSR) alleles and compared with the CE dataset based on 12 reference samples. The comparison showed that HTS SSR genotyping can successfully replace the CE system in further experiments. We believe that, with next-generation sequencing, genotyping can be improved in terms of its speed, accuracy, and price.

Fast sequence-based microsatellite genotyping development workflow

Application of high-throughput sequencing technologies to microsatellite genotyping (SSRseq) has been shown to remove many of the limitations of electrophoresis-based methods and to refine inference of population genetic diversity and structure. We present here a streamlined SSRseq development workflow that includes microsatellite development, multiplexed marker amplification and sequencing, and automated bioinformatics data analysis. We illustrate its application to five groups of species across phyla (fungi, plant, insect and fish) with different levels of genomic resource availability. We found that relying on previously developed microsatellite assay is not optimal and leads to a resulting low number of reliable locus being genotyped. In contrast, de novo ad hoc primer designs gives highly multiplexed microsatellite assays that can be sequenced to produce high quality genotypes for 20-40 loci. We highlight critical upfront development factors to consider for effective SSRseq setup in a wide range of situations. Sequence analysis accounting for all linked polymorphisms along the sequence quickly generates a powerful multi-allelic haplotype-based genotypic dataset, calling to new theoretical and analytical frameworks to extract more information from multi-nucleotide polymorphism marker systems.

Genome- and transcriptome-derived microsatellite loci in lumpfish Cyclopterus lumpus: molecular tools for aquaculture, conservation and fisheries management

The lumpfish Cyclopterus lumpus is commercially exploited in numerous areas of its range in the North Atlantic Ocean, and is important in salmonid aquaculture as a biological agent for controlling sea lice. Despite the economic importance, few genetic resources for downstream applications, such as linkage mapping, parentage analysis, marker-assisted selection (MAS), quantitative trait loci (QTL) analysis, and assessing adaptive genetic diversity are currently available for the species. Here, we identify both genome- and transcriptome-derived microsatellites loci from C. lumpus to facilitate such applications. Across 2,346 genomic contigs, we detected a total of 3,067 microsatellite loci, of which 723 were the most suitable ones for primer design. From 116,555 transcriptomic unigenes, we identified a total of 231,556 microsatellite loci, which may indicate a high coverage of the available STRs. Out of these, primer pairs could only be designed for 6,203 loci. Dinucleotide repeats accounted for 89 percent and 52 percent of the genome- and transcriptome-derived microsatellites, respectively. The genetic composition of the dominant repeat motif types showed differences from other investigated fish species. In the genome-derived microsatellites AC/GT (67.8 percent), followed by AG/CT (15.1 percent) and AT/AT (5.6 percent) were the major motifs. Transcriptome-derived microsatellites showed also most dominantly the AC/GT repeat motif (33 percent), followed by A/T (26.6 percent) and AG/CT (11 percent). Functional annotation of microsatellite-containing transcriptomic sequences showed that the majority of the expressed sequence tags encode proteins involved in cellular and metabolic processes, binding activity and catalytic reactions. Importantly, STRs linked to genes involved in immune system process, growth, locomotion and reproduction were discovered in the present study. The extensive genomic marker information reported here will facilitate molecular ecology studies, conservation initiatives and will benefit many aspects of the breeding programmes of C. lumpus.

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

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

CHIIMP: An automated high-throughput microsatellite genotyping platform reveals greater allelic diversity in wild chimpanzees

Short tandem repeats (STRs), also known as microsatellites, are commonly used to noninvasively genotype wild-living endangered species, including African apes. Until recently, capillary electrophoresis has been the method of choice to determine the length of polymorphic STR loci. However, this technique is labor intensive, difficult to compare across platforms, and notoriously imprecise. Here we developed a MiSeq-based approach and tested its performance using previously genotyped fecal samples from long-term studied chimpanzees in Gombe National Park, Tanzania. Using data from eight microsatellite loci as a reference, we designed a bioinformatics platform that converts raw MiSeq reads into locus-specific files and automatically calls alleles after filtering stutter sequences and other PCR artifacts. Applying this method to the entire Gombe population, we confirmed previously reported genotypes, but also identified 31 new alleles that had been missed due to sequence differences and size homoplasy. The new genotypes, which increased the allelic diversity and heterozygosity in Gombe by 61% and 8%, respectively, were validated by replicate amplification and pedigree analyses. This demonstrated inheritance and resolved one case of an ambiguous paternity. Using both singleplex and multiplex locus amplification, we also genotyped fecal samples from chimpanzees in the Greater Mahale Ecosystem in Tanzania, demonstrating the utility of the MiSeq-based approach for genotyping nonhabituated populations and performing comparative analyses across field sites. The new automated high-throughput analysis platform (available at https://github.com/ShawHahnLab/chiimp) will allow biologists to more accurately and effectively determine wildlife population size and structure, and thus obtain information critical for conservation efforts.