Parentage assignment with genotyping‐by‐sequencing data

Abstract In this paper, we evaluate using genotype‐by‐sequencing (GBS) data to perform parentage assignment in lieu of traditional array data. The use of GBS data raises two issues: First, for low‐coverage (e.g., <2×) GBS data, it may not be possible to call the genotype at many loci, a critical first step for detecting opposing homozygous markers. Second, the amount of sequencing coverage may vary across individuals, making it challenging to directly compare the likelihood scores between putative parents. To address these issues, we extend the probabilistic framework of Huisman (Molecular Ecology Resources, 2017, 17, 1009) and evaluate putative parents by comparing their (potentially noisy) genotypes to a series of proposal distributions. These distributions describe the expected genotype probabilities for the relatives of an individual. We assign putative parents as a parent if they are classified as a parent (as opposed to e.g., an unrelated individual), and if the assignment score passes a threshold. We evaluated this method on simulated data and found that (a) high‐coverage (>2×) GBS data performs similarly to array data and requires only a small number of markers to correctly assign parents and (b) low‐coverage GBS data (as low as 0.1×) can also be used, provided that it is obtained across a large number of markers. When analysing the low‐coverage GBS data, we also found a high number of false positives if the true parent is not contained within the list of candidate parents, but that this false positive rate can be greatly reduced by hand tuning the assignment threshold. We provide this parentage assignment method as a standalone program called AlphaAssign.


| 103
WHALEN Et AL. to arrays, particularly for species that may not have a wellestablished reference genome, or where a suitable array has not been developed. However, the performance of using GBS data for parentage assignment-to our knowledge-is not well understood.
The primary challenge for using GBS data is the potentially high uncertainty in the true genotype of an individual based on the observed genetic data. In a GBS platform, a restriction enzyme is used to cut DNA into fragments that are then sequenced (Baird et al., 2008;Elshire et al., 2011). This means that unlike arrays, which produce called genotypes, GBS produces read counts for the reference and alternative alleles. For high-coverage GBS data, the underlying genotype can easily be called from the read counts. For low-coverage GBS data calling genotypes is more difficult, particularly on loci which only receive a few reads. Distinguishing between heterozygous and homozygous loci is particularly challenging. If GBS produces two reads for the reference alleles and zero reads for the alternative allele, this could indicate that the individual is homozygous for the reference allele, or the individual could be heterozygous and their reference allele was sequenced twice. The difficulty in calling homozygous loci makes parentage assignment particularly difficult because many parentage assignment algorithms, either explicitly or implicitly, rely on finding opposing homozygous loci to filter out putative parents. In addition, the lack of opposing homozygous loci may increase false positive rate of parentage assignment if the true parent is not in the list of putative parents, since full sibs or half sibs of the true parent may appear to be more related to the individual than expected by chance (Meagher & Thompson, 1986).
Likelihood-based methods (e.g., Kalinowski, Taper, & Marshall, 2007;Riester et al., 2009) are one solution to handle genetic data with high uncertainty. In a likelihood-based method, parentage assignment is based on the likelihood of an individual's genotype conditioned on the putative parent's genotype. If the genotypes of either the individual or the putative parent cannot be assessed accurately, this likelihood score can be calculated by marginalizing over possible genotypes. Likelihood methods work well in cases where all individuals have the same amount of genetic data (e.g., same number of markers or sequencing coverage), but may break down when individuals are genotyped at a different number of markers or at different coverage levels. An example of this could be two putative parents with array data. Suppose the first putative parent was genotyped at 50 markers that overlap with the child, and the second was genotyped at 1,000 markers that overlap with the child. If both parents were heterozygous at all loci and we assume that the loci are not linked, then the likelihood value for the first parent would be 0.5 50 (each allele having a 50% chance of being transmitted), whereas the likelihood value for the second parent would be 0.5 1000 . These likelihood values are hard to compare against each other because they are calculated on different sets of markers. This problem can be solved by selecting a subset of markers that are genotyped in all putative parents (which may drastically reduce the amount of information available), or using the population allele frequency for the genotype at missing markers (which disadvantages individuals with missing values).
A more appealing solution for GBS data is to instead change the parentage assignment problem into a relationship classification problem. With this framing, the goal of the algorithm is to classify the relationship between each putative parent and the focal individual (e.g., parent, grandparent, sibling, child). A putative parent is then assigned as the parent, if they are classified as a parent, pass an assignment threshold and are the highest scoring parent out of the list of putative parents (Huisman, 2017;Riester et al., 2009). One of the main advantages of this approach is that the classification task (which is able to filter out most putative parents) only relies on the genetic information available for an individual and a putative parent and does not require direct comparison to other putative parents. This property is particularly appealing for GBS data where the amount of information on each individual may differ depending on the genotyping resources spent and the allele frequency of the loci with sequence reads.
In this paper, we extend the parentage assignment method of Huisman (2017) to explicitly handle GBS data. We then evaluated its performance in a simulated animal breeding population. We found that, similar to array data, it is possible to obtain accurate parent assignment with a fairly small number of sequence reads (e.g., 0.1× coverage), but that ruling out false positives is harder, and that a sizeable number of false positives could occur for medium coverage (0.5-2×) GBS data on a large number of linked markers.

| MATERIALS AND METHODS
Here, we describe our approach for parentage assignment with GBS data. This work builds closely on the probabilistic framework of Huisman (2017), but we present the full model for completeness. To assign parents, we first construct a series of proposal distributions for each putative parent based on the genotypes of a focal individual and it's known relatives. These proposal distributions describe the expected genotypes for a relative as a function of their relationship with the focal individual (e.g., parent, full sib of the parent, unrelated). We then classify each putative parent into one of these relationships, and if it is classified as a parent, and the assignment score passes a threshold, we assign it as the parent. If there are multiple possible parents, the highest scoring individual is assigned. Although this algorithm was originally designed in the context of animals, it also works for diploid and allopolyploid plants.
To simplify the language, we assume that we are attempting to assign the father of a focal individual. For a given focal individual i and its mother m, we calculate the probability that the putative parent f is the true father by: where g x is the genotype of individual x, h is the relationship between the focal individual i and the putative parent f, and the denominator is enumerated over the set of possible relationships h′. In the case where the genotypes of the mother are unknown, we assume that her genotype probabilities are derived from Hardy-Weinberg Equilibrium.
In this paper, we consider four possible relationships: that the putative parent is the true father, a full sib of the true father, a half sib of the true father, or unrelated. The conditional probability distributions for alternative relationships can be constructed via the generative framework we provide below. To simplify calculations, we assume that p(h′) is uniform over all possible relationships. In addition, we assume all markers segregate independently allowing p(g f |g i ,g m ,h) to be calculated as the product of the probability of the putative parent's genotype at each marker k: The assumption of the independent assortment of markers greatly decreases the computational burden of performing these likelihood calculations, and is effectively true for sparse genotype data (e.g., <1,000 markers spread across multiple chromosomes). However, this assumption is violated for denser genotype data (e.g., above 5,000 markers), and may inflate likelihood values if there are long shared haplotype segments between the focal individual and the (potentially unrelated) putative parent.
In the case of array data, and particularly GBS data, our assessment of the true genotypes, g f , g i and g m may be noisy. To account for this noise, we marginalize across possible genotypes based on observed genetic data d = (d i , df, d m ): This model requires the calculation of two terms: (a) the genotype probabilities conditional on the observed data p(g x,k |d x,k ) and (b) the proposal distribution for an individual's genotype based on their relationship with the focal individual p(g f,k |g i,k ,g m,k ,h). We outline how to calculate both terms below.

| Evaluating genotype probabilities conditional on the observed data
In this model, we assume that each marker is biallelic and has four possible phased genotypes, aa, aA, Aa, AA. If we observed array data for marker k, d x,k , the conditional probabilities for each genotype g x,k are: where e is the assumed genotyping error rate. This evaluation of individual genotype probabilities differs from Huisman (2017), where it is assumed that errors can only occur between homozygous and heterozygous states (and not between opposing homozygote states) and distinction is not made between two heterozygous genotypes. The genotype probabilities above correspond more closely to those commonly used in peeling (e.g., Whalen, Ros-Freixedes, Wilson, Gorjanc, & Hickey, 2017) and allow inferences to be made even when the genotyping error rate is high.
With observed GBS data for marker k, d x,k , the conditional genotype probabilities are: where e is the sequencing error rate, n ref is the number of sequence reads supporting the reference allele and n alt is the number of sequence reads supporting the alternative allele. The genotype probabilities in Equation 5 do not sum to one, and so the probabilities need to be normalized for each allele. Equation 4 is consistent with previous work on parentage assignment with array data (Huisman, 2017;Kalinowski et al., 2007), while Equation 5 is consistent with previous work on imputation with GBS-like data (Li, Willer, Ding, Scheet, & Abecasis, 2010;VanRaden, Sun, & O'Connell, 2015;Whalen et al., 2017).

| Generating proposal distributions via single locus peeling
We generate proposal distributions p(g f,k |g i,k ,g m,k ,h) for the genotype probabilities of each relationship via single locus peeling (Elston & Stewart, 1971). Single locus peeling provides a rich generative model for estimating the genotype probabilities of un-genotyped relatives based on the genotypes of an individual and a known parent. Although our presentation differs from Huisman (2017), it results in the same distributions. Under this framework, we calculate the genotype probabilities for three relatives: the father, a full-sib of the father and a half-sib of the father. These probabilities are calculated by first estimating the genotype probabilities for the father, peeling up to the paternal grandparents, and finally peeling down to the full sib and the half sib of the father ( Figure 1). Given genetic data on the focal individual d i and a mother d m , we can construct a proposal distribution for the father via: where p(g i |d i ) is given by Equations 4 or 5 above, and T(g i |g f ,g m ) is the probability that the individual inherited genotype g i conditional on their parents having genotypes g f and g m , e.g., T g i = aA|g f = aA,g m = AA = 0.5 (Marshall, Slate, Kruuk, & Pemberton, 2003). Both the summations in Equation 6, and in the equations below, are over all four phased genotype states.
Using Equation 6, we can peel up to construct a joint distribution for the genotypes of the paternal grandparents (g gf ,g gm ): where p(g f |d i ,df) is given in Equation 6, above. We can then peel down to generate the proposal distributions for a full sib and a half sib of the father. The proposal distributions differ in whether the full joint distribution of both grandparents is used (full sib, fs), or if only one of the grandparents is used and the other parent assumed to have genotypes based on Hardy Weinberg Equilibrium (half sib, hs): where p(g null ) represents the probability of having a genotype if that genotype was drawn at random from the population.
The proposal distribution for an unrelated individual simply assumes that their genotypes are drawn at random from the population according to Hardy Weinberg Equilibrium: To assign a parent, we calculated an assignment score for each putative parent: The score will be close to 0 if the individual is unlikely to be the father, and tends towards positive infinity with increasing evidence that the individual is the father. A putative parent was assigned as the true parent if its assignment score was the highest of the putative parents considered, and was higher than a threshold. In the simulations, we used a threshold value of 10. This corresponds to a greater than 99.99% posterior probability that the selected individual is the true parent (under the assumption that the loci are unlinked), and led to a near zero positive rate in simulations with SNP array data. The exact threshold used should be determined by the relationship between the putative parent and the focal individual, the number of markers, the LD between the markers and the GBS coverage level. We return to this point in more detail in the Section 4.
Although the described process may seem computationally intensive, there are two features which simplify calculations. First, because the proposal distributions depend only on the focal individual and its known parent, the proposal distributions only need to be calculated once and can be re-used for all putative parents of the focal individual. Second, peeling can be performed efficiently as a series of tensor operations on the genotypes of focal individual and its known parent, filtered through the inheritance matrix T, which allows us to take advantage of linear algebra libraries.
of the peeling order for the proposal distributions. The arrows represent the direction in which the peeling operations should be performed. Hardy-Weinberg equilibrium is used to generate the genotype distributions for the unrelated individual, the mother of the half sib, and if unknown, the mother's genotype. Although this graphic assumes the mother is known and the father unknown, a symmetric picture could be constructed when the mother is unknown and father known

| Simulated data
The simulated data modelled a livestock population. We initially sampled a set of genomes with 20 chromosomes using the Markovian coalescent simulator MaCS (Chen, Marjoram, & Wall, 2009). For this, we assumed that each chromosome is 10 8 bp long, a per site mutation rate is 2.5 × 10 −8 , a per site recombination rate is 1.0 × 10 −8 , and that effective population size changed over time, based on estimates for the Holstein cattle population (Villa-Angulo et al., 2009). We set the effective population size to 100 in the final generation of the coalescent simulation and to 1,256, 4,350 and 43,500 at, respectively, 1,000, 10,000, and 100,000 generations ago, with linear changes in between. We then used the sampled chromosomes to initiate a population of 1,000 animals with equal sex proportions. We simulated this population for five generations. In each generation, we selected 10 males and mated them at random to 100 females. Each potential focal individual therefore had one true father, four male full sibs of the father, and 45 male half sibs of the father. All individuals were genotyped at 50,000 markers. Subsets of these markers were used in different simulations as described below. Array data were simulated without any errors, due to the low error rate for modern SNP genotyping arrays (<1%; e.g., Kalinowski et al., 2007). In addition to array data, we generated low-coverage GBS data for the last two generations of individuals. We simulated GBS data on the same sites as the array data. This either models the situation where the GBS platform and the genotyping array capture the same loci, or where they capture different loci and only the loci in common are used. We assumed that GBS was performed at coverage levels between 0.1× to 10×. For each coverage level, the number of sequence reads at a given marker was generated via a Poisson distribution with mean equal to the coverage level. Each read randomly sampled one of the two alleles at a marker. The read sampling process also included a small sequencing error rate of 0.1%. We generated the simulated data using the R package AlphaSimR (Gaynor, Gorjanc, Wilson, Money, & Hickey, 2018), which is available at www.alphagenes.roslin.ed.ac.uk/AlphaSimR.

| Scenarios
We evaluated the accuracy of parent assignment for the last generation of 1,000 individuals across four different scenarios. In the first scenario (a), we analysed the accuracy of performing parent assignment when: • the mother was known and genotyped (although this is later relaxed), • all of the male full-and half-sibs along with 50 other individuals (total of 100 potential parents) were putative parents, • and either both the parents and progeny had array data, the parents had array data and the progeny had GBS data, or both the parents and the progeny had GBS data.
These sub-scenarios span a spectrum of possible practical settings. The sub-scenario where the parents had array data but the progeny had GBS data may represent the case where the progeny are initially genotyped with a low-cost GBS platform and any selected parents are re-genotyped with an array, or may represent the case where the putative parents had GBS data, but then individual markers were called at high accuracy using a pedigree (with already established relationships) or population based imputation method. In the remaining scenarios we focused on the case where both parents and progeny had GBS data and analysed (b) the impact of knowing and genotyping the known alternative parent, (c) the impact of restricting the pool of putative parents to either 100 unrelated individuals, 45 half sibs, or the four full sibs, and (d) examined how the false positive rate changed depending on the threshold used for assignment (see below).
In each scenario, we performed three evaluations. To evaluate the overall accuracy, we assumed the true parent was included in the list of putative parents, and evaluated accuracy by the number of times the top parent was the true parent. To evaluate the true positive rate we included the true parent in the list of putative parents, but assigned the top scoring parent only if it passed an assignment threshold. To evaluate the false positive rate, we excluded the true parent from the list of putative parents and counted the number of times the top scoring parent passed the assignment threshold. The first evaluation represented a case where we know the true parent is included in the list of putative parents (e.g., groups of females cohabitating with multiple males or artificial insemination using polyspermic matings). The second and third evaluations were designed to assess performance when we are not sure whether or not the true parent is included in the list of potential parents (e.g., natural service sires or wild populations).

| Software
Parentage assignment was performed using AlphaAssign (http://www.alphagenes.roslin.ed.ac.uk/alphasuite-softwares/alphaassign/) which, implements the described algorithm. AlphaAssign has three run-time parameters: (a) an assumed genotyping error rate for array data, (b) an assumed sequencing error rate for GBS data, and (c) an assignment threshold to determine the required score to assign a putative parent as a parent. Throughout this paper, we assumed a 1% genotyping error rate, a 0.1% sequencing error rate, an assignment threshold of 10 (corresponding to a 99.99% posterior probability that the individual is the true parent, under the assumption of unlinked loci) although we varied the assignment threshold in the final set of simulations.

| Parent assignment with array and GBS data
First, we examined the number of markers required for accurate parentage assignment when both parents and progeny were genotyped with array data. If the true parent was included in the list of putative parents (and an assignment threshold was used), 100 markers were required to obtain 100% parentage assignment accuracy. If the true parent was excluded from the list of putative parents, the false positive rate was less than 0.1% if there were between 50-350 markers, and there were no false positives when there were more than 500 markers.
Unlike array data where the number of markers can be more easily varied, for GBS data the number of markers is usually determined by the choice of restriction enzymes while the amount of coverage obtained on each individual can be varied. Because of this, we focused on the required coverage level to accurately assign parents based on a fixed number of markers. Figure 2 shows the accuracy and false positive rates based on the amount of coverage allocated to each progeny, stratified by the number of markers that this coverage is spread over. Because performance with array data was nearly identical to that with 10× GBS data we did not include array data in Figure 2.
We evaluated the performance of parentage assignment when the parents were genotyped with array data and the progeny were genotyped with GBS data. If the true parent was included in the list of putative parents, a coverage of 0.4× was required to obtain 100% accuracy when there were 50,000 GBS markers. The required coverage increased to 1× for 5,000 markers, and to 2× for 1,000 markers. If the true parent was excluded from the list of putative parents, we found that the false positive rate was less than 0.2% in all cases.
The accuracy of parentage assignment decreased when both the parents and progeny had GBS data. If the true parent was included in the list of putative parents, a coverage of 0.4× was required to obtain 100% accuracy when there were 50,000 GBS markers. The required coverage increased to 2× for 5,000 markers, and to 5× for 1,000 markers. If the true parent was excluded from the list of putative parents, we found that the false positive rate was as high as a 60%. These false positives were clustered on low to medium coverage GBS data (0.1-3×) with a large number of markers (>1,000). Figure 3 stratifies the false positive rate based on whether unrelated individuals, half-sibs of the true parent, or full-sibs of the true parent were included in the list of putative parents. In all cases, we assume both the parents and the offspring had GBS data. In line with expectations, we found a high false positive rate (as high as 60% in some conditions) when only the full-sibs of the true parent were included as putative parents. This decreased to at most 35% when only the half-sibs of the true parent were included and to under 20% when only unrelated individuals were included. As seen previously, most of the false positives occurred when there were a large number of markers and low to medium coverage GBS data. Figure 4 compares the performance of parentage assignment when one of the parents is known and genotyped compared to when neither parent is known or genotyped. In all cases, we assumed that both the parents and offspring had GBS data. We found that having one parent known and genotyped increased the accuracy of parentage assignment and decreased the number of false positives in all cases. The benefit was largest when both the progeny and parents had high-coverage GBS data.

| Controlling false assignments by
modifying the threshold Figure 5 shows the true positive rates and false positive rates for when sequencing resources were spread over 50,000 markers, as a function of the threshold used to assign a putative parent as the parent. We found that, compared to the results in Figures 2 and 3, it was possible to substantially reduce the false positive rate by increasing the assignment threshold, but that the ideal threshold depends on the total coverage. The relationship between the false positive and true positive rate is given as a receiver operating characteristic in Figure 5c.

| Timing
The algorithm took 3 min and 54 s to assign parents for 1,000 progeny, each with 100 putative parents. The progeny and their parents were genotyped using GBS data across 5,000 markers. The algorithm scales linearly with the number of markers and the number of putative parents per individual.

| DISCUSSION
In this paper, we extended the parentage assignment method of Huisman (2017) to account for low-coverage sequence data and analysed the performance of parentage assignment when genotyping is performed via sequencing instead of the traditional genome-wide arrays. We found that high-coverage GBS data (i.e., 10× or higher) has the same performance as array data. We also found that low-coverage GBS data (as low as 0.1×) can be used to perform parentage assignment as long as it is obtained on a sufficiently large number of markers (e.g., 0.5× GBS on 5,000 markers), but that there may be a large number of false assignments if the true parent is not included in the list of putative parents. The number of false positives could be reduced by modifying the threshold used to call assignments. In light of these results, we will discuss (a) the accuracy of parentage assignment, (b) potential extensions to control the false positive rate, and (c) the use of peeling to construct the proposal distributions in more detail.

GBS data
A goal of this work was to quantify the amount of GBS data required to accurately perform parentage assignment. We found that, similar to array data, the total amount of data required is relatively low. For example, when using high-coverage GBS data between 100-200 markers are required to accurately assign parents. This is in line with previous estimates for array data (Fisher et al., 2009;Strucken et al., 2016;Tortereau et al., 2017), where between 50 and 700 markers were required. The differences in the exact number of markers required (100-200 compared to 50-700) is likely due to the structure of the underlying genetic data (i.e., number of chromosomes, minor allele frequency of the markers), and the assumption in this study that one of the parents was already known and genotyped.
In addition to being able to use high-coverage GBS data to perform parentage assignment, we found that low-coverage GBS F I G U R E 2 Parentage assignment performance when array or GBS data was available for the parents and GBS data was available for the progeny. The left panels give the number of correct assignments (for 1,000 progeny) when the true parent was on the list of putative parents and no assignment threshold was used, and the top scoring parent was always assigned. The middle panels give the number of correct assignments when the true parent was on the list of putative parents and assignment threshold was used. The right panels give the number of incorrect assignments when the true parent was excluded from the list of putative parents [Colour figure can be viewed at wileyonlinelibrary.com] data could also be used, provided it was spread across a larger number of markers. The increase in required number of markers is due to the lower information content at an individual loci for low-coverage GBS data, requiring the data to be pooled across a larger number of markers to achieve the same level of accuracy.
The results of this study suggest that GBS data-either high-coverage data on a small number of markers, or low-coverage data on a large number of markers-is an effective alternative to array data for performing parentage assignment. This result is particularly important given the emerging importance of GBS as an alternative for SNP array data, both in species where SNP arrays are available (e.g., De Donato, Peters, Mitchell, Hussain, & Imumorin, 2013;Brouard, Boyle, Ibeagha-Awemu, & Bissonnette, 2017) and in those where SNP arrays have not been constructed (e.g., Robledo, Palaiokostas, Bargelloni, Martínez, & Houston, 2017;Palaiokostas et al., 2018).

| Controlling the false positive rate
During our analysis of low-coverage data, we found an inflation of false positives when both the parents and the progeny had GBS data. These false positives were likely due to the fact that with between 1-3× coverage GBS data we were able to determine that two animals are genetically similar, but were not able to obtain a sufficient number of loci with precisely inferred genotypes to find opposing homozygous loci. This was particularly the case when the list of putative parents included individuals that were highly related to the true parents (i.e., full or half sibs), and when the alternative parent was unknown or ungenotyped.
Consistent with previous work, we found that using a handtuned assignment threshold could reduce the number of false positives (Huisman, 2017;Riester et al., 2009). An alternative approach would be to adaptively determine the assignment threshold via introspection of the underlying data (Grashei, Ødegård, & Meuwissen, 2018). In the majority of the simulations, a fixed threshold of 10 was used, inspired by requiring 99.99% of the posterior probability supporting the hypothesis that the individual was the true parent. As we demonstrate in Figure 5, substantially raising the threshold for assignment could reduce the false-positive rate even for 50,000 markers and low-coverage sequence data, although at the cost of a decreased true-positive rate. The optimal threshold value for assignment depends on the overall sequencing coverage, making it challenging to use a fixed threshold in cases where individuals are sequenced at different coverages. We believe that automating this process is an area for future research, and may depend on the exact breeding programme structure, the exact GBS system deployed (e.g., Baird et al., 2008;Elshire et al., 2011), and reason that parentage information is required.
The false positive rate can also be decreased by increasing the accuracy of the parent genotypes. As was highlighted in Figure 2, if the putative parent and the alternative parent have SNP array data (or high accuracy information at all loci), the false positive rate is close to zero. This could be achieved by either re-genotyping candidate parents with higher coverage GBS or SNP array. Alternatively, if there are other individuals in the population who are genotyped, pedigree or population-based imputation methods might be used to call the genotypes of the putative parents at high accuracy (Browning & Browning, 2007;Li et al., 2010;Meuwissen & Goddard, 2010;Whalen et al., 2017), although care must be taken that this imputation does not produce a large number of errors (Chan, Hamblin, & Jannink, 2016).
Furthermore, we believe that the issue of false parent assignments may be less of an issue in the context of commercial agricultural populations compared to wild populations for two reasons. First, most of the false assignments that we observed were cases where the true parent was not included in the pedigree and a full-or half-sib of the true parent was included and wrongly assigned as a parent. In the context of many animal breeding programmes, the routine use of pairs of sibs as parents may not commonly arise because of explicit efforts to manage diversity and inbreeding (e.g., Woolliams, Berg, Dagnachew, & Meuwissen, 2015). Second, due to the genetic similarity between the full-sib of the true parent and the true parent, using the full-sib of the true parent as a "proxy" parent for the progeny may have limited impact on downstream applications such as estimation of breeding values. Further research is required to quantify the impact of such false positives in downstream applications.

| Constructing proposal distributions via peeling
In this paper, we closely followed the approach of Huisman (2017) for performing parentage assignment, with two differences. First, we modified the genotype probability function to handle sequence data. Second, we recast the construction of proposal distributions for relatives as a series of peeling operations on artificial pedigrees. We believe the later development is of more interest. Peeling provides a rich and computationally efficient framework for estimating the genotypes of a relative based on the genotypes of individuals in an existing pedigree. In this paper, we focused on a small number of possible relationships, but this framework can be easily extended to consider a wider and potentially complex class of relatives (e.g., siblings of the focal individual, cousins of the parent, or grandparents), or could be altered to assess alternative relationships (e.g., performing grandparent assignment instead of parentage assignment). Use of these additional relationship classes may depend on the purpose of a particular application.

| CONCLUSION
In conclusion, we extended the algorithm of Huisman (2017) to perform parentage assignment with sequence data, and evaluated the performance of using low-coverage GBS data for parentage assignment. We found that low-coverage GBS data could be used for accurate parentage assignment, but that there may be concerns with false positives if the true parent is not included on the list of putative parents. Such false positives might be mitigated on a case-by-case basis by tuning the assignment criteria used. These results suggest that GBS data can be used as an alternative to array data for parentage assignment.