Genotyping‐by‐sequencing for estimating relatedness in nonmodel organisms: Avoiding the trap of precise bias

Spatial population genetic structure of Caquetaia kraussii (Steindachner, 1878) evidenced by species-specific microsatellite loci in the middle and low basin of the Cauca River, Colombia

The adaptative responses and divergent evolution shown in the environments habited by the Cichlidae family allow to understand different biological properties, including fish genetic diversity and structure studies. In a zone that has been historically submitted to different anthropogenic pressures, this study assessed the genetic diversity and population structure of cichlid Caquetaia kraussii, a sedentary species with parental care that has a significant ecological role for its contribution to redistribution and maintenance of sedimentologic processes in its distribution area. This study developed de novo 16 highly polymorphic species-specific microsatellite loci that allowed the estimation of the genetic diversity and differentiation in 319 individuals from natural populations in the area influenced by the Ituango hydroelectric project in the Colombian Cauca River. Caquetaia kraussii exhibits high genetic diversity levels (Ho: 0.562-0.885; He: 0.583-0.884) in relation to the average neotropical cichlids and a three group-spatial structure: two natural groups upstream and downstream the Nechí River mouth, and one group of individuals with high relatedness degree, possibly independently formed by founder effect in the dam zone. The three genetic groups show recent bottlenecks, but only the two natural groups have effective population size that suggest their long-term permanence. The information generated is relevant not only for management programs and species conservation purposes, but also for broadening the available knowledge on the factors influencing neotropical cichlids population genetics.

Species-specific SNP arrays for non-invasive genetic monitoring of a vulnerable bat

Genetic tagging from scats is one of the minimally invasive sampling (MIS) monitoring approaches commonly used to guide management decisions and evaluate conservation efforts. Microsatellite markers have traditionally been used but are prone to genotyping errors. Here, we present a novel method for individual identification in the Threatened ghost bat Macroderma gigas using custom-designed Single Nucleotide Polymorphism (SNP) arrays on the MassARRAY system. We identified 611 informative SNPs from DArTseq data from which three SNP panels (44-50 SNPs per panel) were designed. We applied SNP genotyping and molecular sexing to 209 M. gigas scats collected from seven caves in the Pilbara, Western Australia, employing a two-step genotyping protocol and identifying unique genotypes using a custom-made R package, ScatMatch. Following data cleaning, the average amplification rate was 0.90 ± 0.01 and SNP genotyping errors were low (allelic dropout 0.003 ± 0.000) allowing clustering of scats based on one or fewer allelic mismatches. We identified 19 unique bats (9 confirmed/likely males and 10 confirmed/likely females) from a maternity and multiple transitory roosts, with two male bats detected using roosts, 9 km and 47 m apart. The accuracy of our SNP panels enabled a high level of confidence in the identification of individual bats. Targeted SNP genotyping is a valuable tool for monitoring and tracking of non-model species through a minimally invasive sampling approach.

Social antagonism facilitates supergene expansion in ants

Antagonistic selection has long been considered a major driver of the formation and expansion of sex chromosomes. For example, sexually antagonistic variation on an autosome can select for suppressed recombination between that autosome and the sex chromosome, leading to a neo-sex chromosome. Autosomal supergenes, chromosomal regions containing tightly linked variants affecting the same complex trait, share similarities with sex chromosomes, raising the possibility that sex chromosome evolution models can explain the evolution of genome structure and recombination in other contexts. We tested this premise in a Formica ant species, wherein we identified four supergene haplotypes on chromosome 3 underlying colony social organization and sex ratio. We discovered a novel rearranged supergene variant (9r) on chromosome 9 underlying queen miniaturization. The 9r is in strong linkage disequilibrium with one chromosome 3 haplotype (P2) found in multi-queen (polygyne) colonies. We suggest that queen miniaturization is strongly disfavored in the single-queen (monogyne) background and is thus socially antagonistic. As such, divergent selection experienced by ants living in alternative social "environments" (monogyne and polygyne) may have contributed to the emergence of a genetic polymorphism on chromosome 9 and associated queen-size dimorphism. Consequently, an ancestral polygyne-associated haplotype may have expanded to include the polymorphism on chromosome 9, resulting in a larger region of suppressed recombination spanning two chromosomes. This process is analogous to the formation of neo-sex chromosomes and consistent with models of expanding regions of suppressed recombination. We propose that miniaturized queens, 16%-20% smaller than queens without 9r, could be incipient intraspecific social parasites.

Best practices for genotype imputation from low-coverage sequencing data in natural populations

Monitoring genetic diversity in wild populations is a central goal of ecological and evolutionary genetics and is critical for conservation biology. However, genetic studies of nonmodel organisms generally lack access to species-specific genotyping methods (e.g. array-based genotyping) and must instead use sequencing-based approaches. Although costs are decreasing, high-coverage whole-genome sequencing (WGS), which produces the highest confidence genotypes, remains expensive. More economical reduced representation sequencing approaches fail to capture much of the genome, which can hinder downstream inference. Low-coverage WGS combined with imputation using a high-confidence reference panel is a cost-effective alternative, but the accuracy of genotyping using low-coverage WGS and imputation in nonmodel populations is still largely uncharacterized. Here, we empirically tested the accuracy of low-coverage sequencing (0.1-10×) and imputation in two natural populations, one with a large (n = 741) reference panel, rhesus macaques (Macaca mulatta), and one with a smaller (n = 68) reference panel, gelada monkeys (Theropithecus gelada). Using samples sequenced to coverage as low as 0.5×, we could impute genotypes at >95% of the sites in the reference panel with high accuracy (median r2  ≥ 0.92). We show that low-coverage imputed genotypes can reliably calculate genetic relatedness and population structure. Based on these data, we also provide best practices and recommendations for researchers who wish to deploy this approach in other populations, with all code available on GitHub (https://github.com/mwatowich/LoCSI-for-non-model-species). Our results endorse accurate and effective genotype imputation from low-coverage sequencing, enabling the cost-effective generation of population-scale genetic datasets necessary for tackling many pressing challenges of wildlife conservation.

2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics

Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.

Sex ratio and relatedness in the Griffon vulture (Gyps fulvus) population of Serbia

Background

Once a widespread species across the region of Southeast Europe, the Griffon vulture is now confined to small and isolated populations across the Balkan Peninsula. The population from Serbia represents its biggest and most viable population that can serve as an important reservoir of genetic diversity from which the birds can be used for the region's reintroduction programmes. The available genetic data for this valuable population are scarce and as a protected species that belongs to the highly endangered vulture group, it needs to be well described so that it can be properly managed and used as a restocking population. Considering the serious recent bottleneck event that the Griffon vulture population from Serbia experienced we estimated the overall relatedness among the birds from this population. Sex ratio, another important parameter that shows the vitality and strength of the population was evaluated as well.

Methods

During the annual monitoring that was performed in the period from 2013-2021, we collected blood samples from individual birds that were marked in the nests. In total, 169 samples were collected and each was used for molecular sexing while 58 presumably unrelated birds from different nests were used for inbreeding and relatedness analyses. The relatedness was estimated using both biparentally (10 microsatellite loci) and uniparentally (Cytb and D-loop I of mitochondrial DNA) inherited markers.

Results

The level of inbreeding was relatively high and on average it was 8.3% while the mean number of relatives for each bird was close to three. The sex ratio was close to 1:1 and for the analysed period of 9 years, it didn't demonstrate a statistically significant deviation from the expected ratio of 1:1, suggesting that this is a stable and healthy population. Our data suggest that, even though a relatively high level of inbreeding can be detected among the individual birds, the Griffon vulture population from Serbia can be used as a source population for restocking and reintroduction programmes in the region. These data combined with previously observed genetic differentiation between the populations from the Iberian and Balkan Peninsulas suggest that the introduction of foreign birds should be avoided and that local birds should be used instead.

Seascape genomics of common dolphins (Delphinus delphis) reveals adaptive diversity linked to regional and local oceanography

Background

High levels of standing genomic variation in wide-ranging marine species may enhance prospects for their long-term persistence. Patterns of connectivity and adaptation in such species are often thought to be influenced by spatial factors, environmental heterogeneity, and oceanographic and geomorphological features. Population-level studies that analytically integrate genome-wide data with environmental information (i.e., seascape genomics) have the potential to inform the spatial distribution of adaptive diversity in wide-ranging marine species, such as many marine mammals. We assessed genotype-environment associations (GEAs) in 214 common dolphins (Delphinus delphis) along > 3000 km of the southern coast of Australia.

Results

We identified 747 candidate adaptive SNPs out of a filtered panel of 17,327 SNPs, and five putatively locally-adapted populations with high levels of standing genomic variation were disclosed along environmentally heterogeneous coasts. Current velocity, sea surface temperature, salinity, and primary productivity were the key environmental variables associated with genomic variation. These environmental variables are in turn related to three main oceanographic phenomena that are likely affecting the dispersal of common dolphins: (1) regional oceanographic circulation, (2) localised and seasonal upwellings, and (3) seasonal on-shelf circulation in protected coastal habitats. Signals of selection at exonic gene regions suggest that adaptive divergence is related to important metabolic traits.

Conclusion

To the best of our knowledge, this represents the first seascape genomics study for common dolphins (genus Delphinus). Information from the associations between populations and their environment can assist population management in forecasting the adaptive capacity of common dolphins to climate change and other anthropogenic impacts.

Reduced genetic diversity associated with the northern expansion of an amphibian species with high habitat specialization, Ascaphus truei , resolved using two types of genetic markers

Reconstruction of historical relationships between geographic regions within a species’ range can indicate dispersal patterns and help predict future responses to shifts in climate. Ascaphus truei (coastal tailed frog) is an indicator species of the health of forests and perennial streams in the Coastal and Cascade Mountains of the Pacific Northwest of North America. We used two genetic techniques—microsatellite and genotype‐by‐sequencing (GBS)—to compare the within‐region genetic diversity of populations near the northern extent of the species’ range (British Columbia, Canada) to two geographic regions in British Columbia and two in Washington, USA, moving toward the core of the range. Allelic richness and heterozygosity declined substantially as latitude increased. The northernmost region had the lowest mean expected heterozygosities for both techniques (microsatellite, M = 0.20, SE = 0.080; GBS, M = 0.025, SE = 0.0010) and the southernmost region had the highest (microsatellite, M = 0.88, SE = 0.054; GBS, M = 0.20, SE = 0.0029). The northernmost regions (NC and MC) clustered together in population structure models for both genetic techniques. Our discovery of reduced diversity may have important conservation and management implications for population connectivity and the response of A. truei to climate change.

Early queen joining and long-term queen associations in polygyne colonies of an invasive wasp revealed by longitudinal genetic analysis

Invasive social insects rank among the most damaging of terrestrial species. They are responsible for extensive damage and severely threaten the biodiversity of environments where they are introduced. Variation in colony social form commonly occurs in introduced populations of yellowjacket wasps (genus Vespula). In particular, invasive colonies may contain multiple queens (i.e., polygyne) and persist several years, while in the native range, the colonies are usually annual and harbor a single queen (i.e., monogyne). In this study, we used genome-wide loci obtained by double digest restriction site-associated DNA sequencing (RADseq) to investigate the genetic structure and queen turnover in colonies of the western yellowjacket, Vespula pensylvanica, in their introduced range in Hawaii. Of the 27 colonies monitored over four months (October-January), 19 were polygyne and already contained multiple queens on the first day of sampling. Contrary to previous speculation, this finding suggests that polygyny often arises early in the annual colony cycle, before the production of new queens in the fall. Furthermore, polygyne colonies exhibited a prolonged average lifespan relative to those headed by a single queen. As a result, there is no clear window during which colony eradication efforts would be more effective than upon first discovery. The relatedness among nestmate queens was slightly above zero, indicating that these colonies are generally composed of nonrelatives. The queen turnover within each colony was low, and we detected some full-sibling workers sampled up to four months apart. Finally, we did not detect any population structure among colonies, suggesting that queens disperse up to several kilometers. Taken together, our results provide the first insights into the requeening dynamics in this invasive and incipiently polygyne population and illuminate the early establishment of multiple long-lasting queens in these damaging colonies.

Genome-wide analysis reveals associations between climate and regional patterns of adaptive divergence and dispersal in American pikas

Understanding the role of adaptation in species' responses to climate change is important for evaluating the evolutionary potential of populations and informing conservation efforts. Population genomics provides a useful approach for identifying putative signatures of selection and the underlying environmental factors or biological processes that may be involved. Here, we employed a population genomic approach within a space-for-time study design to investigate the genetic basis of local adaptation and reconstruct patterns of movement across rapidly changing environments in a thermally sensitive mammal, the American pika (Ochotona princeps). Using genotypic data at 49,074 single-nucleotide polymorphisms (SNPs), we analyzed patterns of genome-wide diversity, structure, and migration along three independent elevational transects located at the northern extent (Tweedsmuir South Provincial Park, British Columbia, Canada) and core (North Cascades National Park, Washington, USA) of the Cascades lineage. We identified 899 robust outlier SNPs within- and among-transects. Of those annotated to genes with known function, many were linked with cellular processes related to climate stress including ATP-binding, ATP citrate synthase activity, ATPase activity, hormone activity, metal ion-binding, and protein-binding. Moreover, we detected evidence for contrasting patterns of directional migration along transects across geographic regions that suggest an increased propensity for American pikas to disperse among lower elevation populations at higher latitudes where environments are generally cooler. Ultimately, our data indicate that fine-scale demographic patterns and adaptive processes may vary among populations of American pikas, providing an important context for evaluating biotic responses to climate change in this species and other alpine-adapted mammals.

There Is No 'Rule of Thumb': Genomic Filter Settings for a Small Plant Population to Obtain Unbiased Gene Flow Estimates

The application of high-density polymorphic single-nucleotide polymorphisms (SNP) markers derived from high-throughput sequencing methods has heralded plenty of biological questions about the linkages of processes operating at micro- and macroevolutionary scales. However, the effects of SNP filtering practices on population genetic inference have received much less attention. By performing sensitivity analyses, we empirically investigated how decisions about the percentage of missing data (MD) and the minor allele frequency (MAF) set in bioinformatic processing of genomic data affect direct (i.e., parentage analysis) and indirect (i.e., fine-scale spatial genetic structure - SGS) gene flow estimates. We focus specifically on these manifestations in small plant populations, and particularly, in the rare tropical plant species Dinizia jueirana-facao, where assumptions implicit to analytical procedures for accurate estimates of gene flow may not hold. Avoiding biases in dispersal estimates are essential given this species is facing extinction risks due to habitat loss, and so we also investigate the effects of forest fragmentation on the accuracy of dispersal estimates under different filtering criteria by testing for recent decrease in the scale of gene flow. Our sensitivity analyses demonstrate that gene flow estimates are robust to different setting of MAF (0.05-0.35) and MD (0-20%). Comparing the direct and indirect estimates of dispersal, we find that contemporary estimates of gene dispersal distance (σ _{r t} = 41.8 m) was ∼ fourfold smaller than the historical estimates, supporting the hypothesis of a temporal shift in the scale of gene flow in D. jueirana-facao, which is consistent with predictions based on recent, dramatic forest fragmentation process. While we identified settings for filtering genomic data to avoid biases in gene flow estimates, we stress that there is no 'rule of thumb' for bioinformatic filtering and that relying on default program settings is not advisable. Instead, we suggest that the approach implemented here be applied independently in each separate empirical study to confirm appropriate settings to obtain unbiased population genetics estimates.

The role of a mechanistic host in maintaining arctic rabies variant distributions: Assessment of functional genetic diversity in Alaskan red fox (Vulpes vulpes)

Populations are exposed to different types and strains of pathogens across heterogeneous landscapes, where local interactions between host and pathogen may present reciprocal selective forces leading to correlated patterns of spatial genetic structure. Understanding these coevolutionary patterns provides insight into mechanisms of disease spread and maintenance. Arctic rabies (AR) is a lethal disease with viral variants that occupy distinct geographic distributions across North America and Europe. Red fox (Vulpes vulpes) are a highly susceptible AR host, whose range overlaps both geographically distinct AR strains and regions where AR is absent. It is unclear if genetic structure exists among red fox populations relative to the presence/absence of AR or the spatial distribution of AR variants. Acquiring these data may enhance our understanding of the role of red fox in AR maintenance/spread and inform disease control strategies. Using a genotyping-by-sequencing assay targeting 116 genomic regions of immunogenetic relevance, we screened for sequence variation among red fox populations from Alaska and an outgroup from Ontario, including areas with different AR variants, and regions where the disease was absent. Presumed neutral SNP data from the assay found negligible levels of neutral genetic structure among Alaskan populations. The immunogenetically-associated data identified 30 outlier SNPs supporting weak to moderate genetic structure between regions with and without AR in Alaska. The outliers included SNPs with the potential to cause missense mutations within several toll-like receptor genes that have been associated with AR outcome. In contrast, there was a lack of genetic structure between regions with different AR variants. Combined, we interpret these data to suggest red fox populations respond differently to the presence of AR, but not AR variants. This research increases our understanding of AR dynamics in the Arctic, where host/disease patterns are undergoing flux in a rapidly changing Arctic landscape, including the continued northward expansion of red fox into regions previously predominated by the arctic fox (Vulpes lagopus).

hiphop: Improved paternity assignment among close relatives using a simple exclusion method for biallelic markers

Assignment of parentage with molecular markers is most difficult when the true parents have close relatives in the adult population. Here, we present an efficient solution to that problem by extending simple exclusion approaches to parentage analysis with single nucleotide polymorphic markers (SNPs). We augmented the previously published homozygote opposite test (hot), which counts mismatches due to the offspring and candidate parent having different homozygous genotypes, with an additional test. In this case, parents homozygous for the same SNP are incompatible with heterozygous offspring (i.e., "Homozygous Identical Parents, Heterozygous Offspring are Precluded": hiphop). We tested this approach in a cooperatively breeding bird, the superb fairy-wren, Malurus cyaneus, where rates of extra-pair paternity are exceptionally high, and where paternity assignment is challenging because breeding males typically have first-order adult relatives in their neighbourhood. Combining the tests and conditioning on the maternal genotype with a set of 1376 autosomal SNPs always allowed us to distinguish a single most likely sire from his relatives, and also to identify cases where the true sire must have been unsampled. In contrast, if just the hot test was used, we failed to identify a single most-likely sire in 2.5% of cases. Resampling enabled us to create guidelines for the number of SNPs required when first-order relatives coexist in the mating pool. Our method, implemented in the R package hiphop, therefore provides unambiguous parentage assignments even in systems with complex social organisation. We also identified a suite of Z- and W-linked SNPs that always identified sex correctly.

Genome-wide SNPs detect no evidence of genetic population structure for reef manta rays (Mobula alfredi) in southern Mozambique

Little is known about the extent of genetic connectivity along continuous coastlines in manta rays, or whether site visitation is influenced by relatedness. Such information is pertinent to defining population boundaries and understanding localized dispersal patterns and behaviour. Here, we use 3057 genome-wide single-nucleotide polymorphisms (SNPs) to evaluate population genetic structure and assess the levels of relatedness at aggregation sites of reef manta rays (Mobula alfredi) in southern Mozambique (n = 114). Contrary to indications of limited dispersal along the southern Mozambican coastline inferred from photo-identification and telemetry studies, our results show no evidence of population structure (non-significant F_ST < 0.001) for M. alfredi along this coast. We also found no evidence that individuals sampled at the same site were more related than expected by chance for males, females or across both sexes, suggesting that kinship may not influence visitation patterns at these sites. We estimated the effective population size (N_e) of this population to be 375 (95% CI = 369-380). Comparison to a distant eastern Indian Ocean site (Western Australia, n = 15) revealed strong genetic differentiation between Mozambique and Western Australia (F_ST = 0.377), identifying the Indian Ocean basin as a barrier to dispersal. Our findings show that genetic connectivity in M. alfredi extends for several hundred kilometres along continuous coastlines. We therefore recommend that the population in Mozambique be considered a discrete management unit, and future conservation plans should prioritize integrated strategies along the entire southern coastline.

Genetic assessment of a conservation breeding program of the houbara bustard (Chlamydotis undulata undulata) in Morocco, based on pedigree and molecular analyses

Protection and restoration of species in the wild may require conservation breeding programs under genetic management to minimize deleterious effects of genetic changes that occur in captivity, while preserving populations' genetic diversity and evolutionary resilience. Here, through interannual pedigree analyses, we first assessed the efficiency of a 21-year genetic management, including minimization of mean kinship, inbreeding avoidance, and regular addition of founders, of a conservation breeding program targeting on Houbara bustard (Chlamydotis undulata undulata) in Morocco. Secondly, we compared pedigree analyses, the classical way of assessing and managing genetic diversity in captivity, to molecular analyses based on seven microsatellites. Pedigree-based results indicated an efficient maintenance of the genetic diversity (99% of the initial genetic diversity retained) while molecular-based results indicated an increase in allelic richness and an increase in unbiased expected heterozygosity across time. The pedigree-based average inbreeding coefficient F remained low (between 0.0004 and 0.003 in 2017) while the proportion of highly inbred individuals (F > .1) decreased over time and reached 0.2% in 2017. Furthermore, pedigree-based F and molecular-based individual multilocus heterozygosity were weakly negatively correlated, (Pearson's r = -.061 when considering all genotyped individuals), suggesting that they cannot be considered as alternatives, but rather as complementary sources of information. These findings suggest that a strict genetic monitoring and management, based on both pedigree and molecular tools can help mitigate genetic changes and allow to preserve genetic diversity and evolutionary resilience in conservation breeding programs.

A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Parentage and relatedness reconstruction in Pinus sylvestris using genotyping-by-sequencing

Estimating kinship is fundamental for studies of evolution, conservation, and breeding. Genotyping-by-sequencing (GBS) and other restriction based genotyping methods have become widely applied in these applications in non-model organisms. However, sequencing errors, depth, and reproducibility between library preps could potentially hinder accurate genetic inferences. In this study, we tested different sets of parameters in data filtering, different reference populations and eight estimation methods to obtain a robust procedure for relatedness estimation in Scots pine (Pinus sylvestris L.). We used a seed orchard as our study system, where candidate parents are known and pedigree reconstruction can be compared with theoretical expectations. We found that relatedness estimates were lower than expected for all categories of kinship estimated if the proportion of shared SNPs was low. However, estimates reached expected values if loci showing an excess of heterozygotes were removed and genotyping error rates were considered. The genetic variance-covariance matrix (G-matrix) estimation, however, performed poorly in kinship estimation. The reduced relatedness estimates are likely due to false heterozygosity calls. We analyzed the mating structure in the seed orchard and identified a selfing rate of 3% (including crosses between clone mates) and external pollen contamination of 33.6%. Little genetic structure was observed in the sampled Scots pine natural populations, and the degree of inbreeding in the orchard seed crop is comparable to natural stands. We illustrate that under our optimized data processing procedure, relatedness, and genetic composition, including level of pollen contamination within a seed orchard crop, can be established consistently by different estimators.

RAD-sequencing for estimating genomic relatedness matrix-based heritability in the wild: A case study in roe deer

Estimating the evolutionary potential of quantitative traits and reliably predicting responses to selection in wild populations are important challenges in evolutionary biology. The genomic revolution has opened up opportunities for measuring relatedness among individuals with precision, enabling pedigree-free estimation of trait heritabilities in wild populations. However, until now, most quantitative genetic studies based on a genomic relatedness matrix (GRM) have focused on long-term monitored populations for which traditional pedigrees were also available, and have often had access to knowledge of genome sequence and variability. Here, we investigated the potential of RAD-sequencing for estimating heritability in a free-ranging roe deer (Capreolous capreolus) population for which no prior genomic resources were available. We propose a step-by-step analytical framework to optimize the quality and quantity of the genomic data and explore the impact of the single nucleotide polymorphism (SNP) calling and filtering processes on the GRM structure and GRM-based heritability estimates. As expected, our results show that sequence coverage strongly affects the number of recovered loci, the genotyping error rate and the amount of missing data. Ultimately, this had little effect on heritability estimates and their standard errors, provided that the GRM was built from a minimum number of loci (above 7,000). Genomic relatedness matrix-based heritability estimates thus appear robust to a moderate level of genotyping errors in the SNP data set. We also showed that quality filters, such as the removal of low-frequency variants, affect the relatedness structure of the GRM, generating lower h² estimates. Our work illustrates the huge potential of RAD-sequencing for estimating GRM-based heritability in virtually any natural population.

Quality and quantity of genetic relatedness data affect the analysis of social structure

Kinship plays a fundamental role in the evolution of social systems and is considered a key driver of group living. To understand the role of kinship in the formation and maintenance of social bonds, accurate measures of genetic relatedness are critical. Genotype-by-sequencing technologies are rapidly advancing the accuracy and precision of genetic relatedness estimates for wild populations. The ability to assign kinship from genetic data varies depending on a species' or population's mating system and pattern of dispersal, and empirical data from longitudinal studies are crucial to validate these methods. We use data from a long-term behavioural study of a polygynandrous, bisexually philopatric marine mammal to measure accuracy and precision of parentage and genetic relatedness estimation against a known partial pedigree. We show that with moderate but obtainable sample sizes of approximately 4,235 SNPs and 272 individuals, highly accurate parentage assignments and genetic relatedness coefficients can be obtained. Additionally, we subsample our data to quantify how data availability affects relatedness estimation and kinship assignment. Lastly, we conduct a social network analysis to investigate the extent to which accuracy and precision of relatedness estimation improve statistical power to detect an effect of relatedness on social structure. Our results provide practical guidance for minimum sample sizes and sequencing depth for future studies, as well as thresholds for post hoc interpretation of previous analyses.

Genome-wide association study of an unusual dolphin mortality event reveals candidate genes for susceptibility and resistance to cetacean morbillivirus

Infectious diseases are significant demographic and evolutionary drivers of populations, but studies about the genetic basis of disease resistance and susceptibility are scarce in wildlife populations. Cetacean morbillivirus (CeMV) is a highly contagious disease that is increasing in both geographic distribution and incidence, causing unusual mortality events (UME) and killing tens of thousands of individuals across multiple cetacean species worldwide since the late 1980s. The largest CeMV outbreak in the Southern Hemisphere reported to date occurred in Australia in 2013, where it was a major factor in a UME, killing mainly young Indo-Pacific bottlenose dolphins (Tursiops aduncus). Using cases (nonsurvivors) and controls (putative survivors) from the most affected population, we carried out a genome-wide association study to identify candidate genes for resistance and susceptibility to CeMV. The genomic data set consisted of 278,147,988 sequence reads and 35,493 high-quality SNPs genotyped across 38 individuals. Association analyses found highly significant differences in allele and genotype frequencies among cases and controls at 65 SNPs, and Random Forests conservatively identified eight as candidates. Annotation of these SNPs identified five candidate genes (MAPK8, FBXW11, INADL, ANK3 and ACOX3) with functions associated with stress, pain and immune responses. Our findings provide the first insights into the genetic basis of host defence to this highly contagious disease, enabling the development of an applied evolutionary framework to monitor CeMV resistance across cetacean species. Biomarkers could now be established to assess potential risk factors associated with these genes in other CeMV-affected cetacean populations and species. These results could also possibly aid in the advancement of vaccines against morbilliviruses.

Genetic patterns in peripheral marine populations of the fusilier fish Caesio cuning within the Kuroshio Current

AIM

Mayr's central-peripheral population model (CCPM) describes the marked differences between central and peripheral populations in genetic diversity, gene flow, and census size. When isolation leads to genetic divergence, these peripheral populations have high evolutionary value and can influence biogeographic patterns. In tropical marine species with pelagic larvae, powerful western-boundary currents have great potential to shape the genetic characteristics of peripheral populations at latitudinal extremes. We tested for the genetic patterns expected by the CCPM in peripheral populations that are located within the Kuroshio Current for the Indo-Pacific reef fish, Caesio cuning.

METHODS

We used a panel of 2,677 SNPs generated from restriction site-associated DNA (RAD) sequencing to investigate genetic diversity, relatedness, effective population size, and spatial patterns of population connectivity from central to peripheral populations of C. cuning along the Kuroshio Current.

RESULTS

Principal component and cluster analyses indicated a genetically distinct lineage at the periphery of the C. cuning species range and examination of SNPs putatively under divergent selection suggested potential for local adaptation in this region. We found signatures of isolation-by-distance and significant genetic differences between nearly all sites. Sites closest to the periphery exhibited increased within-population relatedness and decreased effective population size.

MAIN CONCLUSIONS

Despite the potential for homogenizing gene flow along the Kuroshio Current, peripheral populations in C. cuning conform to the predictions of the CCPM. While oceanography, habitat availability, and dispersal ability are all likely to shape the patterns found in C. cuning across this central-peripheral junction, the impacts of genetic drift and natural selection in increasing smaller peripheral populations appear to be probable influences on the lineage divergence found in the Ryukyu Islands.

Whole genome amplification and exome sequencing of archived schistosome miracidia

Adult schistosomes live in the blood vessels and cannot easily be sampled from humans, so archived miracidia larvae hatched from eggs expelled in feces or urine are commonly used for population genetic studies. Large collections of archived miracidia on FTA cards are now available through the Schistosomiasis Collection at the Natural History Museum (SCAN). Here we describe protocols for whole genome amplification of Schistosoma mansoni and Schistosome haematobium miracidia from these cards, as well as real time PCR quantification of amplified schistosome DNA. We used microgram quantities of DNA obtained for exome capture and sequencing of single miracidia, generating dense polymorphism data across the exome. These methods will facilitate the transition from population genetics, using limited numbers of markers to population genomics using genome-wide marker information, maximising the value of collections such as SCAN.

Urban rat races: spatial population genomics of brown rats (Rattus norvegicus) compared across multiple cities

Urbanization often substantially influences animal movement and gene flow. However, few studies to date have examined gene flow of the same species across multiple cities. In this study, we examine brown rats (Rattus norvegicus) to test hypotheses about the repeatability of neutral evolution across four cities: Salvador, Brazil; New Orleans, USA; Vancouver, Canada; and New York City, USA. At least 150 rats were sampled from each city and genotyped for a minimum of 15 000 genome-wide single nucleotide polymorphisms. Levels of genome-wide diversity were similar across cities, but varied across neighbourhoods within cities. All four populations exhibited high spatial autocorrelation at the shortest distance classes (less than 500 m) owing to limited dispersal. Coancestry and evolutionary clustering analyses identified genetic discontinuities within each city that coincided with a resource desert in New York City, major waterways in New Orleans, and roads in Salvador and Vancouver. Such replicated studies are crucial to assessing the generality of predictions from urban evolution, and have practical applications for pest management and public health. Future studies should include a range of global cities in different biomes, incorporate multiple species, and examine the impact of specific characteristics of the built environment and human socioeconomics on gene flow.

From conservation genetics to conservation genomics: a genome-wide assessment of blue whales (Balaenoptera musculus) in Australian feeding aggregations

Genetic datasets of tens of markers have been superseded through next-generation sequencing technology with genome-wide datasets of thousands of markers. Genomic datasets improve our power to detect low population structure and identify adaptive divergence. The increased population-level knowledge can inform the conservation management of endangered species, such as the blue whale (Balaenoptera musculus). In Australia, there are two known feeding aggregations of the pygmy blue whale (B. m. brevicauda) which have shown no evidence of genetic structure based on a small dataset of 10 microsatellites and mtDNA. Here, we develop and implement a high-resolution dataset of 8294 genome-wide filtered single nucleotide polymorphisms, the first of its kind for blue whales. We use these data to assess whether the Australian feeding aggregations constitute one population and to test for the first time whether there is adaptive divergence between the feeding aggregations. We found no evidence of neutral population structure and negligible evidence of adaptive divergence. We propose that individuals likely travel widely between feeding areas and to breeding areas, which would require them to be adapted to a wide range of environmental conditions. This has important implications for their conservation as this blue whale population is likely vulnerable to a range of anthropogenic threats both off Australia and elsewhere.