Development of an SNP‐based parentage verification panel for lovebirds

The genus Agapornis, or lovebirds, are popular pet parrots worldwide. Currently, breeders are dependent on pedigree records as a selection tool as no molecular parentage verification test is available for any of the nine species. The A. roseicollis reference genome was recently assembled. This was followed by the sequencing of the whole genomes of the parents of the reference genome individual at 30× coverage. The parents’ reads were mapped against the reference genome to identify SNPs. Over 1.6 million SNPs, shared between the parents, were discovered using the Genome Analysis Toolkit pipeline. SNPs were filtered to a panel of 480 SNPs based on Genome Analysis Toolkit parameters. The panel of 480 SNPs was genotyped in a population of 960 lovebirds across seven species. A panel of 262 SNPs was compiled that included SNPs successfully amplified across all species. The 262‐SNP panel was reduced based on the observed heterozygosity (H_O) and minor allele frequency (MAF) values per SNP to include the lowest number of SNPs with the highest exclusion power for parentage verification. Two smaller panels consisting of 195 SNPs with MAF and H_O values >0.1 and 40 SNPs with MAF and H_O values >0.3, were constructed. The panels were verified using 43 families from different species with known relationships to evaluate the exclusion power of each panel. The 195 SNP panel with an average exclusion probability of 99.9% and MAF and H_O values >0.1 was proposed as the routine Agapornis parentage verification panel.     

Estimating the inbreeding level and genetic relatedness in an isolated population of critically endangered Sichuan taimen (Hucho Bleekeri) using genome‐wide SNP markers

Sichuan taimen (Hucho bleekeri) is critically endangered fish listed in The Red List of Threatened Species compiled by the International Union for Conservation of Nature (IUCN). Specific locus amplified fragment sequencing (SLAF‐seq)‐based genotyping was performed for Sichuan taimen with 43 yearling individuals from three locations in Taibai River (a tributary of Yangtze River) that has been sequestered from its access to the ocean for more than 30 years since late 1980s. Applying the inbreeding level and genetic relatedness estimation using 15,396 genome‐wide SNP markers, we found that the inbreeding level of this whole isolated population was at a low level (2.6 × 10^?3 ± 0.079), and the means of coancestry coefficients within and between the three sampling locations were all very low (close to 0), too. Genomic differentiation was negatively correlated with the geographical distances between the sampling locations (p < .001), and the 43 individuals could be considered as genetically independent two groups. The low levels of genomic inbreeding and relatedness indicated a relatively large number of sexually mature individuals were involved in reproduction in Taibai River. This study suggested a genomic‐relatedness‐guided breeding and conservation strategy for wild fish species without pedigree information records.     

Validation of genetic polymorphisms on BTA14 associated with carcass trait in a commercial Hanwoo population

The objective of this study was to validate the association of significant SNPs identified from a previous genome‐wide association study with carcass weight (CWT) in a commercial Hanwoo population. We genotyped 13 SNPs located on BTA14 in 867 steers from Korea Hanwoo feedlot bulls. Of these 13 SNPs, five SNPs, namely rs29021868, rs110061498, rs109546980, rs42404006 and rs42303720, were found to be significantly associated (P < 0.001) with CWT. These five significant markers spanned the 24.3 to 29.4 Mb region of BTA14. The most significant marker (rs29021868) for CWT in this study had a 13.07 kg allele substitution effect and accounted for 2.4% of the additive genetic variance in the commercial Hanwoo population. The SNP marker rs109546980 was found to be significantly associated with both CWT (P < 0.001) and eye muscle area (P < 0.001) and could potentially be exploited for marker‐assisted selection in Hanwoo cattle. We also genotyped the ss319607402 variation, which maps to intron2 of PLAG1 gene and which is already reported to be associated with height, to identify any significant association with carcass weight; however, no such association was observed in this Hanwoo commercial population.     

A 34K SNP genotyping array for Populus trichocarpa: Design,application to the study of natural populations and transferability to other Populus species

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost‐effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre‐ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.     

Association between putative functional variants in the PSMB9 gene and risk of melanoma – re‐analysis of published melanoma genome‐wide association studies

To mine possibly hidden causal single‐nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome‐wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer‐free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype–phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.     

Identifying genetic loci controlling neonatal passive transfer of immunity using a hybrid genotyping strategy

Colostrum intake is critical to a piglet's survival and can be measured by precipitating out the γ‐immunoglobulins from serum with ammonium sulfate (immunocrit). Genetic analysis of immunocrits on 5312 piglets indicated that the heritabilities (se) for direct and maternal effects were 0.13 (0.06) and 0.53 (0.08) respectively. To identify QTL for direct genetic effects, piglets with the highest and lowest immunocrits from 470 litters were selected. Six sets of DNA pools were created based on sire of the litter. These 12 DNA pools were applied to Illumina Porcine SNP60 BeadChips. Normalized X and Y values were analyzed. Three different SNP selection methods were used: deviation of the mean from high vs. low pools, the deviation adjusted for variance based on binomial theory and ANOVA. The 25 highest ranking SNPs were selected from each evaluation for further study along with 12 regions selected based on a five‐SNP window approach. Selected SNPs were individually genotyped in the 988 piglets included in pools as well as in 524 piglets that had intermediate immunocrits. Association analyses were conducted fitting an animal model using the estimated genetic parameters. Nineteen SNPs were nominally associated (P < 0.01) with immunocrit values, of which nine remained significant (P < 0.05) after Bonferroni correction, located in 16 genomic regions on 13 chromosomes. In conclusion, the pooling strategy reduced the cost to scan the genome by more than 80% and identified genomic regions associated with a piglet's ability to acquire γ‐immunoglobulin from colostrum. Each method to rank SNPs from the pooled analyses contributed unique validated markers, suggesting that multiple analyses will reveal more QTL than a single analysis.     

Heritability estimates from genomewide relatedness matrices in wild populations: Application to a passerine,using a small sample size

Genomic developments have empowered the investigation of heritability in wild populations directly from genomewide relatedness matrices (GRM). Such GRM‐based approaches can in particular be used to improve or substitute approaches based on social pedigree (PED‐social). However, measuring heritability from GRM in the wild has not been widely applied yet, especially using small samples and in nonmodel species. Here, we estimated heritability for four quantitative traits (tarsus length, wing length, bill length and body mass), using PED‐social, a pedigree corrected by genetic data (PED‐corrected) and a GRM from a small sample (n = 494) of blue tits from natural populations in Corsica genotyped at nearly 50,000 filtered SNPs derived from RAD‐seq. We also measured genetic correlations among traits, and we performed chromosome partitioning. Heritability estimates were slightly higher when using GRM compared to PED‐social, and PED‐corrected yielded intermediate values, suggesting a minor underestimation of heritability in PED‐social due to incorrect pedigree links, including extra‐pair paternity, and to lower information content than the GRM. Genetic correlations among traits were similar between PED‐social and GRM but credible intervals were very large in both cases, suggesting a lack of power for this small data set. Although a positive linear relationship was found between the number of genes per chromosome and the chromosome heritability for tarsus length, chromosome partitioning similarly showed a lack of power for the three other traits. We discuss the usefulness and limitations of the quantitative genetic inferences based on genomic data in small samples from wild populations.     

SNP‐based heritability and genetic architecture of tarsal osteochondrosis in North American Standardbred horses

Osteochondrosis is a common developmental orthopedic disease characterized by a failure of endochondral ossification. Standardbred horses are recognized as being predisposed to tarsal osteochondrosis. Prior heritability estimates for tarsal osteochondrosis in European Standardbreds and related trotting breeds have been based on pedigree data and range from 17–29%. Here, we report on genetic architecture and heritability based on high‐density genotyping data in a cohort of North American Standardbreds (n = 479) stringently phenotyped for tarsal osteochondrosis. Whole‐genome array genotyping data were imputed to ~2 million single nucleotide polymorphisms (SNPs). SNP‐based heritability of osteochondrosis in this population was explained by 2326 SNPs. The majority of these SNPs (86.6%) had small effects, whereas fewer SNPs had moderate or large effects (10% and 2.9% respectively), which is consistent with a polygenic/complex disease. Heritability was estimated at 0.24 ± 0.16 using two methods of restricted maximum likelihood analysis, as implemented in gcta (with and without a weighted relatedness matrix) and ldak software. Estimates were validated using bootstrapping. Heritability estimates were within the range previously reported and suggest that osteochondrosis is moderately heritable but that a significant portion of disease risk is due to environmental factors and/or genotype × environment interactions. Future identification of the genes/variants that have the most impact on disease risk may allow early recognition of high‐risk individuals.     

netview p: a network visualization tool to unravel complex population structure using genome‐wide SNPs

Network‐based approaches are emerging as valuable tools for the analysis of complex genetic structure in wild and captive populations. netview p combines data quality control with the construction of population networks through mutual k‐nearest neighbours thresholds applied to genome‐wide SNPs. The program is cross‐platform compatible, open‐source and efficiently operates on data ranging from hundreds to hundreds of thousands of SNPs. The pipeline was used for the analysis of pedigree data from simulated (n = 750, SNPs = 1279) and captive silver‐lipped pearl oysters (n = 415, SNPs = 1107), wild populations of the European hake from the Atlantic and Mediterranean (n = 834, SNPs = 380) and grey wolves from North America (n = 239, SNPs = 78 255). The population networks effectively visualize large‐ and fine‐scale genetic structure within and between populations, including family‐level structure and relationships. netview p comprises a network‐based addition to other population analysis tools and provides user‐friendly access to a complex network analysis pipeline through implementation in python .     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

SNPs in SNCA,MCCC1, DLG2, GBF1 and MBNL2 are associated with Parkinson's disease in southern Chinese population

Numerous single nucleotide polymorphisms (SNPs), which have been identified as susceptibility factors for Parkinson's disease (PD) as per genome‐wide association studies, have not been fully characterized for PD patients in China. This study aimed to replicate the relationship between 12 novel SNPs of 12 genes and PD risk in southern Chinese population. Twelve SNPs of 12 genes were detected in 231 PD patients and 249 controls, using the SNaPshot technique. Meta‐analysis was used to assess heterogeneity of effect sizes between this study and published data. The impact of SNPs on gene expression was investigated by analysing the SNP‐gene association in the expression quantitative trait loci (eQTL) data sets. rs8180209 of SNCA (allele model: P = .047, OR = 0.77; additive model: P = .047, OR = 0.77), rs2270968 of MCCC1 (dominant model: P = .024, OR = 1.52), rs7479949 of DLG2 (recessive model; P = .019, OR = 1.52), rs10748818 of GBF1 (additive model: P < .001, OR = 0.37), and rs4771268 of MBNL2 (recessive model: P = .003, OR = 0.48) were replicated to be significantly associated with the increased risk of PD. Noteworthy, a meta‐analysis of previous studies suggested rs8180209, rs2270968, rs7479949 and rs4771268 were in line with those of our cohort. Our study replicated five novel functional SNPs in SNCA, MCCC1, DLG2, GBF1 and MBNL2 could be associated with increased risk of PD in southern Chinese population.     

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.     

Evaluation of PAX3 genetic variants and nevus number

The presence of a high nevus number is the strongest phenotypic predictor of melanoma risk. Here, we describe the results of a three‐stage study directed at identifying risk variants for the high nevus phenotype. At the first stage, 263 melanoma cases from Barcelona were genotyped for 223 single‐nucleotide polymorphisms (SNPs) in 39 candidate genes. Seven SNPs in the PAX3 gene were found to be significantly associated with nevus number under the additive model. Next, the associations for seven PAX3 variants were evaluated in 1217 melanoma cases and 475 controls from Leeds; and in 3054 healthy twins from TwinsUK. Associations with high nevus number were detected for rs6754024 (P values < 0.01) in the Barcelona and Leeds datasets and for rs2855268 (P values < 0.01) in the Barcelona and the TwinsUK sets. Associations (P values < 0.001) in the opposite direction were detected for rs7600206 and rs12995399 in the Barcelona and TwinsUK sets. This study suggests that SNPs in PAX3 are associated with nevus number, providing support for PAX3 as a candidate nevus gene. Further studies are needed to examine the role of PAX3 in melanoma susceptibility.     

Estimating population size using single‐nucleotide polymorphism‐based pedigree data

Reliable population estimates are an important aspect of sustainable wildlife management and conservation but can be difficult to obtain for rare and elusive species. Here, we test a new census method based on pedigree reconstruction recently developed by Creel and Rosenblatt (2013). Using a panel of 96 single‐nucleotide polymorphisms (SNPs), we genotyped fecal samples from two Swedish brown bear populations for pedigree reconstruction. Based on 433 genotypes from central Sweden (CS) and 265 from northern Sweden (NS), the population estimates (N = 630 for CS, N = 408 for NS) fell within the 95% CI of the official estimates. The precision and accuracy improved with increasing sampling intensity. Like genetic capture–mark–recapture methods, this method can be applied to data from a single sampling session. Pedigree reconstruction combined with noninvasive genetic sampling may thus augment population estimates, particularly for rare and elusive species for which sampling may be challenging.     

Development of SNP‐genotyping arrays in two shellfish species

Use of SNPs has been favoured due to their abundance in plant and animal genomes, accompanied by the falling cost and rising throughput capacity for detection and genotyping. Here, we present in vitro (obtained from targeted sequencing) and in silico discovery of SNPs, and the design of medium‐throughput genotyping arrays for two oyster species, the Pacific oyster, Crassostrea gigas, and European flat oyster, Ostrea edulis. Two sets of 384 SNP markers were designed for two Illumina GoldenGate arrays and genotyped on more than 1000 samples for each species. In each case, oyster samples were obtained from wild and selected populations and from three‐generation families segregating for traits of interest in aquaculture. The rate of successfully genotyped polymorphic SNPs was about 60% for each species. Effects of SNP origin and quality on genotyping success (Illumina functionality Score) were analysed and compared with other model and nonmodel species. Furthermore, a simulation was made based on a subset of the C. gigas SNP array with a minor allele frequency of 0.3 and typical crosses used in shellfish hatcheries. This simulation indicated that at least 150 markers were needed to perform an accurate parental assignment. Such panels might provide valuable tools to improve our understanding of the connectivity between wild (and selected) populations and could contribute to future selective breeding programmes.     

RAD genotyping reveals fine‐scale genetic structuring and provides powerful population assignment in a widely distributed marine species,the American lobster (Homarus americanus)

Deciphering genetic structure and inferring connectivity in marine species have been challenging due to weak genetic differentiation and limited resolution offered by traditional genotypic methods. The main goal of this study was to assess how a population genomics framework could help delineate the genetic structure of the American lobster (Homarus americanus) throughout much of the species’ range and increase the assignment success of individuals to their location of origin. We genotyped 10 156 filtered SNPs using RAD sequencing to delineate genetic structure and perform population assignment for 586 American lobsters collected in 17 locations distributed across a large portion of the species’ natural distribution range. Our results revealed the existence of a hierarchical genetic structure, first separating lobsters from the northern and southern part of the range (F_CT = 0.0011; P‐value = 0.0002) and then revealing a total of 11 genetically distinguishable populations (mean F_ST = 0.00185; CI: 0.0007–0.0021, P‐value < 0.0002), providing strong evidence for weak, albeit fine‐scale population structuring within each region. A resampling procedure showed that assignment success was highest with a subset of 3000 SNPs having the highest F_ST. Applying Anderson's (Molecular Ecology Resources, 2010, 10, 701) method to avoid ‘high‐grading bias’, 94.2% and 80.8% of individuals were correctly assigned to their region and location of origin, respectively. Lastly, we showed that assignment success was positively associated with sample size. These results demonstrate that using a large number of SNPs improves fine‐scale population structure delineation and population assignment success in a context of weak genetic structure. We discuss the implications of these findings for the conservation and management of highly connected marine species, particularly regarding the geographic scale of demographic independence.     

State‐dependent cooperation in burying beetles: parents adjust their contribution towards care based on both their own and their partner's size

Handicapping experiments on species with biparental care show that a focal parent increases its contribution when its partner is handicapped. Such results are interpreted as evidence for negotiation, whereby each parent adjusts its amount of care to that of its partner. However, it is currently unclear whether the focal parent responds to a change in its handicapped partner's behaviour or state. To address this gap, we conducted an experiment on the burying beetle Nicrophorus vespilloides where we first generated different‐sized males and females by varying the duration of larval development. We then used a 2 × 2 factorial design in which a small or large male was paired with a small or large female. Small females provided less direct care (food provisioning and interactions with larvae) than large females, and both males and females provided less direct care when paired with a small partner. Thus, the focal parent adjusted its contribution towards care based on both its own state and that of its partner. There was also evidence for negotiation between the two parents as the focal parent adjusted its contribution based on the amount of care by its partner. However, there was no evidence that negotiation accounted for how the focal parent responded to its partner's size. Our results have important implications for our understanding of biparental cooperation as they show that each parent adjusts its contribution not only based on the amount of care provided by its partner but also based on its own state and its partner's state.     

The aunt and uncle effect: an empirical evaluation of the confounding influence of full sibs of parents on pedigree reconstruction.

This study used simulations and a known two-generation pedigree of chinook salmon (Oncorhynchus tshawytscha) to evaluate the effect of full sibs of parents on pedigree reconstruction. Parentage analysis was conducted on 100 parent pair-offspring relationships from pedigrees with unrelated (simulation) and related (chinook salmon) candidate parents. Parentage assignment success for the chinook salmon was lower than in the simulated populations. For example, the six most variable loci (mean H(E) = 0.87) provided a mean of 97% unambiguous assignments in the simulated population and 67% unambiguous assignments for the chinook salmon. Estimates of the pairwise relatedness coefficient ((xy)) for most nonexcluded false parents and true parents of chinook salmon offspring exceeded 0.50. These results support the conclusion that closely related candidate parents decrease the power of genetic markers for pedigree reconstruction based on exclusion. Ambiguous parentage may be resolved using single parent- and parent pair-offspring likelihood analysis, however, these methods should be used with caution and they are not replacements for using more loci when many candidate parents are full sibs.     

Relationship type affects the reliability of dispersal distance estimated using pedigree inferences in partially sampled populations: A case study involving invasive American mink in Scotland

Estimating dispersal—a key parameter for population ecology and management—is notoriously difficult. The use of pedigree assignments, aided by likelihood‐based software, has become popular to estimate dispersal rate and distance. However, the partial sampling of populations may produce false assignments. Further, it is unknown how the accuracy of assignment is affected by the genealogical relationships of individuals and is reflected by software‐derived assignment probabilities. Inspired by a project managing invasive American mink (Neovison vison), we estimated individual dispersal distances using inferred pairwise relationships of culled individuals. Additionally, we simulated scenarios to investigate the accuracy of pairwise inferences. Estimates of dispersal distance varied greatly when derived from different inferred pairwise relationships, with mother–offspring relationship being the shortest (average = 21 km) and the most accurate. Pairs assigned as maternal half‐siblings were inaccurate, with 64%–97% falsely assigned, implying that estimates for these relationships in the wild population were unreliable. The false assignment rate was unrelated to the software‐derived assignment probabilities at high dispersal rates. Assignments were more accurate when the inferred parents were older and immigrants and when dispersal rates between subpopulations were low (1% and 2%). Using 30 instead of 15 loci increased pairwise reliability, but half‐sibling assignments were still inaccurate (>59% falsely assigned). The most reliable approach when using inferred pairwise relationships in polygamous species would be not to use half‐sibling relationship types. Our simulation approach provides guidance for the application of pedigree inferences under partial sampling and is applicable to other systems where pedigree assignments are used for ecological inference.