A practical guide to methods of parentage analysis

The use of molecular techniques for parentage analysis has been a booming science for over a decade. The most important technological breakthrough was the introduction of microsatellite markers to molecular ecology, an advance that was accompanied by a proliferation and refinement of statistical techniques for the analysis of parentage data. Over the last several years, we have seen steady progress in a number of areas related to parentage analysis, and the prospects for successful studies continue to improve. Here, we provide an updated guide for scientists interested in embarking on parentage analysis in natural or artificial populations of organisms, with a particular focus on computer software packages that implement various methods of analysis. Our survey of the literature shows that there are a few established methods that perform extremely well in the analysis of most types of parentage studies. However, particular experimental designs or study systems can benefit from some of the less well-known computer packages available. Overall, we find that parentage analysis is feasible and satisfying in most systems, and we try to provide a simple roadmap to help other scientists navigate the confusing topography of statistical techniques.     

New SNPs for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus)

Sustainable management of sea mullet (Mugil cephalus) fisheries needs to account for recent observations of regional‐scale differentiation. Population genetic analysis is sought to assess the situation of this ecologically and economically important fish species in eastern Australian waters. Here, we report (i) new population genetic markers [single nucleotide polymorphisms (SNPs) and potential microsatellites], (ii) first estimates of spatial genetic differentiation and (iii) prospective power tests for designing more comprehensive studies. Six DNA samples from three sampling regions (North Queensland, South Queensland and central New South Wales) on the eastern coast of Australia were used to prepare restriction site associated DNA (RAD) tag libraries from genomic DNA digested with EcoRI and MseI. A pooled sample of regional RAD tag libraries was sequenced using the Roche GS‐FLX Titanium platform. A total of 172 837 raw reads (17.4 Mbp) were retrieved, 95 500 of which were used to discover 1267 SNPs and 1417 microsatellites. A subset of 161 SNPs was validated based on 63 additional DNA samples genotyped using the Sequenom MassArray (iPLEX Gold chemistry). Altogether 92 SNPs (57%) were confirmed, with 40% of these marking fixed variants between northern and southern sampling regions. Our preliminary findings indicate a multispecies fishery stock of M. cephalus in eastern Australian waters, but suggest that strong genetic differentiation occurs north of major fishing grounds. Low potential differentiation within major fishing grounds (e.g. F_ST = 0.0025) can be resolved with a likely power ≥67% by using standard sample sizes of 50 and validated subsets of available markers.     

Strategies for determining kinship in wild populations using genetic data

Knowledge of kin relationships between members of wild animal populations has broad application in ecology and evolution research by allowing the investigation of dispersal dynamics, mating systems, inbreeding avoidance, kin recognition, and kin selection as well as aiding the management of endangered populations. However, the assessment of kinship among members of wild animal populations is difficult in the absence of detailed multigenerational pedigrees. Here, we first review the distinction between genetic relatedness and kinship derived from pedigrees and how this makes the identification of kin using genetic data inherently challenging. We then describe useful approaches to kinship classification, such as parentage analysis and sibship reconstruction, and explain how the combined use of marker systems with biparental and uniparental inheritance, demographic information, likelihood analyses, relatedness coefficients, and estimation of misclassification rates can yield reliable classifications of kinship in groups with complex kin structures. We outline alternative approaches for cases in which explicit knowledge of dyadic kinship is not necessary, but indirect inferences about kinship on a group‐ or population‐wide scale suffice, such as whether more highly related dyads are in closer spatial proximity. Although analysis of highly variable microsatellite loci is still the dominant approach for studies on wild populations, we describe how the long‐awaited use of large‐scale single‐nucleotide polymorphism and sequencing data derived from noninvasive low‐quality samples may eventually lead to highly accurate assessments of varying degrees of kinship in wild populations.     

Unexpected ancestry of Populus seedlings from a hybrid zone implies a large role for postzygotic selection in the maintenance of species

In the context of potential interspecific gene flow, the integrity of species will be maintained by reproductive barriers that reduce genetic exchange, including traits associated with prezygotic isolation or poor performance of hybrids. Hybrid zones can be used to study the importance of different reproductive barriers, particularly when both parental species and hybrids occur in close spatial proximity. We investigated the importance of barriers to gene flow that act early vs. late in the life cycle of European Populus by quantifying the prevalence of homospecific and hybrid matings within a mosaic hybrid zone. We obtained genotypic data for 11 976 loci from progeny and their maternal parents and constructed a Bayesian model to estimate individual admixture proportions and hybrid classes for sampled trees and for the unsampled pollen parent. Matings that included one or two hybrid parents were common, resulting in admixture proportions of progeny that spanned the whole range of potential ancestries between the two parental species. This result contrasts strongly with the distribution of admixture proportions in adult trees, where intermediate hybrids and each of the parental species are separated into three discrete ancestry clusters. The existence of the full range of hybrids in seedlings is consistent with weak reproductive isolation early in the life cycle of Populus. Instead, a considerable amount of selection must take place between the seedling stage and maturity to remove many hybrid seedlings. Our results highlight that high hybridization rates and appreciable hybrid fitness do not necessarily conflict with the maintenance of species integrity.     

Pedigree reconstruction from SNP data: parentage assignment,sibship clustering and beyond

Data on hundreds or thousands of single nucleotide polymorphisms (SNPs) provide detailed information about the relationships between individuals, but currently few tools can turn this information into a multigenerational pedigree. I present the r package sequoia , which assigns parents, clusters half‐siblings sharing an unsampled parent and assigns grandparents to half‐sibships. Assignments are made after consideration of the likelihoods of all possible first‐, second‐ and third‐degree relationships between the focal individuals, as well as the traditional alternative of being unrelated. This careful exploration of the local likelihood surface is implemented in a fast, heuristic hill‐climbing algorithm. Distinction between the various categories of second‐degree relatives is possible when likelihoods are calculated conditional on at least one parent of each focal individual. Performance was tested on simulated data sets with realistic genotyping error rate and missingness, based on three different large pedigrees (N = 1000–2000). This included a complex pedigree with overlapping generations, occasional close inbreeding and some unknown birth years. Parentage assignment was highly accurate down to about 100 independent SNPs (error rate <0.1%) and fast (<1 min) as most pairs can be excluded from being parent–offspring based on opposite homozygosity. For full pedigree reconstruction, 40% of parents were assumed nongenotyped. Reconstruction resulted in low error rates (<0.3%), high assignment rates (>99%) in limited computation time (typically <1 h) when at least 200 independent SNPs were used. In three empirical data sets, relatedness estimated from the inferred pedigree was strongly correlated to genomic relatedness.     

Identification of high utility SNPs for population assignment and traceability purposes in the pig using high-throughput sequencing

The objectives of this study were to develop breed-specific single nucleotide polymorphisms (SNPs) in five pig breeds sequenced with Illumina's Genome Analyzer and to investigate their usefulness for breed assignment purposes. DNA pools were prepared for Duroc, Landrace, Large White, Pietrain and Wild Boar. The total number of animals used for sequencing was 153. SNP discovery was performed by aligning the filtered reads against Build 7 of the pig genome. A total of 313,964 high confidence SNPs were identified and analysed for the presence of breed-specific SNPs (defined in this context as SNPs for which one of the alleles was detected in only one breed). There were 29,146 putative breed-specific SNPs identified, of which 4441 were included in the PorcineSNP60 beadchip. Upon re-examining the genotypes obtained using the beadchip, 193 SNPs were confirmed as being breed specific. These 193 SNPs were subsequently used to assign an additional 490 individuals from the same breeds, using the sequenced individuals as reference populations. In total, four breed assignment tests were performed. Results showed that for all methods tested 99% of the animals were correctly assigned, with an average probability of assignment of at least 99.2%, indicating the high utility of breed-specific markers for breed assignment and traceability. This study provides a blueprint for the way next-generation sequencing technologies can be used for the identification of breed-specific SNPs, as well as evidence that these SNPs may be a powerful tool for breed assignment and traceability of animal products to their breeds of origin.     

From microsatellites to single nucleotide polymorphisms for the genetic monitoring of a critically endangered sturgeon

The use of genetic information is crucial in conservation programs for the establishment of breeding plans and for the evaluation of restocking success. Short tandem repeats (STRs) have been the most widely used molecular markers in such programs, but next‐generation sequencing approaches have prompted the transition to genome‐wide markers such as single nucleotide polymorphisms (SNPs). Until now, most sturgeon species have been monitored using STRs. The low diversity found in the critically endangered European sturgeon (Acipenser sturio), however, makes its future genetic monitoring challenging, and the current resolution needs to be increased. Here, we describe the discovery of a highly informative set of 79 SNPs using double‐digest restriction‐associated DNA (ddRAD) sequencing and its validation by genotyping using the MassARRAY system. Comparing with STRs, the SNP panel proved to be highly efficient and reproducible, allowing for more accurate parentage and kinship assignments' on 192 juveniles of known pedigree and 40 wild‐born adults. We explore the effectiveness of both markers to estimated relatedness and inbreeding, using simulated and empirical datasets. Interestingly, we found significant correlations between STRs and SNPs at individual heterozygosity and inbreeding that give support to a reasonable representation of whole genome diversity for both markers. These results are useful for the conservation program of A. sturio in building a comprehensive studbook, which will optimize conservation strategies. This approach also proves suitable for other case studies in which highly discriminatory genetic markers are needed to assess parentage and kinship.     

Reviewing the somatic genetics of melanoma: from current to future analytical approaches

Metastatic melanoma has traditionally been difficult to treat, and although molecularly based targeted therapies have shown promising results, they have yet to show consistent improvements in overall survival rates. Thus, identifying the key mutation events underlying the etiology of metastatic melanoma will no doubt lead to the improvement of existing therapeutic approaches and the development of new treatment strategies. Significant advances toward understanding the complexity of the melanoma genome have recently been achieved using next-generation sequencing (NGS) technologies. However, identifying those mutations driving tumorigenesis will continue to be a challenge for researchers, in part because of the high rates of mutation compared to other cancers. This article will review the catalog of mutations identified in melanoma through a variety of approaches, including the use of unbiased exome and whole-genome NGS platforms, as well discuss complementary strategies for identifying driver mutations. The promise of personalized medicine afforded by better understanding these mutation events should provide impetus for increased activity and rapid advances in this field.     

A bayesian framework for parentage analysis: the value of genetic and other biological data

We develop fractional allocation models and confidence statistics for parentage analysis in mating systems. The models can be used, for example, to estimate the paternities of candidate males when the genetic mother is known or to calculate the parentage of candidate parent pairs when neither is known. The models do not require two implicit assumptions made by previous models, assumptions that are potentially erroneous. First, we provide formulas to calculate the expected parentage, as opposed to using a maximum likelihood algorithm to calculate the most likely parentage. The expected parentage is superior as it does not assume a symmetrical probability distribution of parentage and therefore, unlike the most likely parentage, will be unbiased. Second, we provide a mathematical framework for incorporating additional biological data to estimate the prior probability distribution of parentage. This additional biological data might include behavioral observations during mating or morphological measurements known to correlate with parentage. The value of multiple sources of information is increased accuracy of the estimates. We show that when the prior probability of parentage is known, and the expected parentage is calculated, fractional allocation provides unbiased estimates of the variance in reproductive success, thereby correcting a problem that has previously plagued parentage analyses. We also develop formulas to calculate the confidence interval in the parentage estimates, thus enabling the assessment of precision. These confidence statistics have not previously been available for fractional models. We demonstrate our models with several biological examples based on data from two fish species that we study, coho salmon (Oncorhychus kisutch) and bluegill sunfish (Lepomis macrochirus). In coho, multiple males compete to fertilize a single female's eggs. We show how behavioral observations taken during spawning can be combined with genetic data to provide an accurate calculation of each male's paternity. In bluegill, multiple males and multiple females may mate in a single nest. For a nest, we calculate the fertilization success and the 95% confidence interval of each candidate parent pair.     

DNA degradation in fish: Practical solutions and guidelines to improve DNA preservation for genomic research

The more demanding requirements of DNA preservation for genomic research can be difficult to meet when field conditions limit the methodological approaches that can be used or cause samples to be stored in suboptimal conditions. Such limitations may increase rates of DNA degradation, potentially rendering samples unusable for applications such as genome‐wide sequencing. Nonetheless, little is known about the impact of suboptimal sampling conditions. We evaluated the performance of two widely used preservation solutions (1. DESS: 20% DMSO, 0.25 M EDTA, NaCl saturated solution, and 2. Ethanol >99.5%) under a range of storage conditions over a three‐month period (sampling at 1 day, 1 week, 2 weeks, 1 month, and 3 months) to provide practical guidelines for DNA preservation. DNA degradation was quantified as the reduction in average DNA fragment size over time (DNA fragmentation) because the size distribution of DNA segments plays a key role in generating genomic datasets. Tissues were collected from a marine teleost species, the Australasian snapper, Chrysophrys auratus. We found that the storage solution has a strong effect on DNA preservation. In DESS, DNA was only moderately degraded after three months of storage while DNA stored in ethanol showed high levels of DNA degradation already within 24 hr, making samples unsuitable for next‐generation sequencing. Here, we conclude that DESS was the most promising solution when storing samples for genomic applications. We recognize that the best preservation protocol is highly dependent on the organism, tissue type, and study design. We highly recommend performing similar experiments before beginning a study. This study highlights the importance of testing sample preservation protocols and provides both practical and economical advice to improve DNA preservation when sampling for genome‐wide applications.     

Extending RAD tag analysis to microbial ecology: a comparison between MultiLocus Sequence Typing and 2b‐RAD to investigate Listeria monocytogenes genetic structure

The advent of next‐generation sequencing (NGS) has dramatically changed bacterial typing technologies, increasing our ability to differentiate bacterial isolates. Despite it is now possible to sequence a bacterial genome in a few days and at reasonable costs, most genetic analyses do not require whole‐genome sequencing, which also remains impractical for large population samples due to the cost of individual library preparation and bioinformatics. More traditional sequencing approaches, however, such as MultiLocus Sequence Typing (mlst ) are quite laborious and time‐consuming, especially for large‐scale analyses. In this study, a genotyping approach based on restriction site‐associated (RAD) tag sequencing, 2b‐RAD, was applied to characterize Listeria monocytogenes strains. To verify the feasibility of the method, an in silico analysis was performed on 30 available complete genomes. For the same set of strains, in silico mlst analysis was conducted as well. Subsequently, 2b‐RAD and mlst analyses were experimentally carried out on 58 isolates collected from food samples or food‐processing sites. The obtained results demonstrate that 2b‐RAD predicts mlst types and often provides more detailed information on population structure than mlst . Moreover, the majority of variants differentiating identical sequence type isolates mapped against accessory fragments, thus providing additional information to characterize strains. Although mlst still represents a reliable typing method, large‐scale studies on molecular epidemiology and public health, as well as bacterial phylogenetics, population genetics and biosafety could benefit of a low cost and fast turnaround time approach such as the 2b‐RAD analysis proposed here.     

Challenges in studying genomic structural variant formation mechanisms: the short-read dilemma and beyond

Next-generation sequencing (NGS) technologies have revolutionised the analysis of genomic structural variants (SVs), providing significant insights into SV de novo formation based on analyses of rearrangement breakpoint junctions. The short DNA reads generated by NGS, however, have also created novel obstacles by biasing the ascertainment of SVs, an aspect that we refer to as the 'short-read dilemma'. For example, recent studies have found that SVs are often complex, with SV formation generating large numbers of breakpoints in a single event (multi-breakpoint SVs) or structurally polymorphic loci having multiple allelic states (multi-allelic SVs). This complexity may be obscured in short reads, unless the data is analysed and interpreted within its wider genomic context. We discuss how novel approaches will help to overcome the short-read dilemma, and how integration of other sources of information, including the structure of chromatin, may help in the future to deepen the understanding of SV formation processes.     

From barcodes to genomes: extending the concept of DNA barcoding

DNA barcoding has had a major impact on biodiversity science. The elegant simplicity of establishing massive scale databases for a few barcode loci is continuing to change our understanding of species diversity patterns, and continues to enhance human abilities to distinguish among species. Capitalizing on the developments of next generation sequencing technologies and decreasing costs of genome sequencing, there is now the opportunity for the DNA barcoding concept to be extended to new kinds of genomic data. We illustrate the benefits and capacity to do this, and also note the constraints and barriers to overcome before it is truly scalable. We advocate a twin track approach: (i) continuation and acceleration of global efforts to build the DNA barcode reference library of life on earth using standard DNA barcodes and (ii) active development and application of extended DNA barcodes using genome skimming to augment the standard barcoding approach.     

Assignment of paternity groups without access to parental genotypes: multiple mating and developmental plasticity in squid

We present a novel approach to investigating sibling relationships and reconstructing parental genotypes from a progeny array. The Bayesian method we have employed is flexible and may be applicable to a variety of situations in addition to the one presented here. While mutation rates and breeding population allele frequencies can be taken into account, the model requires relatively few loci and makes few assumptions. Paternity of 270 veined squid (Loligo forbesi) hatchlings from three egg strings collected from one location was assigned using five microsatellite loci. Paternal and maternal genotypes reconstructed for each of the three strings were identical, strongly indicating that a single female produced the strings that were fertilized by the same four males. The proportion of eggs fertilized was not equal between males in all three strings, with male 1 siring most offspring (up to 68% in string 1), through to male 4 siring the least (as low as 2.4% in string 1). Although temperature had a profound effect on incubation time, paternity did not affect this trait at 12 degrees C or 8 degrees C.