Application of SNPs for population genetics of nonmodel organisms: new opportunities and challenges

Collecting baleen whale blow samples by drone: A minimally intrusive tool for conservation genetics

In coastal British Columbia, Canada, marine megafauna such as humpback whales (Megaptera novaeangliae) and fin whales (Balaenoptera physalus velifera) have been subject to a history of exploitation and near extirpation. While their populations have been in recovery, significant threats are posed to these vulnerable species by proposed natural resource ventures in this region, in addition to the compounding effects of anthropogenic climate change. Genetic tools play a vital role in informing conservation efforts, but the associated collection of tissue biopsy samples can be challenging for the investigators and disruptive to the ongoing behaviour of the targeted whales. Here, we evaluate a minimally intrusive approach based on collecting exhaled breath condensate, or respiratory 'blow' samples, from baleen whales using an unoccupied aerial system (UAS), within Gitga'at First Nation territory for conservation genetics. Minimal behavioural responses to the sampling technique were observed, with no response detected 87% of the time (of 112 UAS deployments). DNA from whale blow (n = 88 samples) was extracted, and DNA profiles consisting of 10 nuclear microsatellite loci, sex identification and mitochondrial (mt) DNA haplotypes were constructed. An average of 7.5 microsatellite loci per individual were successfully genotyped. The success rates for mtDNA and sex assignment were 80% and 89% respectively. Thus, this minimally intrusive sampling method can be used to describe genetic diversity and generate genetic profiles for individual identification. The results of this research demonstrate the potential of UAS-collected whale blow for conservation genetics from a remote location.

Genetic characterization of cashmere goat (Capra hircus) populations in Mongolia

Objective

Characterization studies of the phenotypic and genetic diversity of Mongolian goats are limited, despite several goat breeds being registered in the country. This study aimed to evaluate the phenotypic and genetic diversity of 14 cashmere goat populations in Mongolia, consisting largely of identified goat breeds.

Methods

Body weight, cashmere quality, and coat color were the phenotypic traits considered in this study. A linear model was used to fit body weight and cashmere traits, and least squares means (LSMs) were estimated for the region and location classes. Genetic diversity and structure were assessed using a goat 50K SNP array.

Results

The studied populations exhibited greater phenotypic diversity at the regional level. A very small overall differentiation index (Fst: 0.017) was revealed by Wright's Fst and a very small overall inbreeding index (F ROH1 :0.019) was revealed based on runs of homozygosity. Genetic clustering of populations by principal components showed large variances for the two goat populations of the Russian admixture (Gobi Gurvan Saikhan and Uuliin Bor), and smaller but differentiated clusters for the remaining populations. Similar results were observed in the admixture analysis, which identified populations with the highest (Govi Gurvan Saikhan and Uuliin Bor) and lowest (Tsagaan Ovoo Khar) exotic admixtures. A genomewide association study (GWAS) of body weight and cashmere traits identified a few significant variants on chromosomes 2, 4, 5, 9, and 15, with the strongest variant for cashmere yield on chromosome 4. The GWAS on coat color yielded nine significant variants, with the strongest variants located on chromosomes 6, 13, and 18 and potential associations with KIT, ASIP, and MC1R genes. These signals were also found in other studies on coat color and patterns in goats.

Conclusion

Mongolian cashmere goats showed relatively low genetic differentiation and low inbreeding levels, possibly caused by the traditional pastoral livestock management system and the practice of trading breeding bucks across provinces, along with a recent increase in the goat population. Further investigation of cashmere traits using larger samples and alternative methods may help identify the genes or genomic regions underlying cashmere quality in goats.

New insights into mycotoxin risk management through fungal population genetics and genomics

Mycotoxin contamination of food and feed is a major global concern. Chronic or acute dietary exposure to contaminated food and feed can negatively affect both human and animal health. Contamination occurs through plant infection by toxigenic fungi, primarily Aspergillus and Fusarium spp., either before or after harvest. Despite the application of various management strategies, controlling these pathogens remains a major challenge primarily because of their ability to adapt to environmental changes and selection pressures. Understanding the genetic structure of plant pathogen populations is pivotal for gaining new insights into their biology and epidemiology, as well as for understanding the mechanisms behind their adaptability. Such deeper understanding is crucial for developing effective and preemptive management strategies tailored to the evolving nature of pathogenic populations. This review focuses on the population-level variations within the two most economically significant toxigenic fungal genera according to space, host, and pathogenicity. Outcomes in terms of migration patterns, gene flow within populations, mating abilities, and the potential for host jumps are examined. We also discuss effective yet often underutilized applications of population genetics and genomics to address practical challenges in the epidemiology and disease control of toxigenic fungi.

RADseq datasets of native beans from Mexico

Forty-five accessions of the genus Phaseolus from the orthodox seed collection of the National Center for Genetic Resources (CNRG) of the National Institute of Forestry, Agricultural, and Livestock Research (INIFAP) of Mexico were sequenced using RADseq. The species utilized were: P. acutifolius (14), P. coccineus (12), P. lunatus (8), P. dumosus (6), P. leptostachyus (2), P. filiformis (2), and P. vulgaris (1). A variant call file (VCF) was generated using GATK with the P. vulgaris reference genome GCF_000499845.1, identifying 97,103 shared SNPs among the species. These data have the potential to be used for studies of genetic diversity intra and interspecies, phylogeny, evolution, genetic resource conservation, and agricultural improvement.

Internal transcribed spacer as effective molecular marker for the detection of natural hybridization between the bivalves Pinna nobilis and Pinna rudis

The Pinna nobilis, a Mediterranean mollusc, has suffered population declines due to a massive mortality event associated with various factors including the parasite Haplosporidium pinnae. Some populations show resilience, possibly due to local environmental conditions. In this study, a molecular multiplex PCR method was developed using species-specific primers targeting Internal Transcribed Spacer (ITS) regions of P. nobilis and P. rudis, allowing accurate species identification and hybrid detection. Samples from Mediterranean areas were analysed, including putative hybrids and individuals from five other bivalve species. DNA was isolated, ITS regions were amplified and sequenced, and phylogenetic analyses confirmed species differentiation and primer specificity. The multiplex-PCR successfully identified P. nobilis, P. rudis, and their hybrids based on distinct amplicon patterns. This study highlights the value of molecular tools in species conservation, especially for monitoring and managing hybridization, supporting effective biodiversity conservation strategies.

Genetic diversity and population structure of Piper nigrum (black pepper) accessions based on next-generation SNP markers

Despite the economic importance of Piper nigrum (black pepper), a highly valued crop worldwide, development and utilization of genomic resources have remained limited, with diversity assessments often relying on only a few samples or DNA markers. Here we employed restriction-site associated DNA sequencing to analyze 175 P. nigrum accessions from eight main black pepper growing regions in Sri Lanka. The sequencing effort resulted in 1,976 million raw reads, averaging 11.3 million reads per accession, revealing 150,356 high-quality single nucleotide polymorphisms (SNPs) distributed across 26 chromosomes. Population structure analysis revealed two subpopulations (K = 2): a dominant group consisting of 152 accessions sourced from both home gardens and large-scale cultivations, and a smaller group comprising 23 accessions exclusively from native collections in home gardens. This clustering was further supported by principal component analysis, with the first two principal components explaining 35.2 and 12.1% of the total variation. Genetic diversity analysis indicated substantial gene flow (Nm = 342.21) and a low fixation index (FST = 0.00073) between the two subpopulations, with no clear genetic differentiation among accessions from different agro-climatic regions. These findings demonstrate that most current black pepper genotypes grown in Sri Lanka share a common genetic background, emphasizing the necessity to broaden the genetic base to enhance resilience to biotic and abiotic stresses. This study represents the first attempt at analyzing black pepper genetic diversity using high-resolution SNP markers, laying the foundation for future genome-wide association studies for SNP-based gene discovery and breeding.

Evaluation of DNA yield from various tissue and sampling sources for use in single nucleotide polymorphism panels

Genetics studies are used by wildlife managers and researchers to gain inference into a population of a species of interest. To gain these insights, microsatellites have been the primary method; however, there currently is a shift from microsatellites to single nucleotide polymorphisms (SNPs). With the different DNA requirements between microsatellites and SNPs, an investigation into which samples can provide adequate DNA yield is warranted. Using samples that were collected from previous genetic projects from regions in the USA from 2014 to 2021, we investigated the DNA yield of eight sample categories to gain insights into which provided adequate DNA to be used in ddRADseq or already developed high- or medium-density SNP panels. We found seven sample categories that met the DNA requirements for use in all three panels, and one sample category that did not meet any of the three panels requirements; however, DNA integrity was highly variable and not all sample categories that met panel DNA requirements could be considered high quality DNA. Additionally, we used linear random-effects models to determine which covariates would have the greatest influence on DNA yield. We determined that all covariates (tissue type, storage method, preservative, DNA quality, time until DNA extraction and time after DNA extraction) could influence DNA yield.

Global and Local Ancestry and its Importance: A Review

The fastest way to significantly change the composition of a population is through admixture, an evolutionary mechanism. In animal breeding history, genetic admixture has provided both short-term and long-term advantages by utilizing the phenomenon of complementarity and heterosis in several traits and genetic diversity, respectively. The traditional method of admixture analysis by pedigree records has now been replaced greatly by genome-wide marker data that enables more precise estimations. Among these markers, SNPs have been the popular choice since they are cost-effective, not so laborious, and automation of genotyping is easy. Certain markers can suggest the possibility of a population's origin from a sample of DNA where the source individual is unknown or unwilling to disclose their lineage, which are called Ancestry-Informative Markers (AIMs). Revealing admixture level at the locus-specific level is termed as local ancestry and can be exploited to identify signs of recent selective response and can account for genetic drift. Considering the importance of genetic admixture and local ancestry, in this mini-review, both concepts are illustrated, encompassing basics, their estimation/identification methods, tools/software used and their applications.

FishSNP: a high quality cross-species SNP database of fishes

The progress of aquaculture heavily depends on the efficient utilization of diverse genetic resources to enhance production efficiency and maximize profitability. Single nucleotide polymorphisms (SNPs) have been widely used in the study of aquaculture genomics, genetics, and breeding research since they are the most prevalent molecular markers on the genome. Currently, a large number of SNP markers from cultured fish species are scattered in individual studies, making querying complicated and data reuse problematic. We compiled relevant SNP data from literature and public databases to create a fish SNP database, FishSNP ( http://bioinfo.ihb.ac.cn/fishsnp ), and also used a unified analysis pipeline to process raw data that the author of the literature did not perform SNP calling on to obtain SNPs with high reliability. This database presently contains 45,690,243 (45 million) nonredundant SNP data for 13 fish species, with 30,288,958 (30 million) of those being high-quality SNPs. The main function of FishSNP is to search, browse, annotate and download SNPs, which provide researchers various and comprehensive associated information.

Population structure and genetic diversity in Eucalyptus pellita based on SNP markers

Eucalyptus pellita has the characteristics of rapid growth and high resistance. However, there is little research on molecular breeding of E. pellita, which is essential to shortening breeding life and selecting quality varieties. Therefore, a crucial step before selective breeding can be carried out to increase the wood quality of E. pellita is identifying genetic diversity and population structure using single nucleotide polymorphism (SNP) markers. In this study, the genetic diversity of 1st generation 196 E. pellita families from 23 geographically defined was assessed using 1,677,732 SNP markers identified by whole genome resequencing. SNP annotation showed that the ratio of non-synonymous to synonymous coding mutations was 0.83. Principal component analysis (PCA), phylogenetic tree, and population structure analysis permitted the families to be categorized into three groups, one of which (G2) contains most of the Indonesian (IDN) and Papua New Guinea (PNG) families. Genetic relationship analysis showed that IDN was closely related to PNG. Genetic diversity analysis showed that He, PIC, I, and H mean values were 0.2502, 0.2027, 0.3815, and 0.2680, respectively. PCA analysis classified various provenances in QLD into two categories (G1 and G3). The genetic diversity of G3 was higher than that of G2. The results of genetic differentiation (Fst) showed that PNG region was divided into two groups (PNG1 and PNG2), the Fst (0.172) between QLD and PNG2 region was higher than QLD and PNG1, and the Fst (0.024) between IDN and PNG1 is smaller than IDN and PNG2. A Mantel test revealed a positive correlation between the genetic and geographic distance of E. pellita. This study has a certain reference value for genetic identification, germplasm preservation, and breeding of E. pellita. Also, it provides a basis for subsequent association analysis to explore excellent alleles and introduction.

Unmasking microsatellite deceptiveness and debunking hybridization with SNPs in four marine copepod species of Calanus

Interspecific hybridization events are on the rise in natural systems due to climate change disrupting species barriers. Across taxa, microsatellites have long been the molecular markers of choice to identify admixed individuals. However, with the advent of high-throughput sequencing easing the generation of genome-wide datasets, incorrect reports of hybridization resulting from microsatellite technical artefacts have been uncovered in a growing number of taxa. In the marine zooplankton genus Calanus (Copepoda), whose species are used as climate change indicators, microsatellite markers have suggested hybridization between C. finmarchicus and C. glacialis, while other nuclear markers (InDels) never detected any admixed individuals, leaving the scientific community divided. Here, for the first time, we investigated the potential for hybridization among C. finmarchicus, C. glacialis, C. helgolandicus and C. hyperboreus using two large and independent SNP datasets. These were derived firstly from a protocol of target-capture applied to 179 individuals collected from 17 sites across the North Atlantic and Arctic Oceans, including sympatric areas, and second from published RNA sequences. All SNP-based analyses were congruent in showing that Calanus species are distinct and do not appear to hybridize. We then thoroughly re-assessed the microsatellites showing hybrids, with the support of published transcriptomes, and identified technical issues plaguing eight out of 10 microsatellites, including size homoplasy, paralogy, potential for null alleles and even two primer pairs targeting the same locus. Our study illustrates how deceptive microsatellites can be when applied to the investigation of hybridization.

Parentage exclusion of close relatives in haplodiploid species

Parentage exclusion probability is usually calculated to evaluate the informativeness of a set of markers for, and the statistical power of, a parentage analysis. Equations for parentage exclusion probability have been derived in various scenarios such as paternity exclusion when maternity is known or unknown or when candidate males are unrelated or loosely related (being from the same subpopulation) to the father. All previous work assumes a diploid species. Although marker-based parentage analyses have been conducted in haploidiploid species (such as ants, bees and wasps) for diploid offspring at the individual level or haploid offspring at the class level, rigorously derived formulations of parentage exclusion probability for haploid offspring at the individual level are lacking, which prevents the precise evaluation of the informativeness for and the statistical power of a parentage analysis. In this study we derive equations for the exclusion probability of maternity of a haploid male when multiple mother candidates (workers or queens) are unrelated or fullsibs to the mother. The usefulness of the equations is exemplified by numerical examples, and the results are discussed in the context of the study of worker reproductivity in eusocial haplodiploid species. The results are especially valuable for an optimal experimental design in determining sampling intensities (e.g. number of markers and number of individuals) to achieve satisfactory statistical power of a parentage analysis in investigating workers' reproductivity in eusocial haplodiploid species.

Genetic diversity and population structure of maize inbred lines using phenotypic traits and single nucleotide polymorphism (SNP) markers

Understanding germplasm's genetic diversity is essential for developing new and improved cultivars with stable yields under diverse environments. The objective of this study was to determine the genetic diversity and population structure of 128 maize inbred lines sourced from the International Institute of Tropical Agriculture (IITA), the International Maize and Wheat Improvement Centre (CIMMYT), and the University of KwaZulu-Natal (UKZN) using 11,450 informative single nucleotide polymorphism (SNP) markers. The inbred lines revealed highly significant (p < 0.001) levels of variability for the key phenotypic traits. The SNP markers had a mean gene diversity (GD) and polymorphic information content (PIC) of 0.40 and 0.31, respectively, indicating the existence of substantial genetic variation across the germplasm panel. The model-based population structure analysis identified three subpopulations (K = 3) among the inbred lines. This corroborated the phylogenetic analysis using phenotypic traits and molecular markers which classified the inbred lines into three groups. The findings of this study identified considerable genetic diversity for the selection of inbred lines with favourable alleles for multiple traits and could be useful to initiate marker-assisted selection (MAS) to identify significant loci associated with agronomic performance and multiple-stress tolerance.

Genome Variation Map of Domestic Qinghai-Tibet Plateau Yaks by SLAF-Seq Reveals Genetic Footprint during Artificial Selection

The yak (Bos grunniens) was domesticated in the high-altitude QTP. Research about their genetic diversity and population structure is limited. In this study, we resequenced the genome of 494 domestic yaks using Specific-Locus Amplified Fragment Sequencing (SLAF-seq). The survey was conducted on six populations sampled from isolated locations in China in order to analyze their structure and genetic diversity. These six domestic populations were clearly grouped into two independent clusters, with Jinchuan, Changtai, and Jiulong showing a tight genetic relationship with the wild yak. Nerve development pathways were enriched with GO enrichment analysis of 334 domesticated genes. Major genomic regions associated with the differentiation of domestic yaks were detected. These findings provide preliminary information on the yak genome variability, useful to understand the genomic characteristics of different populations in QTP.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.

Applying genomic approaches to delineate conservation strategies using the freshwater mussel Margaritifera margaritifera in the Iberian Peninsula as a model

Effective conservation actions to counteract the current decline of populations and species require a deep knowledge on their genetic structure. We used Single Nucleotide Polymorphisms (SNPs) to infer the population structure of the highly threatened freshwater pearl mussel Margaritifera margaritifera in the Iberian Peninsula. A total of 130 individuals were collected from 26 locations belonging to 16 basins. We obtained 31,692 SNPs through Genotyping by Sequencing (GBS) and used this dataset to infer population structure. Genetic diversity given as observed heterozygosity was low. Pairwise F_ST comparisons revealed low levels of genetic differentiation among geographically close populations. Up to 3 major genetic lineages were determined: Atlantic, Cantabrian and Douro. This structure suggests a close co-evolutionary process with brown trout (Salmo trutta), the primordial fish host of this mussel in the studied area. Some sub-basins showed some genetic structuring, whereas in others no intrapopulation differentiation was found. Our results confirm that genetic conservation units do not match individual basins, and that knowledge about the genetic structure is necessary before planning recovery plans that may involve relocation or restocking. The same reasoning should be applied to strictly freshwater species that are sessile or have restricted dispersal abilities and are currently imperiled worldwide.

Introgressive hybridisation between domestic pigs (Sus scrofa domesticus) and endemic Corsican wild boars (S. s. meridionalis): effects of human-mediated interventions

Owing to the intensified domestication process with artificial trait selection, introgressive hybridisation between domestic and wild species poses a management problem. Traditional free-range livestock husbandry, as practiced in Corsica and Sardinia, is known to facilitate hybridisation between wild boars and domestic pigs (Sus scrofa). Here, we assessed the genetic distinctness and genome-wide domestic pig ancestry levels of the Corsican wild boar subspecies S. s. meridionalis, with reference to its Sardinian conspecifics, employing a genome-wide single nucleotide polymorphism (SNP) assay and mitochondrial control region (mtCR) haplotypes. We also assessed the reliance of morphological criteria and the melanocortin-1 receptor (MC1R) coat colour gene to identify individuals with domestic introgression. While Corsican wild boars showed closest affinity to Sardinian and Italian wild boars compared to other European populations based on principal component analysis, the observation of previously undescribed mtCR haplotypes and high levels of nuclear divergence (Weir’s θ > 0.14) highlighted the genetic distinctness of Corsican S. s. meridionalis. Across three complementary analyses of mixed ancestry (i.e., STRUCTURE, PCADMIX, and ELAI), proportions of domestic pig ancestry were estimated at 9.5% in Corsican wild boars, which was significantly higher than in wild boars in Sardinia, where free-range pig keeping was banned in 2012. Comparison of morphologically pure- and hybrid-looking Corsican wild boars suggested a weak correlation between morphological criteria and genome-wide domestic pig ancestry. The study highlights the usefulness of molecular markers to assess the direct impacts of management practices on gene flow between domestic and wild species.

Genome-wide SSR markers in bottle gourd: development, characterization, utilization in assessment of genetic diversity of National Genebank of India and synteny with other related cucurbits

Lagenaria siceraria (Molina) Standley is an important cultivated crop with its immense importance in pharmaceutical industry and as vegetable. Its seed, root, stem, leaves, flower, and fruit are used as an ointment for ailment of various diseases throughout Asia. Despite its worldwide importance, informative co-dominant microsatellite markers in the bottle gourd crop are very restricted, impeding genetic improvement, cultivar identification, and phylogenetic studies. Next-generation sequencing has revolutionized the approaches for discovery, assessment, and validation of molecular markers. We conducted a genome-wide analysis, for developing SSR markers by utilizing restriction site-associated DNA sequencing (RAD-Seq) data obtained from NCBI. By performing in silico mining of microsatellite repeat motifs, we developed 45,066 perfect SSR markers. Of which 207 markers were successfully validated and 120 (57.97%) polymorphic primer pairs were utilized for an in-depth genetic diversity and population structure analysis of 96 accessions from the National Genebank of India. Tetranucleotide repeats (∼34.3%) were the most prevalent followed by trinucleotide repeats (∼30.73%), further 21.03%, 9.6%, and 4.3% of di-, penta-, and hexa-nucleotide repeats in the bottle gourd genome, respectively. Synteny of SSR markers on 11 bottle gourd linkage groups was correlated with the 7 chromosomes of cucumber (93.2%), 12 chromosomes of melon (87.4%), and 11 of watermelon (90.8%). The generated SSR markers provide a valuable tool for germplasm characterization, genetic linkage map construction, studying synteny, gene discovery, and for breeding in bottle gourd and other cucurbits species. KEY MESSAGE: Development of 45,066 perfect microsatellite markers as a valuable tool for marker assisted selection (MAS) in plant breeding.

Comparison of seven SNP calling pipelines for the next-generation sequencing data of chickens

Single nucleotide polymorphisms (SNPs) are widely used in genome-wide association studies and population genetics analyses. Next-generation sequencing (NGS) has become convenient, and many SNP-calling pipelines have been developed for human NGS data. We took advantage of a gap knowledge in selecting the appropriated SNP calling pipeline to handle with high-throughput NGS data. To fill this gap, we studied and compared seven SNP calling pipelines, which include 16GT, genome analysis toolkit (GATK), Bcftools-single (Bcftools single sample mode), Bcftools-multiple (Bcftools multiple sample mode), VarScan2-single (VarScan2 single sample mode), VarScan2-multiple (VarScan2 multiple sample mode) and Freebayes pipelines, using 96 NGS data with the different depth gradients of approximately 5X, 10X, 20X, 30X, 40X, and 50X coverage from 16 Rhode Island Red chickens. The sixteen chickens were also genotyped with a 50K SNP array, and the sensitivity and specificity of each pipeline were assessed by comparison to the results of SNP arrays. For each pipeline, except Freebayes, the number of detected SNPs increased as the input read depth increased. In comparison with other pipelines, 16GT, followed by Bcftools-multiple, obtained the most SNPs when the input coverage exceeded 10X, and Bcftools-multiple obtained the most when the input was 5X and 10X. The sensitivity and specificity of each pipeline increased with increasing input. Bcftools-multiple had the highest sensitivity numerically when the input ranged from 5X to 30X, and 16GT showed the highest sensitivity when the input was 40X and 50X. Bcftools-multiple also had the highest specificity, followed by GATK, at almost all input levels. For most calling pipelines, there were no obvious changes in SNP numbers, sensitivities or specificities beyond 20X. In conclusion, (1) if only SNPs were detected, the sequencing depth did not need to exceed 20X; (2) the Bcftools-multiple may be the best choice for detecting SNPs from chicken NGS data, but for a single sample or sequencing depth greater than 20X, 16GT was recommended. Our findings provide a reference for researchers to select suitable pipelines to obtain SNPs from the NGS data of chickens or nonhuman animals.

Automating microsatellite screening and primer design from multi-individual libraries using Micro-Primers

Analysis of intra- and inter-population diversity has become important for defining the genetic status and distribution patterns of a species and a powerful tool for conservation programs, as high levels of inbreeding could lead into whole population extinction in few generations. Microsatellites (SSR) are commonly used in population studies but discovering highly variable regions across species' genomes requires demanding computation and laboratorial optimization. In this work, we combine next generation sequencing (NGS) with automatic computing to develop a genomic-oriented tool for characterizing SSRs at the population level. Herein, we describe a new Python pipeline, named Micro-Primers, designed to identify, and design PCR primers for amplification of SSR loci from a multi-individual microsatellite library. By combining commonly used programs for data cleaning and microsatellite mining, this pipeline easily generates, from a fastq file produced by high-throughput sequencing, standard information about the selected microsatellite loci, including the number of alleles in the population subset, and the melting temperature and respective PCR product of each primer set. Additionally, potential polymorphic loci can be identified based on the allele ranges observed in the population, to easily guide the selection of optimal markers for the species. Experimental results show that Micro-Primers significantly reduces processing time in comparison to manual analysis while keeping the same quality of the results. The elapsed times at each step can be longer depending on the number of sequences to analyze and, if not assisted, the selection of polymorphic loci from multiple individuals can represent a major bottleneck in population studies.

Genotyping-in-Thousands by sequencing panel development and application to inform kokanee salmon (Oncorhynchus nerka) fisheries management at multiple scales

The ability to differentiate life history variants is vital for estimating fisheries management parameters, yet traditional survey methods can be inaccurate in mixed-stock fisheries. Such is the case for kokanee, the freshwater resident form of sockeye salmon (Oncorhynchus nerka), which exhibits various reproductive ecotypes (stream-, shore-, deep-spawning) that co-occur with each other and/or anadromous O. nerka in some systems across their pan-Pacific distribution. Here, we developed a multi-purpose Genotyping-in-Thousands by sequencing (GT-seq) panel of 288 targeted single nucleotide polymorphisms (SNPs) to enable accurate kokanee stock identification by geographic basin, migratory form, and reproductive ecotype across British Columbia, Canada. The GT-seq panel exhibited high self-assignment accuracy (93.3%) and perfect assignment of individuals not included in the baseline to their geographic basin, migratory form, and reproductive ecotype of origin. The GT-seq panel was subsequently applied to Wood Lake, a valuable mixed-stock fishery, revealing high concordance (>98%) with previous assignments to ecotype using microsatellites and TaqMan® SNP genotyping assays, while improving resolution, extending a long-term time-series, and demonstrating the scalability of this approach for this system and others.

Population analysis of retrotransposons in giraffe genomes supports RTE decline and widespread LINE1 activity in Giraffidae

BACKGROUND

The majority of structural variation in genomes is caused by insertions of transposable elements (TEs). In mammalian genomes, the main TE fraction is made up of autonomous and non-autonomous non-LTR retrotransposons commonly known as LINEs and SINEs (Long and Short Interspersed Nuclear Elements). Here we present one of the first population-level analysis of TE insertions in a non-model organism, the giraffe. Giraffes are ruminant artiodactyls, one of the few mammalian groups with genomes that are colonized by putatively active LINEs of two different clades of non-LTR retrotransposons, namely the LINE1 and RTE/BovB LINEs as well as their associated SINEs. We analyzed TE insertions of both types, and their associated SINEs in three giraffe genome assemblies, as well as across a population level sampling of 48 individuals covering all extant giraffe species.

RESULTS

The comparative genome screen identified 139,525 recent LINE1 and RTE insertions in the sampled giraffe population. The analysis revealed a drastically reduced RTE activity in giraffes, whereas LINE1 is still actively propagating in the genomes of extant (sub)-species. In concert with the extremely low activity of the giraffe RTE, we also found that RTE-dependent SINEs, namely Bov-tA and Bov-A2, have been virtually immobile in the last 2 million years. Despite the high current activity of the giraffe LINE1, we did not find evidence for the presence of currently active LINE1-dependent SINEs. TE insertion heterozygosity rates differ among the different (sub)-species, likely due to divergent population histories.

CONCLUSIONS

The horizontally transferred RTE/BovB and its derived SINEs appear to be close to inactivation and subsequent extinction in the genomes of extant giraffe species. This is the first time that the decline of a TE family has been meticulously analyzed from a population genetics perspective. Our study shows how detailed information about past and present TE activity can be obtained by analyzing large-scale population-level genomic data sets.

High levels of connectivity over large distances in the diadematid sea urchin Centrostephanus sylviae

Most benthic marine invertebrates with sedentary benthic adult phases have planktonic larvae that permit connectivity between geographically isolated populations. Planktonic larval duration and oceanographic processes are vital to connecting populations of species inhabiting remote and distant islands. In the present study, we analyzed the population genetic structure of the sea urchin Centrostephanus sylviae, which inhabits only the Juan Fernández Archipelago and the Desventuradas islands, separated by more than 800 km. For 92 individuals collected from Robinson Crusoe and Selkirk Islands (Juan Fernández Archipelago) and San Ambrosio Island (Desventuradas Islands), 7,067 single nucleotide polymorphisms (SNPs) were obtained. The results did not show a spatial genetic structure for C. sylviae; relative high migration rates were revealed between the islands. An analysis of the water circulation pattern in the area described a predominant northward water flow with periods of inverted flow, suggesting that larvae could move in both directions. Overall, this evidence suggests that C. sylviae comprises a single large population composed of individuals separated by more than 800 km.

Characterization of the genetic basis of local adaptation of wheat landraces from Iran and Pakistan using genome-wide association study

Characterization of genomic regions underlying adaptation of landraces can reveal a quantitative genetics framework for local wheat (Triticum aestivum L.) adaptability. A collection of 512 wheat landraces from the eastern edge of the Fertile Crescent in Iran and Pakistan were genotyped using genome-wide single nucleotide polymorphism markers generated by genotyping-by-sequencing. The minor allele frequency (MAF) and the heterozygosity (H) of Pakistani wheat landraces (MAF = 0.19, H = 0.008) were slightly higher than the Iranian wheat landraces (MAF = 0.17, H = 0.005), indicating that Pakistani landraces were slightly more genetically diverse. Population structure analysis clearly separated the Pakistani landraces from Iranian landraces, which indicates two separate adaptability trajectories. The large-scale agro-climatic data of seven variables were quite dissimilar between Iran and Pakistan as revealed by the correlation coefficients. Genome-wide association study identified 91 and 58 loci using agroclimatic data, which likely underpin local adaptability of the wheat landraces from Iran and Pakistan, respectively. Selective sweep analysis identified significant hits on chromosomes 4A, 4B, 6B, 7B, 2D, and 6D, which were colocalized with the loci associated with local adaptability and with some known genes related to flowering time and grain size. This study provides insight into the genetic diversity with emphasis on the genetic architecture of loci involved in adaptation to local environments, which has breeding implications.

Waste not, want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single-nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black-capped vireo source-sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.

A Semi-Automated SNP-Based Approach for Contaminant Identification in Biparental Polyploid Populations of Tropical Forage Grasses

Artificial hybridization plays a fundamental role in plant breeding programs since it generates new genotypic combinations that can result in desirable phenotypes. Depending on the species and mode of reproduction, controlled crosses may be challenging, and contaminating individuals can be introduced accidentally. In this context, the identification of such contaminants is important to avoid compromising further selection cycles, as well as genetic and genomic studies. The main objective of this work was to propose an automated multivariate methodology for the detection and classification of putative contaminants, including apomictic clones (ACs), self-fertilized individuals, half-siblings (HSs), and full contaminants (FCs), in biparental polyploid progenies of tropical forage grasses. We established a pipeline to identify contaminants in genotyping-by-sequencing (GBS) data encoded as allele dosages of single nucleotide polymorphism (SNP) markers by integrating principal component analysis (PCA), genotypic analysis (GA) measures based on Mendelian segregation, and clustering analysis (CA). The combination of these methods allowed for the correct identification of all contaminants in all simulated progenies and the detection of putative contaminants in three real progenies of tropical forage grasses, providing an easy and promising methodology for the identification of contaminants in biparental progenies of tetraploid and hexaploid species. The proposed pipeline was made available through the polyCID Shiny app and can be easily coupled with traditional genetic approaches, such as linkage map construction, thereby increasing the efficiency of breeding programs.

Conservation implications of genetic structure in the narrowest endemic quillwort from the Eastern Amazon

The quillwort Isoëtes cangae is a critically endangered species occurring in a single lake in Serra dos Carajás, Eastern Amazon. Low genetic diversity and small effective population sizes (N e) are expected for narrow endemic species (NES). Conservation biology studies centered in a single species show some limitations, but they are still useful considering the limited time and resources available for protection of species at risk of extinction. Here, we evaluated the genetic diversity, population structure, N e, and minimum viable population (MVP) of I. cangae to provide information for effective conservation programs. Our analyses were based on 55 individuals collected from the Amendoim Lake and 35,638 neutral SNPs. Our results indicated a single panmictic population, moderate levels of genetic diversity, and N e in the order of thousands, contrasting the expected for NES. Negative FIS values were also found, suggesting that I. cangae is not under risk of inbreeding depression. Our findings imply that I. cangae contains enough genetic diversity to ensure evolutionary potential and that all individuals should be treated as one demographic unit. These results provide essential information to optimize ex situ conservation efforts and genetic diversity monitoring, which are currently applied to guide I. cangae conservation plans.

A SNP assay for assessing diversity in immune genes in the honey bee (Apis mellifera L.)

With a growing number of parasites and pathogens experiencing large-scale range expansions, monitoring diversity in immune genes of host populations has never been so important because it can inform on the adaptive potential to resist the invaders. Population surveys of immune genes are becoming common in many organisms, yet they are missing in the honey bee (Apis mellifera L.), a key managed pollinator species that has been severely affected by biological invasions. To fill the gap, here we identified single nucleotide polymorphisms (SNPs) in a wide range of honey bee immune genes and developed a medium-density assay targeting a subset of these genes. Using a discovery panel of 123 whole-genomes, representing seven A. mellifera subspecies and three evolutionary lineages, 180 immune genes were scanned for SNPs in exons, introns (< 4 bp from exons), 3' and 5´UTR, and < 1 kb upstream of the transcription start site. After application of multiple filtering criteria and validation, the final medium-density assay combines 91 quality-proved functional SNPs marking 89 innate immune genes and these can be readily typed using the high-sample-throughput iPLEX MassARRAY system. This medium-density-SNP assay was applied to 156 samples from four countries and the admixture analysis clustered the samples according to their lineage and subspecies, suggesting that honey bee ancestry can be delineated from functional variation. In addition to allowing analysis of immunogenetic variation, this newly-developed SNP assay can be used for inferring genetic structure and admixture in the honey bee.

Genome-wide SNP analysis to assess the genetic population structure and diversity of Acrocomia species

Acrocomia (Arecaceae) is a genus widely distributed in tropical and subtropical America that has been achieving economic interest due to the great potential of oil production of some of its species. In particular A. aculeata, due to its vocation to supply oil with the same productive capacity as the oil palm (Elaeis guineenses) even in areas with water deficit. Although eight species are recognized in the genus, the taxonomic classification based on morphology and geographic distribution is still controversial. Knowledge about the genetic diversity and population structure of the species is limited, which has limited the understanding of the genetic relationships and the orientation of management, conservation, and genetic improvement activities of species of the genus. In the present study, we analyzed the genomic diversity and population structure of Acrocomia genus, including 172 samples from seven species, with a focus on A. aculeata with 117 samples covering a wide geographical area of occurrence of the species, using Single Nucleotide Polymorphism (SNP) markers originated from Genotyping By Sequencing (GBS).The genetic structure of the Acrocomia species were partially congruent with the current taxonomic classification based on morphological characters, recovering the separation of the species A. aculeata, A. totai, A. crispa and A. intumescens as distinct taxonomic groups. However, the species A. media was attributed to the cluster of A. aculeata while A. hassleri and A. glauscescens were grouped together with A. totai. The species that showed the highest and lowest genetic diversity were A. totai and A. media, respectively. When analyzed separately, the species A. aculeata showed a strong genetic structure, forming two genetic groups, the first represented mainly by genotypes from Brazil and the second by accessions from Central and North American countries. Greater genetic diversity was found in Brazil when compared to the other countries. Our results on the genetic diversity of the genus are unprecedented, as is also establishes new insights on the genomic relationships between Acrocomia species. It is also the first study to provide a more global view of the genomic diversity of A. aculeata. We also highlight the applicability of genomic data as a reference for future studies on genetic diversity, taxonomy, evolution and phylogeny of the Acrocomia genus, as well as to support strategies for the conservation, exploration and breeding of Acrocomia species and in particular A. aculeata.

Genetic hypervariability of a Northeastern Atlantic venomous rockfish

BACKGROUND

Understanding the interplay between climate and current and historical factors shaping genetic diversity is pivotal to infer changes in marine species range and communities' composition. A phylogeographical break between the Atlantic and the Mediterranean has been documented for several marine organisms, translating into limited dispersal between the two basins.

METHODS

In this study, we screened the intraspecific diversity of 150 individuals of the Madeira rockfish (Scorpaena maderensis) across its distributional range (seven sampling locations in the Atlantic and Mediterranean basins) using the mitochondrial control region and the nuclear S7 first intron.

RESULTS

The present work is the most comprehensive study done for this species, yielding no genetic structure across sampled locations and no detectable Atlantic-Mediterranean break in connectivity. Our results reveal deep and hyper-diverse bush-like genealogies with large numbers of singletons and very few shared haplotypes. The genetic hyper-diversity found for the Madeira rockfish is relatively uncommon in rocky coastal species, whose dispersal capability is limited by local oceanographic patterns. The effect of climate warming on the distribution of the species is discussed.

Whole genome sequencing reveals high differentiation, low levels of genetic diversity and short runs of homozygosity among Swedish wels catfish

The use of genetic markers in the context of conservation is largely being outcompeted by whole-genome data. Comparative studies between the two are sparse, and the knowledge about potential effects of this methodology shift is limited. Here, we used whole-genome sequencing data to assess the genetic status of peripheral populations of the wels catfish (Silurus glanis), and discuss the results in light of a recent microsatellite study of the same populations. The Swedish populations of the wels catfish have suffered from severe declines during the last centuries and persists in only a few isolated water systems. Fragmented populations generally are at greater risk of extinction, for example due to loss of genetic diversity, and may thus require conservation actions. We sequenced individuals from the three remaining native populations (Båven, Emån, and Möckeln) and one reintroduced population of admixed origin (Helge å), and found that genetic diversity was highest in Emån but low overall, with strong differentiation among the populations. No signature of recent inbreeding was found, but a considerable number of short runs of homozygosity were present in all populations, likely linked to historically small population sizes and bottleneck events. Genetic substructure within any of the native populations was at best weak. Individuals from the admixed population Helge å shared most genetic ancestry with the Båven population (72%). Our results are largely in agreement with the microsatellite study, and stresses the need to protect these isolated populations at the northern edge of the distribution of the species.

Mytilus trossulus introgression and consequences for shell traits in longline cultivated mussels

Mussels belonging to the Mytilus species complex (M. edulis, ME; M. galloprovincialis, MG; and M. trossulus, MT) often occur in sympatry, facilitating introgressive hybridization. This may be further promoted by mussel aquaculture practices, with MT introgression often resulting in commercially unfavourable traits such as low meat yield and weak shells. To investigate the relationship between genotype and shell phenotype, genetic and morphological variability was quantified across depth (1 m to 7 m) along a cultivation rope at a mussel farm on the West coast of Scotland. A single nuclear marker (Me15/16) and a novel panel of 33 MT-diagnostic single nucleotide polymorphisms were used to evaluate stock structure and the extent of MT introgression across depth. Variation in shell strength, determined as the maximum compression force for shell puncture, and shell shape using geometric morphometric analysis were evaluated in relation to cultivation depth and the genetic profiles of the mussels. Overall, ME was the dominant genotype across depth, followed by ME × MG hybrids and smaller quantities of ME × MT hybrids and pure MT individuals. In parallel, we identified multiple individuals that were either predominantly homozygous or heterozygous for MT-diagnostic alleles, likely representing pure MT and first-generation ME × MT hybrids, respectively. Both the proportion of individuals carrying MT alleles and MT allele frequency declined with depth. Furthermore, MT-introgressed individuals had significantly weaker and more elongate shells than nonintrogressed individuals. This study provides detailed insights into stock structure along a cultivation rope and suggests that practical methods to assess shell strength and shape of cultivated mussels may facilitate the rapid identification of MT, limiting the impact of this commercially damaging species.

Genetic association with boldness and maternal performance in a free-ranging population of grey seals (Halichoerus grypus)

Individual variation in quantitative traits clearly influence many ecological and evolutionary processes. Moderate to high heritability estimates of personality and life-history traits suggest some level of genetic control over these traits. Yet, we know very little of the underlying genetic architecture of phenotypic variation in the wild. In this study, we used a candidate gene approach to investigate the association of genetic variants with repeated measures of boldness and maternal performance traits (weaning mass and lactation duration) collected over an 11- and 28-year period, respectively, in a free-ranging population of grey seals on Sable Island National Park Reserve, Canada. We isolated and re-sequenced five genes: dopamine receptor D4 (DRD4), serotonin transporter (SERT), oxytocin receptor (OXTR), and melanocortin receptors 1 (MC1R) and 5 (MC5R). We discovered single nucleotide polymorphisms (SNPs) in each gene; and, after accounting for loci in linkage disequilibrium and filtering due to missing data, we were able to test for genotype-phenotype relationships at seven loci in three genes (DRD4, SERT, and MC1R). We tested for association between these loci and traits of 180 females having extreme shy-bold phenotypes using mixed-effects models. One locus within SERT was significantly associated with boldness (effect size = 0.189) and a second locus within DRD4 with weaning mass (effect size = 0.232). Altogether, genotypes explained 6.52-13.66% of total trait variation. Our study substantiates SERT and DRD4 as important determinants of behaviour, and provides unique insight into the molecular mechanisms underlying maternal performance variation in a marine predator.

Population Genetics of the Invasive Red Fox, Vulpes vulpes, in South-Eastern Australia

The use of genetic information in conservation biology has become more widespread with genetic information more readily available for non-model organisms. It has also been recognized that genetic information from invasive species can inform their management and control. The red fox poses a significant threat to Australian native fauna and the agricultural industry. Despite this, there are few recently published studies investigating the population genetics of foxes in Australia. This study investigated the population genetics of 94 foxes across the Illawarra and Shoalhaven regions of New South Wales, Australia. Diversity Array sequencing technology was used to genotype a large number of single nucleotide polymorphisms (N = 33,375). Moderate genetic diversity and relatedness were observed across the foxes sampled. Low to moderate levels of inbreeding, high-levels of identity-by-state values, as well as high identity-by-descent values were also found. There was limited evidence for population genetic structure among the foxes across the landscape sampled, supporting the presence of a single population across the study area. This indicates that there may be no barriers hindering fox dispersal across the landscape.

Genome-wide association analysis of bean fly resistance and agro-morphological traits in common bean

The bean fly (Ophiomyia spp) is a key insect pest causing significant crop damage and yield loss in common bean (Phaseolus vulgaris L., 2n = 2x = 22). Development and deployment of agronomic superior and bean fly resistant common bean varieties aredependent on genetic variation and the identification of genes and genomic regions controlling economic traits. This study's objective was to determine the population structure of a diverse panel of common bean genotypes and deduce associations between bean fly resistance and agronomic traits based on single nucleotide polymorphism (SNP) markers. Ninety-nine common bean genotypes were phenotyped in two seasons at two locations and genotyped with 16 565 SNP markers. The genotypes exhibited significant variation for bean fly damage severity (BDS), plant mortality rate (PMR), and pupa count (PC). Likewise, the genotypes showed significant variation for agro-morphological traits such as days to flowering (DTF), days to maturity (DTM), number of pods per plant (NPP), number of seeds per pod (NSP), and grain yield (GYD). The genotypes were delineated into two populations, which were based on the Andean and Mesoamerican gene pools. The genotypes exhibited a minimum membership coefficient of 0.60 to their respective populations. Eighty-three significant (P<0.01) markers were identified with an average linkage disequilibrium of 0.20 at 12Mb across the 11 chromosomes. Three markers were identified, each having pleiotropic effects on two traits: M100049197 (BDS and NPP), M3379537 (DTF and PC), and M13122571 (NPP and GYD). The identified markers are useful for marker-assisted selection in the breeding program to develop common bean genotypes with resistance to bean fly damage.

Evolutionary genomics of endangered Hawaiian tree snails (Achatinellidae: Achatinellinae) for conservation of adaptive capacity

Phylogenomic studies can provide insights into speciation, adaptation, and extinction, while providing a roadmap for conservation. Hawaiian tree snails are a model system for an adaptive radiation facing an extinction crisis. In the last 5 years, nearly all populations of Hawaiian tree snails across the 30 remaining species in the subfamily Achatinellinae (Achatinellidae) have declined from hundreds or thousands in the wild down to undetectable levels. Nearly 100 species historically occurred across dramatic environmental gradients on five of the Hawaiian Islands, but habitat loss, overcollection, and predation by invasive species have decimated populations. As such, this system offers the opportunity to integrate efforts to conserve evolutionary potential into conservation planning for a rapidly declining subfamily. Here, we used genome-wide, restriction-site associated DNA sequencing (RADseq), along with mitochondrial genome reconstruction, to resolve evolutionary relationships to inform conservation efforts. Phylogenetic analysis of nearly 400k genome-wide SNPs from 59 populations and 25 species across six genera in the family Achatinellidae, was generally concordant with taxonomy, geography, and mtDNA with several notable exceptions; mtDNA was unable to resolve some deeper nodes (e.g., the monophyly of Achatinella), while SNP data did not resolve as many shallow nodes. Both phylogenetic and coalescent analysis revealed deep divergences between populations within Achatinella mustelina that were consistent with species-level differences. Given cryptic species-level divergence within populations that are geographically proximate, they are at higher risk of extirpation from invasive predators and climate change than previously assumed. This study clarifies evolutionary relationships within this model system for adaptive radiation, forming the basis for conservation strategies such as translocation, captive rearing, and hybridization trials to prevent the loss of capacity to adapt to rapidly changing environmental conditions.

Sample identification and pedigree reconstruction in Wolverine (Gulo gulo) using SNP genotyping of non-invasive samples

For conservation genetic studies using non-invasively collected samples, genome-wide data may be hard to acquire. Until now, such studies have instead mostly relied on analyses of traditional genetic markers such as microsatellites (SSRs). Recently, high throughput genotyping of single nucleotide polymorphisms (SNPs) has become available, expanding the use of genomic methods to include non-model species of conservation concern. We have developed a 96-marker SNP array for use in applied conservation monitoring of the Scandinavian wolverine (Gulo gulo) population. By genotyping more than a thousand non-invasively collected samples, we were able to obtain precise estimates of different types of genotyping errors and sample dropout rates. The SNP panel significantly outperforms the SSR markers (and DBY intron markers for sexing) both in terms of precision in genotyping, sex assignment and individual identification, as well as in the proportion of samples successfully genotyped. Furthermore, SNP genotyping offers a simplified laboratory and analysis pipeline with fewer samples needed to be repeatedly genotyped in order to obtain reliable consensus data. In addition, we utilised a unique opportunity to successfully demonstrate the application of SNP genotype data for reconstructing pedigrees in wild populations, by validating the method with samples from wild individuals with known relatedness. By offering a simplified workflow with improved performance, we anticipate this methodology will facilitate the use of non-invasive samples to improve genetic management of many different types of populations that have previously been challenging to survey.

Population Genetic Analyses of Botrytis cinerea Isolates From Michigan Vineyards Using a High-Throughput Marker System Approach

As sequencing costs continue to decrease, new tools are being developed for assessing pathogen diversity and population structure. Traditional marker types, such as microsatellites, are often more cost effective than single-nucleotide polymorphism (SNP) panels when working with small numbers of individuals, but may not allow for fine scale evaluation of low or moderate structure in populations. Botrytis cinerea is a necrotrophic plant pathogen with high genetic variability that can infect more than 200 plant species worldwide. A panel of 52 amplicons were sequenced for 82 isolates collected from four Michigan vineyards representing 2 years of collection and varying fungicide resistance. A panel of nine microsatellite markers previously described was also tested across 74 isolates from the same population. A microsatellite and SNP marker analysis of B. cinerea populations was performed to assess the genetic diversity and population structure of Michigan vineyards, and the results from both marker types were compared. Both methods were able to detect population structure associated with resistance to the individual fungicides thiabendazole and boscalid, and multiple fungicide resistance (MFR). Microsatellites were also able to differentiate population structure associated with another fungicide, fluopyram, while SNPs were able to additionally differentiate structure based on year. For both methods, AMOVA results were similar, with microsatellite results explaining a smaller portion of the variation compared with the SNP results. The SNP-based markers presented here were able to successfully differentiate population structure similar to microsatellite results. These SNP markers represent new tools to discriminate B. cinerea isolates within closely related populations using multiple targeted sequences.

Soybean anthracnose caused by Colletotrichum species: Current status and future prospects

Soybean (Glycine max) is one of the most important cultivated plants worldwide as a source of protein-rich foods and animal feeds. Anthracnose, caused by different lineages of the hemibiotrophic fungus Colletotrichum, is one of the main limiting factors to soybean production. Losses due to anthracnose have been neglected, but their impact may threaten up to 50% of the grain production.

TAXONOMY

While C. truncatum is considered the main species associated with soybean anthracnose, recently other species have been reported as pathogenic on this host. Until now, it has not been clear whether the association of new Colletotrichum species with the disease is related to emerging species or whether it is due to the undergoing changes in the taxonomy of the genus.

DISEASE SYMPTOMS

Typical anthracnose symptoms are pre- and postemergence damping-off; dark, depressed, and irregular spots on cotyledons, stems, petioles, and pods; and necrotic laminar veins on leaves that can result in premature defoliation. Symptoms may evolve to pod rot, immature opening of pods, and premature germination of grains.

CHALLENGES

As accurate species identification of the causal agent is decisive for disease control and prevention, in this work we review the taxonomic designation of Colletotrichum isolated from soybean to understand which lineages are pathogenic on this host. We also present a comprehensive literature review of soybean anthracnose, focusing on distribution, symptomatology, epidemiology, disease management, identification, and diagnosis. We consider the knowledge emerging from population studies and comparative genomics of Colletotrichum spp. associated with soybean providing future perspectives in the identification of molecular factors involved in the pathogenicity process.

USEFUL WEBSITE

Updates on Colletotrichum can be found at http://www.colletotrichum.org/. All available Colletotrichum genomes on GenBank can be viewed at http://www.colletotrichum.org/genomics/.

Spontaneous hybridization and introgression between walleye (Sander vitreus) and sauger (Sander canadensis) in two large reservoirs: Insights from genotyping by sequencing

Anthropogenic activities may facilitate undesirable hybridization and genomic introgression between fish species. Walleye (Sander vitreus) and sauger (Sander canadensis) are economically valuable freshwater species that can spontaneously hybridize in areas of sympatry. Levels of genomic introgression between walleye and sauger may be increased by modifications to waterbodies (e.g., reservoir development) and inadvertent propagation of hybrids in stocking programs. We used genotyping by sequencing (GBS) to examine 217 fish from two large reservoirs with mixed populations of walleye and sauger in Saskatchewan, Canada (Lake Diefenbaker, Tobin Lake). Analyses with 20,038 (r90) and 478 (r100) single nucleotide polymorphisms clearly resolved walleye and sauger, and classified hybrids with high confidence. F1, F2, and multigeneration hybrids were detected in Lake Diefenbaker, indicating potentially high levels of genomic introgression. In contrast, only F1 hybrids were detected in Tobin Lake. Field classification of fish was unreliable; 7% of fish were misidentified based on broad species categories. Important for activities such as brood stock selection, 12 of 173 (7%) fish field identified as pure walleye, and one of 24 (4%) identified as pure sauger were actually hybrids. In addition, two of 15 (13%) field-identified hybrids were actually pure walleye or sauger. We conclude that hybridization and introgression are occurring in Saskatchewan reservoirs and that caution is warranted when using these populations in stocking programs. GBS offers a powerful and flexible tool for examining hybridization without preidentification of informative loci, eliminating some of the key challenges associated with other marker types.

Primate landscape genetics: A review and practical guide

Landscape genetics is an emerging field that integrates population genetics, landscape ecology, and spatial statistics to investigate how geographical and environmental features and evolutionary processes such as gene flow, genetic drift, and selection structure genetic variation at both the population and individual levels, with implications for ecology, evolution, and conservation biology. Despite being particularly well suited for primatologists, this method is currently underutilized. Here, we synthesize the current state of research on landscape genetics in primates. We begin by outlining how landscape genetics has been used to disentangle the drivers of diversity, followed by a review of how landscape genetic methods have been applied to primates. This is followed by a section highlighting special considerations when applying the methods to primates, and a practical guide to facilitate further landscape genetics studies using both existing and de novo datasets. We conclude by exploring future avenues of inquiry that could be facilitated by recent developments as well as underdeveloped applications of landscape genetics to primates.

Thermal adaptation rather than demographic history drives genetic structure inferred by copy number variants in a marine fish

Increasing evidence shows that structural variants represent an overlooked aspect of genetic variation with consequential evolutionary roles. Among those, copy number variants (CNVs), including duplicated genomic regions and transposable elements (TEs), may contribute to local adaptation and/or reproductive isolation among divergent populations. Those mechanisms suppose that CNVs could be used to infer neutral and/or adaptive population genetic structure, whose study has been restricted to microsatellites, mitochondrial DNA and Amplified fragment length polymorphism markers in the past and more recently the use of single nucleotide polymorphisms (SNPs). Taking advantage of recent developments allowing CNV analysis from RAD-seq data, we investigated how variation in fitness-related traits, local environmental conditions and demographic history are associated with CNVs, and how subsequent copy number variation drives population genetic structure in a marine fish, the capelin (Mallotus villosus). We collected 1538 DNA samples from 35 sampling sites in the north Atlantic Ocean and identified 6620 putative CNVs. We found associations between CNVs and the gonadosomatic index, suggesting that six duplicated regions could affect female fitness by modulating oocyte production. We also detected 105 CNV candidates associated with water temperature, among which 20% corresponded to genomic regions located within the sequence of protein-coding genes, suggesting local adaptation to cold water by means of gene sequence amplification. We also identified 175 CNVs associated with the divergence of three previously defined parapatric glacial lineages, of which 24% were located within protein-coding genes, making those loci potential candidates for reproductive isolation. Lastly, our analyses unveiled a hierarchical, complex CNV population structure determined by temperature and local geography, which was in stark contrast to that inferred based on SNPs in a previous study. Our findings underline the complementarity of those two types of genomic variation in population genomics studies.