ABSTRACT
BACKGROUND: This study explored the genetic variability in the New Zealand sheep population for economically important skin traits. Skins were collected at slaughter from two progeny test flocks, resulting in 725 skins evaluated for grain strain, flatness, crust leather strength and overall suitability for shoe leather. DNA profiles collected from skins post-slaughter were matched to individual animals using previously collected high-density genotypes. RESULTS: Considerable phenotypic variation for skin traits was observed, with around 40% of the skins being identified as suitable for high-value shoe leather production. Several key traits associated with leather production, including flatness, tear strength, grain strength and grain strain were found to be moderate to highly heritable (h2 = 0.28-0.82). There were no major significant genome-wide association study (GWAS) peaks associated with many of the traits examined, however, one single-nucleotide polymorphism (SNP) reached significance for the flatness of the skin over the hindquarters. CONCLUSION: This research confirms that suitable lamb skins can be bred for use as high-value shoe leather. While moderately to highly heritable, skin traits in New Zealand lambs appear to be polygenic with no genes of major effect underlaying the traits of interest. Given the complex nature of these traits, the identification and selection of animals with higher-value skins may be enabled by geomic selection. © 2022 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.
Subject(s)
Genome-Wide Association Study , Plant Breeding , Animals , New Zealand , Polymorphism, Single Nucleotide , Sheep/genetics , SkinABSTRACT
The goal of this study was to assess the feasibility of across-country genomic predictions in Norwegian White Sheep (NWS) and New Zealand Composite (NZC) sheep populations with similar development history. Different training populations were evaluated (i.e., including only NWS or NZC, or combining both populations). Predictions were performed using the actual phenotypes (normalized) and the single-step GBLUP via Bayesian inference. Genotyped NWS animals born in 2016 (N = 267) were used to assess the accuracy and bias of genomic estimated breeding values (GEBVs) predicted for birth weight (BW), weaning weight (WW), carcass weight (CW), EUROP carcass classification (EUC), and EUROP fat grading (EUF). The accuracy and bias of GEBVs differed across traits and training population used. For instance, the GEBV accuracies ranged from 0.13 (BW) to 0.44 (EUC) for GEBVs predicted including only NWS, from 0.06 (BW) to 0.15 (CW) when including only NZC, and from 0.10 (BW) to 0.41 (EUC) when including both NWS and NZC animals in the training population. The regression coefficients used to assess the spread of GEBVs (bias) ranged from 0.26 (BW) to 0.64 (EUF) for only NWS, 0.10 (EUC) to 0.52 (CW) for only NZC, and from 0.42 (WW) to 2.23 (EUC) for both NWS and NZC in the training population. Our findings suggest that across-country genomic predictions based on ssGBLUP might be possible for NWS and NZC, especially for novel traits.
Subject(s)
Genome , Genomics , Animals , Bayes Theorem , Genotype , Models, Genetic , New Zealand , Phenotype , Polymorphism, Single Nucleotide , Sheep/geneticsABSTRACT
Facial eczema (FE) is a significant metabolic disease that affects New Zealand ruminants. Ingestion of the mycotoxin sporidesmin leads to liver and bile duct damage, which can result in photosensitisation, reduced productivity and death. Strategies used to manage the incidence and severity of the disease include breeding. In sheep, there is considerable genetic variation in the response to FE. A commercial testing program is available for ram breeders who aim to increase tolerance, determined by the concentration of the serum enzyme, gamma-glutamyltransferase 21 days after a measured sporidesmin challenge (GGT21). Genome-wide association studies were carried out to determine regions of the genome associated with GGT21. Two regions on chromosomes 15 and 24 are reported, which explain 5% and 1% of the phenotypic variance in the response to FE, respectively. The region on chromosome 15 contains the ß-globin locus. Of the significant SNPs in the region, one is a missense variant within the haemoglobin subunit ß (HBB) gene. Mass spectrometry of haemoglobin from animals with differing genotypes at this locus indicated that genotypes are associated with different forms of adult ß-globin. Haemoglobin haplotypes have previously been associated with variation in several health-related traits in sheep and warrant further investigation regarding their role in tolerance to FE in sheep. We show a strategic approach to the identification of regions of importance for commercial breeding programs with a combination of discovery, statistical and biological validation. This study highlights the power of using increased density genotyping for the identification of influential genomic regions, combined with subsequent inclusion on lower density genotyping platforms.
Subject(s)
Eczema/genetics , Genome-Wide Association Study/veterinary , Polymorphism, Single Nucleotide , Quantitative Trait Loci , Sheep Diseases/genetics , Animals , Eczema/blood , Eczema/etiology , Eczema/veterinary , Genome-Wide Association Study/methods , Hemoglobins/genetics , Sheep , Sheep Diseases/blood , Sheep Diseases/etiology , Sporidesmins/toxicity , gamma-Glutamyltransferase/bloodABSTRACT
Robust biomarkers of chronological age have been developed in humans and model mammalian species such as rats and mice using DNA methylation data. The concept of these so-called "epigenetic clocks" has emerged from a large body of literature describing the relationship between genome-wide methylation levels and age. Epigenetic clocks exploit this phenomenon and use small panels of differentially methylated cytosine (CpG) sites to make robust predictions of chronological age, independent of tissue type. Here, we present highly accurate livestock epigenetic clocks for which we have used the custom mammalian methylation array "HorvathMammalMethyl40" to construct the first epigenetic clock for domesticated goat (Capra hircus), cattle (Bos taurus), Red (Cervus elaphus) and Wapiti deer (Cervus canadensis) and composite-breed sheep (Ovis aries). Additionally, we have constructed a 'farm animal clock' for all species included in the study, which will allow for robust predictions to be extended to various breeds/strains. The farm animal clock shows similarly high accuracy to the individual species' clocks (r > 0.97), utilizing only 217 CpG sites to estimate age (relative to the maximum lifespan of the species) with a single mathematical model. We hypothesise that the applications of this livestock clock could extend well beyond the scope of chronological age estimates. Many independent studies have demonstrated that a deviation between true age and clock derived molecular age is indicative of past and/or present health (including stress) status. There is, therefore, untapped potential to utilize livestock clocks in breeding programs as a predictor for age-related production traits.
Subject(s)
Aging , Biomarkers/analysis , CpG Islands , DNA Methylation , Epigenesis, Genetic , Longevity , Stress, Physiological , Animals , Cattle , Deer , Goats , Mice , Phenotype , Rats , SheepABSTRACT
The overall aim of the Ovine FAANG project is to provide a comprehensive annotation of the new highly contiguous sheep reference genome sequence (Oar rambouillet v1.0). Mapping of transcription start sites (TSS) is a key first step in understanding transcript regulation and diversity. Using 56 tissue samples collected from the reference ewe Benz2616, we have performed a global analysis of TSS and TSS-Enhancer clusters using Cap Analysis Gene Expression (CAGE) sequencing. CAGE measures RNA expression by 5' cap-trapping and has been specifically designed to allow the characterization of TSS within promoters to single-nucleotide resolution. We have adapted an analysis pipeline that uses TagDust2 for clean-up and trimming, Bowtie2 for mapping, CAGEfightR for clustering, and the Integrative Genomics Viewer (IGV) for visualization. Mapping of CAGE tags indicated that the expression levels of CAGE tag clusters varied across tissues. Expression profiles across tissues were validated using corresponding polyA+ mRNA-Seq data from the same samples. After removal of CAGE tags with <10 read counts, 39.3% of TSS overlapped with 5' ends of 31,113 transcripts that had been previously annotated by NCBI (out of a total of 56,308 from the NCBI annotation). For 25,195 of the transcripts, previously annotated by NCBI, no TSS meeting stringent criteria were identified. A further 14.7% of TSS mapped to within 50 bp of annotated promoter regions. Intersecting these predicted TSS regions with annotated promoter regions (±50 bp) revealed 46% of the predicted TSS were "novel" and previously un-annotated. Using whole-genome bisulfite sequencing data from the same tissues, we were able to determine that a proportion of these "novel" TSS were hypo-methylated (32.2%) indicating that they are likely to be reproducible rather than "noise". This global analysis of TSS in sheep will significantly enhance the annotation of gene models in the new ovine reference assembly. Our analyses provide one of the highest resolution annotations of transcript regulation and diversity in a livestock species to date.
ABSTRACT
Arctic charr (Salvelinus alpinus) is a species of high economic value for the aquaculture industry, and of high ecological value due to its Holarctic distribution in both marine and freshwater environments. Novel genome sequencing approaches enable the study of population and quantitative genetic parameters even on species with limited or no prior genomic resources. Low coverage genotyping by sequencing (GBS) was applied in a selected strain of Arctic charr in Sweden originating from a landlocked freshwater population. For the needs of the current study, animals from year classes 2013 (171 animals, parental population) and 2017 (759 animals; 13 full sib families) were used as a template for identifying genome wide single nucleotide polymorphisms (SNPs). GBS libraries were constructed using the PstI and MspI restriction enzymes. Approximately 14.5K SNPs passed quality control and were used for estimating a genomic relationship matrix. Thereafter a wide range of analyses were conducted in order to gain insights regarding genetic diversity and investigate the efficiency of the genomic information for parentage assignment and breeding value estimation. Heterozygosity estimates for both year classes suggested a slight excess of heterozygotes. Furthermore, FST estimates among the families of year class 2017 ranged between 0.009 - 0.066. Principal components analysis (PCA) and discriminant analysis of principal components (DAPC) were applied aiming to identify the existence of genetic clusters among the studied population. Results obtained were in accordance with pedigree records allowing the identification of individual families. Additionally, DNA parentage verification was performed, with results in accordance with the pedigree records with the exception of a putative dam where full sib genotypes suggested a potential recording error. Breeding value estimation for juvenile growth through the usage of the estimated genomic relationship matrix clearly outperformed the pedigree equivalent in terms of prediction accuracy (0.51 opposed to 0.31). Overall, low coverage GBS has proven to be a cost-effective genotyping platform that is expected to boost the selection efficiency of the Arctic charr breeding program.
Subject(s)
Fresh Water , Trout , Animals , Chromosome Mapping , Genotype , Sweden , Trout/geneticsABSTRACT
Genotypes are often used to assign parentage in agricultural and ecological settings. Sequencing can be used to obtain genotypes but does not provide unambiguous genotype calls, especially when sequencing depth is low in order to reduce costs. In that case, standard parentage analysis methods no longer apply. A strategy for using low-depth sequencing data for parentage assignment is developed here. It entails the use of relatedness estimates along with a metric termed excess mismatch rate which, for parent-offspring pairs or trios, is the difference between the observed mismatch rate and the rate expected under a model of inheritance and allele reads without error. When more than one putative parent has similar statistics, bootstrapping can provide a measure of the relatedness similarity. Putative parent-offspring trios can be further checked for consistency by comparing the offspring's estimated inbreeding to half the parent relatedness. Suitable thresholds are required for each metric. These methods were applied to a deer breeding operation consisting of two herds of different breeds. Relatedness estimates were more in line with expectation when the herds were analyzed separately than when combined, although this did not alter which parents were the best matches with each offspring. Parentage results were largely consistent with those based on a microsatellite parentage panel with three discordant parent assignments out of 1561. Two models are investigated to allow the parentage metrics to be calculated with non-random selection of alleles. The tools and strategies given here allow parentage to be assigned from low-depth sequencing data.
Subject(s)
Genomics , Genotype , Genotyping Techniques , Pedigree , Algorithms , Alleles , Breeding , Databases, Genetic , Family , Gene Frequency , Genomics/methods , Microsatellite Repeats , Models, Genetic , Sequence Analysis, DNAABSTRACT
Over the past 40 years New Zealand (NZ) aquaculture has grown into a significant primary industry. Tonnage is small on a global scale, but the industry has built an international reputation for the supply of high quality seafood to many overseas markets. Since the early 1990s the industry has recognized the potential gains from selective breeding and the challenge has been to develop programs that can overcome biological obstacles (such as larval rearing and mortality) and operate cost-effectively on a relatively small scale while still providing significant gains in multiple traits of economic value. This paper provides an overview of the current status, and a perspective on genomic technology implementation, for the family based genetic improvement programs established for the two main species farmed in NZ: Chinook (king) salmon (Oncorhynchus tshawytscha) and GreenshellTM mussel (Perna canaliculus). These programs have provided significant benefit to the industry in which we are now developing genomic resources based on genotyping-by-sequencing to complement the breeding programs, enable evaluation of the genetic diversity and identify the potential benefits of genomic selection. This represents an opportunity to increase genetic gain and more effectively utilize the potential for within family selection.
ABSTRACT
Next-generation reduced representation sequencing (RRS) approaches show great potential for resolving the structure of wild populations. However, the population structure of species that have shown rapid demographic recovery following severe population bottlenecks may still prove difficult to resolve due to high gene flow between subpopulations. Here, we tested the effectiveness of the RRS method Genotyping-By-Sequencing (GBS) for describing the population structure of the New Zealand fur seal (NZFS, Arctocephalus forsteri), a species that was heavily exploited by the 19th century commercial sealing industry and has since rapidly recolonized most of its former range from a few isolated colonies. Using 26,026 neutral single nucleotide polymorphisms (SNPs), we assessed genetic variation within and between NZFS colonies. We identified low levels of population differentiation across the species range (<1% of variation explained by regional differences) suggesting a state of near panmixia. Nonetheless, we observed subtle population substructure between West Coast and Southern East Coast colonies and a weak, but significant (p = 0.01), isolation-by-distance pattern among the eight colonies studied. Furthermore, our demographic reconstructions supported severe bottlenecks with potential 10-fold and 250-fold declines in response to Polynesian and European hunting, respectively. Finally, we were able to assign individuals treated as unknowns to their regions of origin with high confidence (96%) using our SNP data. Our results indicate that while it may be difficult to detect population structure in species that have experienced rapid recovery, next-generation markers and methods are powerful tools for resolving fine-scale structure and informing conservation and management efforts.
ABSTRACT
Pneumonia is an important issue for sheep production, leading to reduced growth rate and a predisposition to pleurisy. The objective of this study was to identify loci associated with pneumonic lesions and pleurisy in New Zealand progeny test lambs. The lungs from 3,572 progeny-test lambs were scored for presence and severity of pneumonic lesions and pleurisy at slaughter. Animals were genotyped using the Illumina Ovine Infinium HD SNP BeadChip (606,006 markers). The heritability of lung lesion score and pleurisy were calculated using the genomic relationship matrix, and genome-wide association analyses were conducted using EMMAX and haplotype trend regression. At slaughter, 35% of lambs had pneumonic lesions, with 9% showing lesions on more than half of any individual lobe. The number of lambs recorded as having pleurisy by the processing plants was 9%. Heritability estimates for pneumonic lesions and pleurisy scores adjusted for heteroscedasticity (CPSa and PLEURa) were 0.16 (± 0.03) and 0.05 (± 0.02), respectively. Five single-nucleotide polymorphisms (SNPs) were significantly associated with pneumonic lesions at the genome-wide level, and additional 37 SNPs were suggestively significant. Four SNPs were significantly associated with pleurisy, with an additional 11 SNPs reaching the suggestive level of significance. There were no regions that overlapped between the 2 traits. Multiple SNPs were in regions that contained genes involved in either the DNA damage response or the innate immune response, including several that had previously been reported to have associations with respiratory disease. Both EMMAX and HTR analyses of pleurisy data showed a significant peak on chromosome 2, located downstream from the transcription factor SP3. SP3 activates or suppresses the expression of numerous genes, including several genes with known functions in the immune system. This study identified several SNPs associated with genes involved in both the innate immune response and the response to DNA damage that are associated with pneumonic lesions and pleurisy in lambs at slaughter. Additionally, the identification in sheep of several SNPs within genes that have previously been associated with the respiratory system in cattle, pigs, rats, and mice indicates that there may be common pathways that underlie the response to invasion by respiratory pathogens in multiple species.
Subject(s)
Genome-Wide Association Study/veterinary , Pleurisy/veterinary , Polymorphism, Single Nucleotide/genetics , Sheep Diseases/genetics , Animals , Genetic Predisposition to Disease , Genotype , Haplotypes , Lung/pathology , New Zealand , Phenotype , Pleurisy/genetics , SheepABSTRACT
High-throughput sequencing methods that multiplex a large number of individuals have provided a cost-effective approach for discovering genome-wide genetic variation in large populations. These sequencing methods are increasingly being utilized in population genetic studies across a diverse range of species. Two side-effects of these methods, however, are (1) sequencing errors and (2) heterozygous genotypes called as homozygous due to only one allele at a particular locus being sequenced, which occurs when the sequencing depth is insufficient. Both of these errors have a profound effect on the estimation of linkage disequilibrium (LD) and, if not taken into account, lead to inaccurate estimates. We developed a new likelihood method, GUS-LD, to estimate pairwise linkage disequilibrium using low coverage sequencing data that accounts for undercalled heterozygous genotypes and sequencing errors. Our findings show that accurate estimates were obtained using GUS-LD, whereas underestimation of LD results if no adjustment is made for the errors.
Subject(s)
Algorithms , Genome-Wide Association Study/methods , Linkage Disequilibrium , Sequence Analysis, DNA/standards , Animals , Data Accuracy , Deer/genetics , Genome-Wide Association Study/standards , Genotype , HeterozygoteABSTRACT
BACKGROUND: Investments in genetic selection have played a major role in the New Zealand sheep industry competitiveness. Selection may erode genetic diversity, which is a crucial factor for the success of breeding programs. Better understanding of linkage disequilibrium (LD) and ancestral effective population size (Ne) through quantifying this diversity and comparison between populations allows for more informed decisions with regards to selective breeding taking population genetic diversity into account. The estimation of N e can be determined via genetic markers and requires knowledge of genetic distances between these markers. Single nucleotide polymorphisms (SNP) data from a sample of 12,597 New Zealand crossbred and purebred sheep genotyped with the Illumina Ovine SNP50 BeadChip was used to perform a genome-wide scan of LD and N e . Three methods to estimate genetic distances were investigated: 1) M1: a ratio fixed across the whole genome of one Megabase per centiMorgan; 2) M2: the ratios of genetic distance (using M3, below) over physical distance fixed for each chromosome; and, 3) M3: a genetic map of inter-SNP distances estimated using CRIMAP software (v2.503). RESULTS: The estimates obtained with M2 and M3 showed much less variability between autosomes than those with M1, which tended to give lower N e results and higher LD decay. The results suggest that N e has decreased since the development of sheep breeds in Europe and this reduction in Ne has been accelerated in the last three decades. The N e estimated for five generations in the past ranged from 71 to 237 for Texel and Romney breeds, respectively. A low level of genetic kinship and inbreeding was estimated in those breeds suggesting avoidance of mating close relatives. CONCLUSIONS: M3 was considered the most accurate method to create genetic maps for the estimation of LD and Ne. The findings of this study highlight the history of genetic selection in New Zealand crossbred and purebred sheep and these results will be very useful to understand genetic diversity of the population with respect to genetic selection. In addition, it will help geneticists to identify genomic regions which have been preferentially selected within a variety of breeds and populations.
Subject(s)
Chromosome Mapping , Linkage Disequilibrium , Polymorphism, Single Nucleotide , Sheep/genetics , Animals , Breeding , Genetic Markers , Genome , Population DensityABSTRACT
BACKGROUND: Knowledge about the genetic diversity of a population is a crucial parameter for the implementation of successful genomic selection and conservation of genetic resources. The aim of this research was to establish the scientific basis for the implementation of genomic selection in a composite Terminal sheep breeding scheme by providing consolidated linkage disequilibrium (LD) measures across SNP markers, estimating consistency of gametic phase between breed-groups, and assessing genetic diversity measures, such as effective population size (Ne), and population structure parameters, using a large number of animals (n = 14,845) genotyped with a high density SNP chip (606,006 markers). Information generated in this research will be useful for optimizing molecular breeding values predictions and managing the available genetic resources. RESULTS: Overall, as expected, levels of pairwise LD decreased with increasing distance between SNP pairs. The mean LD r2 between adjacent SNP was 0.26 ± 0.10. The most recent effective population size for all animals (687) and separately per breed-groups: Primera (974), Lamb Supreme (380), Texel (227) and Dual-Purpose (125) was quite variable. The genotyped animals were outbred or had an average low level of inbreeding. Consistency of gametic phase was higher than 0.94 for all breed pairs at the average distance between SNP on the chip (~4.74 kb). Moreover, there was not a clear separation between the breed-groups based on principal component analysis, suggesting that a mixed-breed training population for calculation of molecular breeding values would be beneficial. CONCLUSIONS: This study reports, for the first time, estimates of linkage disequilibrium, genetic diversity and population structure parameters from a genome-wide perspective in New Zealand Terminal Sire composite sheep breeds. The levels of linkage disequilibrium indicate that genomic selection could be implemented with the high density SNP panel. The moderate to high consistency of gametic phase between breed-groups and overlapping population structure support the pooling of the animals in a mixed training population for genomic predictions. In addition, the moderate to high Ne highlights the need to genotype and phenotype a large training population in order to capture most of the haplotype diversity and increase accuracies of genomic predictions. The results reported herein are a first step toward understanding the genomic architecture of a Terminal Sire composite sheep population and for the optimal implementation of genomic selection and genome-wide association studies in this sheep population.
Subject(s)
Oligonucleotide Array Sequence Analysis , Polymorphism, Single Nucleotide , Sheep/genetics , Animals , Female , Genetic Markers/genetics , Genotype , Linkage Disequilibrium , Male , Population DensityABSTRACT
BACKGROUND: New Zealand has some unique Terminal Sire composite sheep breeds, which were developed in the last three decades to meet commercial needs. These composite breeds were developed based on crossing various Terminal Sire and Maternal breeds and, therefore, present high genetic diversity compared to other sheep breeds. Their breeding programs are focused on improving carcass and meat quality traits. There is an interest from the industry to implement genomic selection in this population to increase the rates of genetic gain. Therefore, the main objectives of this study were to determine the accuracy of predicted genomic breeding values for various growth, carcass and meat quality traits using a HD SNP chip and to evaluate alternative genomic relationship matrices, validation designs and genomic prediction scenarios. A large multi-breed population (n = 14,845) was genotyped with the HD SNP chip (600 K) and phenotypes were collected for a variety of traits. RESULTS: The average observed accuracies (± SD) for traits measured in the live animal, carcass, and, meat quality traits ranged from 0.18 ± 0.07 to 0.33 ± 0.10, 0.28 ± 0.09 to 0.55 ± 0.05 and 0.21 ± 0.07 to 0.36 ± 0.08, respectively, depending on the scenario/method used in the genomic predictions. When accounting for population stratification by adjusting for 2, 4 or 6 principal components (PCs) the observed accuracies of molecular breeding values (mBVs) decreased or kept constant for all traits. The mBVs observed accuracies when fitting both G and A matrices were similar to fitting only G matrix. The lowest accuracies were observed for k-means cross-validation and forward validation performed within each k-means cluster. CONCLUSIONS: The accuracies observed in this study support the feasibility of genomic selection for growth, carcass and meat quality traits in New Zealand Terminal Sire breeds using the Ovine HD SNP chip. There was a clear advantage on using a mixed training population instead of performing analyzes per genomic clusters. In order to perform genomic predictions per breed group, genotyping more animals is recommended to increase the size of the training population within each group and the genetic relationship between training and validation populations. The different scenarios evaluated in this study will help geneticists and breeders to make wiser decisions in their breeding programs.
Subject(s)
Breeding , Genomics , Meat , Oligonucleotide Array Sequence Analysis , Polymorphism, Single Nucleotide , Sheep/growth & development , Sheep/genetics , Animals , Female , Genotype , MaleABSTRACT
BACKGROUND: Genotype imputation is a key element of the implementation of genomic selection within the New Zealand sheep industry, but many factors can influence imputation accuracy. Our objective was to provide practical directions on the implementation of imputation strategies in a multi-breed sheep population genotyped with three single nucleotide polymorphism (SNP) panels: 5K, 50K and HD (600K SNPs). RESULTS: Imputation from 5K to HD was slightly better (0.6 %) than imputation from 5K to 50K. Two-step imputation from 5K to 50K and then from 50K to HD outperformed direct imputation from 5K to HD. A slight loss in imputation accuracy was observed when a large fixed reference population was used compared to a smaller within-breed reference (including all 50K genotypes on animals from different breeds excluding those in the validation set i.e. to be imputed), but only for a few animals across all imputation scenarios from 5K to 50K. However, a major gain in imputation accuracy for a large proportion of animals (purebred and crossbred), justified the use of a fixed and large reference dataset for all situations. This study also investigated the loss in imputation accuracy specifically for SNPs located at the ends of each chromosome, and showed that only chromosome 26 had an overall imputation (5K to 50K) accuracy for 100 SNPs at each end higher than 60 % (r2). Most of the chromosomes displayed reduced imputation accuracy at least at one of their ends. Prediction of imputation accuracy based on the relatedness of low-density genotypes to those of the reference dataset, before imputation (without running an imputation software) was also investigated. FIMPUTE V2.2 outperformed BEAGLE 3.3.2 across all imputation scenarios. CONCLUSIONS: Imputation accuracy in sheep breeds can be improved by following a set of recommendations on SNP panels, software, strategies of imputation (one- or two-step imputation), and choice of the animals to be genotyped using both high- and low-density SNP panels. We present a method that predicts imputation accuracy for individual animals at the low-density level, before running imputation, which can be used to restrict genomic prediction only to the animals that can be imputed with sufficient accuracy.
ABSTRACT
BACKGROUND: Genotyping-by-sequencing (GBS) is becoming an attractive alternative to array-based methods for genotyping individuals for a large number of single nucleotide polymorphisms (SNPs). Costs can be lowered by reducing the mean sequencing depth, but this results in genotype calls of lower quality. A common analysis strategy is to filter SNPs to just those with sufficient depth, thereby greatly reducing the number of SNPs available. We investigate methods for estimating relatedness using GBS data, including results of low depth, using theoretical calculation, simulation and application to a real data set. RESULTS: We show that unbiased estimates of relatedness can be obtained by using only those SNPs with genotype calls in both individuals. The expected value of this estimator is independent of the SNP depth in each individual, under a model of genotype calling that includes the special case of the two alleles being read at random. In contrast, the estimator of self-relatedness does depend on the SNP depth, and we provide a modification to provide unbiased estimates of self-relatedness. We refer to these methods of estimation as kinship using GBS with depth adjustment (KGD). The estimators can be calculated using matrix methods, which allow efficient computation. Simulation results were consistent with the methods being unbiased, and suggest that the optimal sequencing depth is around 2-4 for relatedness between individuals and 5-10 for self-relatedness. Application to a real data set revealed that some SNP filtering may still be necessary, for the exclusion of SNPs which did not behave in a Mendelian fashion. A simple graphical method (a 'fin plot') is given to illustrate this issue and to guide filtering parameters. CONCLUSION: We provide a method which gives unbiased estimates of relatedness, based on SNPs assayed by GBS, which accounts for the depth (including zero depth) of the genotype calls. This allows GBS to be applied at read depths which can be chosen to optimise the information obtained. SNPs with excess heterozygosity, often due to (partial) polyploidy or other duplications can be filtered based on a simple graphical method.
Subject(s)
Genotyping Techniques/methods , High-Throughput Nucleotide Sequencing/methods , Sequence Analysis, DNA/methods , Algorithms , Animals , Genotype , Humans , Polymorphism, Single NucleotideABSTRACT
Accurate pedigree information is critical to animal breeding systems to ensure the highest rate of genetic gain and management of inbreeding. The abundance of available genomic data, together with development of high throughput genotyping platforms, means that single nucleotide polymorphisms (SNPs) are now the DNA marker of choice for genomic selection studies. Furthermore the superior qualities of SNPs compared to microsatellite markers allows for standardization between laboratories; a property that is crucial for developing an international set of markers for traceability studies. The objective of this study was to develop a high throughput SNP assay for use in the New Zealand sheep industry that gives accurate pedigree assignment and will allow a reduction in breeder input over lambing. This required two phases of development--firstly, a method of extracting quality DNA from ear-punch tissue performed in a high throughput cost efficient manner and secondly a SNP assay that has the ability to assign paternity to progeny resulting from mob mating. A likelihood based approach to infer paternity was used where sires with the highest LOD score (log of the ratio of the likelihood given parentage to likelihood given non-parentage) are assigned. An 84 "parentage SNP panel" was developed that assigned, on average, 99% of progeny to a sire in a problem where there were 3,000 progeny from 120 mob mated sires that included numerous half sib sires. In only 6% of those cases was there another sire with at least a 0.02 probability of paternity. Furthermore dam information (either recorded, or by genotyping possible dams) was absent, highlighting the SNP test's suitability for paternity testing. Utilization of this parentage SNP assay will allow implementation of progeny testing into large commercial farms where the improved accuracy of sire assignment and genetic evaluations will increase genetic gain in the sheep industry.
Subject(s)
Multiplex Polymerase Chain Reaction/methods , Pedigree , Polymorphism, Single Nucleotide , Sheep/genetics , Animals , New ZealandABSTRACT
DNA-based parentage determination accelerates genetic improvement in sheep by increasing pedigree accuracy. Single nucleotide polymorphism (SNP) markers can be used for determining parentage and to provide unique molecular identifiers for tracing sheep products to their source. However, the utility of a particular "parentage SNP" varies by breed depending on its minor allele frequency (MAF) and its sequence context. Our aims were to identify parentage SNPs with exceptional qualities for use in globally diverse breeds and to develop a subset for use in North American sheep. Starting with genotypes from 2,915 sheep and 74 breed groups provided by the International Sheep Genomics Consortium (ISGC), we analyzed 47,693 autosomal SNPs by multiple criteria and selected 163 with desirable properties for parentage testing. On average, each of the 163 SNPs was highly informative (MAF≥0.3) in 48±5 breed groups. Nearby polymorphisms that could otherwise confound genetic testing were identified by whole genome and Sanger sequencing of 166 sheep from 54 breed groups. A genetic test with 109 of the 163 parentage SNPs was developed for matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. The scoring rates and accuracies for these 109 SNPs were greater than 99% in a panel of North American sheep. In a blinded set of 96 families (sire, dam, and non-identical twin lambs), each parent of every lamb was identified without using the other parent's genotype. In 74 ISGC breed groups, the median estimates for probability of a coincidental match between two animals (PI), and the fraction of potential adults excluded from parentage (PE) were 1.1×10(-39) and 0.999987, respectively, for the 109 SNPs combined. The availability of a well-characterized set of 163 parentage SNPs facilitates the development of high-throughput genetic technologies for implementing accurate and economical parentage testing and traceability in many of the world's sheep breeds.
Subject(s)
Breeding/methods , Pedigree , Polymorphism, Single Nucleotide , Sheep/genetics , Animals , Female , Gene Frequency , Genetic Testing/methods , Genotype , Male , Reproducibility of Results , Sequence Analysis, DNA , Spectrometry, Mass, Matrix-Assisted Laser Desorption-IonizationABSTRACT
* The cpr5-1 Arabidopsis thaliana mutant exhibits constitutive activation of salicylic acid (SA), jasmonic acid (JA) and ethylene (ET) signalling pathways and displays enhanced tolerance of heat stress (HS). * cpr5-1 crossed with jar1-1 (a JA-amino acid synthetase) was compromised in basal thermotolerance, as were the mutants opr3 (mutated in OPDA reductase3) and coi1-1 (affected in an E3 ubiquitin ligase F-box; a key JA-signalling component). In addition, heating wild-type Arabidopsis led to the accumulation of a range of jasmonates: JA, 12-oxophytodienoic acid (OPDA) and a JA-isoleucine (JA-Ile) conjugate. Exogenous application of methyl jasmonate protected wild-type Arabidopsis from HS. * Ethylene was rapidly produced during HS, with levels being modulated by both JA and SA. By contrast, the ethylene mutant ein2-1 conferred greater thermotolerance. * These data suggest that JA acts with SA, conferring basal thermotolerance while ET may act to promote cell death.
Subject(s)
Adaptation, Physiological/drug effects , Arabidopsis/drug effects , Arabidopsis/physiology , Cyclopentanes/pharmacology , Oxylipins/pharmacology , Salicylic Acid/pharmacology , Temperature , Arabidopsis/genetics , Arabidopsis Proteins/metabolism , Cyclopentanes/metabolism , Ethylenes/biosynthesis , Fatty Acids, Unsaturated/metabolism , Gene Expression Regulation, Plant/drug effects , Heat-Shock Response/drug effects , Membrane Proteins/metabolism , Mutation/genetics , Nucleotidyltransferases/metabolism , Oxylipins/metabolism , Phenotype , RNA, Messenger/genetics , RNA, Messenger/metabolism , Signal Transduction/drug effectsABSTRACT
The growth of Arabidopsis plants in chilling conditions could be related to their levels of salicylic acid (SA). Plants with the SA hydroxylase NahG transgene grew at similar rates to Col-0 wild types at 23 degrees C, and growth of both genotypes was slowed by transfer to 5 degrees C. However, at 5 degrees C, NahG plants displayed relative growth rates about one-third greater than Col-0, so that by 2 months NahG plants were typically 2.7-fold larger. This resulted primarily from greater cell expansion in NahG rosette leaves. Specific leaf areas and leaf area ratios remained similar in both genotypes. Net assimilation rates were similar in both genotypes at 23 degrees C, but higher in NahG at 5 degrees C. Chlorophyll fluorescence measurements revealed no PSII photodamage in chilled leaves of either genotype. Col-0 shoots at 5 degrees C accumulated SA, particularly in glucosylated form. SA in NahG shoots showed similar tendencies at 5 degrees C, but at greatly depleted levels. Catechol was not detected as a metabolite of the NahG transgene product. We also examined growth and SA levels in SA signaling and metabolism mutants at 5 degrees C. The partially SA-insensitive npr1 mutant displayed growth intermediate between NahG and Col-0, while the SA-deficient eds5 mutant behaved like NahG. In contrast, the cpr1 mutant at 5 degrees C accumulated very high levels of SA and its growth was much more inhibited than wild type. At both temperatures, cpr1 was the only SA-responsive genotype in which oxidative damage (measured as thiobarbituric acid-reactive substances) was significantly different from wild type.