Evaluation of sequencing strategies for whole-genome imputation with hybrid peeling.

BACKGROUND: For assembling large whole-genome sequence datasets for routine use in research and breeding, the sequencing strategy should be adapted to the methods that will be used later for variant discovery and imputation. In this study, we used simulation to explore the impact that the sequencing strategy and level of sequencing investment have on the overall accuracy of imputation using hybrid peeling, a pedigree-based imputation method that is well suited for large livestock populations. METHODS: We simulated marker array and whole-genome sequence data for 15 populations with simulated or real pedigrees that had different structures. In these populations, we evaluated the effect on imputation accuracy of seven methods for selecting which individuals to sequence, the generation of the pedigree to which the sequenced individuals belonged, the use of variable or uniform coverage, and the trade-off between the number of sequenced individuals and their sequencing coverage. For each population, we considered four levels of investment in sequencing that were proportional to the size of the population. RESULTS: Imputation accuracy depended greatly on pedigree depth. The distribution of the sequenced individuals across the generations of the pedigree underlay the performance of the different methods used to select individuals to sequence and it was critical for achieving high imputation accuracy in both early and late generations. Imputation accuracy was highest with a uniform coverage across the sequenced individuals of 2× rather than variable coverage. An investment equivalent to the cost of sequencing 2% of the population at 2× provided high imputation accuracy. The gain in imputation accuracy from additional investment decreased with larger populations and higher levels of investment. However, to achieve the same imputation accuracy, a proportionally greater investment must be used in the smaller populations compared to the larger ones. CONCLUSIONS: Suitable sequencing strategies for subsequent imputation with hybrid peeling involve sequencing ~2% of the population at a uniform coverage 2×, distributed preferably across all generations of the pedigree, except for the few earliest generations that lack genotyped ancestors. Such sequencing strategies are beneficial for generating whole-genome sequence data in populations with deep pedigrees of closely related individuals.

Accuracy of whole-genome sequence imputation using hybrid peeling in large pedigreed livestock populations.

BACKGROUND: The coupling of appropriate sequencing strategies and imputation methods is critical for assembling large whole-genome sequence datasets from livestock populations for research and breeding. In this paper, we describe and validate the coupling of a sequencing strategy with the imputation method hybrid peeling in real animal breeding settings. METHODS: We used data from four pig populations of different size (18,349 to 107,815 individuals) that were widely genotyped at densities between 15,000 and 75,000 markers genome-wide. Around 2% of the individuals in each population were sequenced (most of them at 1× or 2× and 37-92 individuals per population, totalling 284, at 15-30×). We imputed whole-genome sequence data with hybrid peeling. We evaluated the imputation accuracy by removing the sequence data of the 284 individuals with high coverage, using a leave-one-out design. We simulated data that mimicked the sequencing strategy used in the real populations to quantify the factors that affected the individual-wise and variant-wise imputation accuracies using regression trees. RESULTS: Imputation accuracy was high for the majority of individuals in all four populations (median individual-wise dosage correlation: 0.97). Imputation accuracy was lower for individuals in the earliest generations of each population than for the rest, due to the lack of marker array data for themselves and their ancestors. The main factors that determined the individual-wise imputation accuracy were the genotyping status, the availability of marker array data for immediate ancestors, and the degree of connectedness to the rest of the population, but sequencing coverage of the relatives had no effect. The main factors that determined variant-wise imputation accuracy were the minor allele frequency and the number of individuals with sequencing coverage at each variant site. Results were validated with the empirical observations. CONCLUSIONS: We demonstrate that the coupling of an appropriate sequencing strategy and hybrid peeling is a powerful strategy for generating whole-genome sequence data with high accuracy in large pedigreed populations where only a small fraction of individuals (2%) had been sequenced, mostly at low coverage. This is a critical step for the successful implementation of whole-genome sequence data for genomic prediction and fine-mapping of causal variants.

Pigs Lacking the Scavenger Receptor Cysteine-Rich Domain 5 of CD163 Are Resistant to Porcine Reproductive and Respiratory Syndrome Virus 1 Infection.

Porcine reproductive and respiratory syndrome virus (PRRSV) has a narrow host cell tropism, limited to cells of the monocyte/macrophage lineage. CD163 protein is expressed at high levels on the surface of specific macrophage types, and a soluble form is circulating in blood. CD163 has been described as a fusion receptor for PRRSV, with the scavenger receptor cysteine-rich domain 5 (SRCR5) region having been shown to be the interaction site for the virus. As reported previously, we have generated pigs in which exon 7 of the CD163 gene has been deleted using CRISPR/Cas9 editing in pig zygotes. These pigs express CD163 protein lacking SRCR5 (ΔSRCR5 CD163) and show no adverse effects when maintained under standard husbandry conditions. Not only was ΔSRCR5 CD163 detected on the surface of macrophage subsets, but the secreted, soluble protein can also be detected in the serum of the edited pigs, as shown here by a porcine soluble CD163-specific enzyme-linked immunosorbent assay (ELISA). Previous results showed that primary macrophage cells from ΔSRCR5 CD163 animals are resistant to PRRSV-1 subtype 1, 2, and 3 as well as PRRSV-2 infection in vitro Here, ΔSRCR5 pigs were challenged with a highly virulent PRRSV-1 subtype 2 strain. In contrast to the wild-type control group, ΔSRCR5 pigs showed no signs of infection and no viremia or antibody response indicative of a productive infection. Histopathological analysis of lung and lymph node tissue showed no presence of virus-replicating cells in either tissue. This shows that ΔSRCR5 pigs are fully resistant to infection by the virus.IMPORTANCE Porcine reproductive and respiratory syndrome (PRRS) virus (PRRSV) is the etiological agent of PRRS, causing late-term abortions, stillbirths, and respiratory disease in pigs, incurring major economic losses to the worldwide pig industry. The virus is highly mutagenic and can be divided into two species, PRRSV-1 and PRRSV-2, each containing several subtypes. Current control strategies mainly involve biosecurity measures, depopulation, and vaccination. Vaccines are at best only partially protective against infection with heterologous subtypes and sublineages, and modified live vaccines have frequently been reported to revert to virulence. Here, we demonstrate that a genetic-control approach results in complete resistance to PRRSV infection in vivo CD163 is edited so as to remove the viral interaction domain while maintaining protein expression and biological function, averting any potential adverse effect associated with protein knockout. This research demonstrates a genetic-control approach with potential benefits in animal welfare as well as to the pork industry.

Precision engineering for PRRSV resistance in pigs: Macrophages from genome edited pigs lacking CD163 SRCR5 domain are fully resistant to both PRRSV genotypes while maintaining biological function.

Porcine Reproductive and Respiratory Syndrome (PRRS) is a panzootic infectious disease of pigs, causing major economic losses to the world-wide pig industry. PRRS manifests differently in pigs of all ages but primarily causes late-term abortions and stillbirths in sows and respiratory disease in piglets. The causative agent of the disease is the positive-strand RNA PRRS virus (PRRSV). PRRSV has a narrow host cell tropism, limited to cells of the monocyte/macrophage lineage. CD163 has been described as a fusion receptor for PRRSV, whereby the scavenger receptor cysteine-rich domain 5 (SRCR5) region was shown to be an interaction site for the virus in vitro. CD163 is expressed at high levels on the surface of macrophages, particularly in the respiratory system. Here we describe the application of CRISPR/Cas9 to pig zygotes, resulting in the generation of pigs with a deletion of Exon 7 of the CD163 gene, encoding SRCR5. Deletion of SRCR5 showed no adverse effects in pigs maintained under standard husbandry conditions with normal growth rates and complete blood counts observed. Pulmonary alveolar macrophages (PAMs) and peripheral blood monocytes (PBMCs) were isolated from the animals and assessed in vitro. Both PAMs and macrophages obtained from PBMCs by CSF1 stimulation (PMMs) show the characteristic differentiation and cell surface marker expression of macrophages of the respective origin. Expression and correct folding of the SRCR5 deletion CD163 on the surface of macrophages and biological activity of the protein as hemoglobin-haptoglobin scavenger was confirmed. Challenge of both PAMs and PMMs with PRRSV genotype 1, subtypes 1, 2, and 3 and PMMs with PRRSV genotype 2 showed complete resistance to viral infections assessed by replication. Confocal microscopy revealed the absence of replication structures in the SRCR5 CD163 deletion macrophages, indicating an inhibition of infection prior to gene expression, i.e. at entry/fusion or unpacking stages.

Impact of index hopping and bias towards the reference allele on accuracy of genotype calls from low-coverage sequencing.

BACKGROUND: Inherent sources of error and bias that affect the quality of sequence data include index hopping and bias towards the reference allele. The impact of these artefacts is likely greater for low-coverage data than for high-coverage data because low-coverage data has scant information and many standard tools for processing sequence data were designed for high-coverage data. With the proliferation of cost-effective low-coverage sequencing, there is a need to understand the impact of these errors and bias on resulting genotype calls from low-coverage sequencing. RESULTS: We used a dataset of 26 pigs sequenced both at 2× with multiplexing and at 30× without multiplexing to show that index hopping and bias towards the reference allele due to alignment had little impact on genotype calls. However, pruning of alternative haplotypes supported by a number of reads below a predefined threshold, which is a default and desired step of some variant callers for removing potential sequencing errors in high-coverage data, introduced an unexpected bias towards the reference allele when applied to low-coverage sequence data. This bias reduced best-guess genotype concordance of low-coverage sequence data by 19.0 absolute percentage points. CONCLUSIONS: We propose a simple pipeline to correct the preferential bias towards the reference allele that can occur during variant discovery and we recommend that users of low-coverage sequence data be wary of unexpected biases that may be produced by bioinformatic tools that were designed for high-coverage sequence data.

Potential of gene drives with genome editing to increase genetic gain in livestock breeding programs.

BACKGROUND: This paper uses simulation to explore how gene drives can increase genetic gain in livestock breeding programs. Gene drives are naturally occurring phenomena that cause a mutation on one chromosome to copy itself onto its homologous chromosome. METHODS: We simulated nine different breeding and editing scenarios with a common overall structure. Each scenario began with 21 generations of selection, followed by 20 generations of selection based on true breeding values where the breeder used selection alone, selection in combination with genome editing, or selection with genome editing and gene drives. In the scenarios that used gene drives, we varied the probability of successfully incorporating the gene drive. For each scenario, we evaluated genetic gain, genetic variance [Formula: see text], rate of change in inbreeding ([Formula: see text]), number of distinct quantitative trait nucleotides (QTN) edited, rate of increase in favourable allele frequencies of edited QTN and the time to fix favourable alleles. RESULTS: Gene drives enhanced the benefits of genome editing in seven ways: (1) they amplified the increase in genetic gain brought about by genome editing; (2) they amplified the rate of increase in the frequency of favourable alleles and reduced the time it took to fix them; (3) they enabled more rapid targeting of QTN with lesser effect for genome editing; (4) they distributed fixed editing resources across a larger number of distinct QTN across generations; (5) they focussed editing on a smaller number of QTN within a given generation; (6) they reduced the level of inbreeding when editing a subset of the sires; and (7) they increased the efficiency of converting genetic variation into genetic gain. CONCLUSIONS: Genome editing in livestock breeding results in short-, medium- and long-term increases in genetic gain. The increase in genetic gain occurs because editing increases the frequency of favourable alleles in the population. Gene drives accelerate the increase in allele frequency caused by editing, which results in even higher genetic gain over a shorter period of time with no impact on inbreeding.

Genome edited sheep and cattle.

Genome editing tools enable efficient and accurate genome manipulation. An enhanced ability to modify the genomes of livestock species could be utilized to improve disease resistance, productivity or breeding capability as well as the generation of new biomedical models. To date, with respect to the direct injection of genome editor mRNA into livestock zygotes, this technology has been limited to the generation of pigs with edited genomes. To capture the far-reaching applications of gene-editing, from disease modelling to agricultural improvement, the technology must be easily applied to a number of species using a variety of approaches. In this study, we demonstrate zygote injection of TALEN mRNA can also produce gene-edited cattle and sheep. In both species we have targeted the myostatin (MSTN) gene. In addition, we report a critical innovation for application of gene-editing to the cattle industry whereby gene-edited calves can be produced with specified genetics by ovum pickup, in vitro fertilization and zygote microinjection (OPU-IVF-ZM). This provides a practical alternative to somatic cell nuclear transfer for gene knockout or introgression of desirable alleles into a target breed/genetic line.

Highly efficient targeted chromosome deletions using CRISPR/Cas9.

The CRISPR/Cas9 system has emerged as an intriguing new technology for genome engineering. It utilizes the bacterial endonuclease Cas9 which, when delivered to eukaryotic cells in conjunction with a user-specified small guide RNA (gRNA), cleaves the chromosomal DNA at the target site. Here we show that concurrent delivery of gRNAs designed to target two different sites in a human chromosome introduce DNA double-strand breaks in the chromosome and give rise to targeted deletions of the intervening genomic segment. Predetermined genomic DNA segments ranging from several-hundred base pairs to 1 Mbp can be precisely deleted at frequencies of 1-10%, with no apparent correlation between the size of the deleted fragment and the deletion frequency. The high efficiency of this technique holds promise for large genomic deletions that could be useful in generation of cell and animal models with engineered chromosomes.

Substitution of warthog NF-κB motifs into RELA of domestic pigs is not sufficient to confer resilience to African swine fever virus.

African swine fever virus (ASFV) causes a lethal, haemorrhagic disease in domestic swine that threatens pig production across the globe. Unlike domestic pigs, warthogs, which are wildlife hosts of the virus, do not succumb to the lethal effects of infection. There are three amino acid differences between the sequence of the warthog and domestic pig RELA protein; a subunit of the NF-κB transcription factor that plays a key role in regulating the immune response to infections. Domestic pigs with all 3 or 2 of the amino acids from the warthog RELA orthologue have been generated by gene editing. To assess if these variations confer resilience to ASF we established an intranasal challenge model with a moderately virulent ASFV. No difference in clinical, virological or pathological parameters were observed in domestic pigs with the 2 amino acid substitution. Domestic pigs with all 3 amino acids found in warthog RELA were not resilient to ASF but a delay in onset of clinical signs and less viral DNA in blood samples and nasal secretions was observed in some animals. Inclusion of these and additional warthog genetic traits into domestic pigs may be one way to assist in combating the devastating impact of ASFV.

A Strategy To Exploit Surrogate Sire Technology in Livestock Breeding Programs.

In this work, we performed simulations to develop and test a strategy for exploiting surrogate sire technology in animal breeding programs. Surrogate sire technology allows the creation of males that lack their own germline cells, but have transplanted spermatogonial stem cells from donor males. With this technology, a single elite male donor could give rise to huge numbers of progeny, potentially as much as all the production animals in a particular time period. One hundred replicates of various scenarios were performed. Scenarios followed a common overall structure but differed in the strategy used to identify elite donors and how these donors were used in the product development part. The results of this study showed that using surrogate sire technology would significantly increase the genetic merit of commercial sires, by as much as 6.5 to 9.2 years' worth of genetic gain compared to a conventional breeding program. The simulations suggested that a strategy involving three stages (an initial genomic test followed by two subsequent progeny tests) was the most effective of all the strategies tested. The use of one or a handful of elite donors to generate the production animals would be very different to current practice. While the results demonstrate the great potential of surrogate sire technology there are considerable risks but also other opportunities. Practical implementation of surrogate sire technology would need to account for these.

Mammalian interspecies substitution of immune modulatory alleles by genome editing.

We describe a fundamentally novel feat of animal genetic engineering: the precise and efficient substitution of an agronomic haplotype into a domesticated species. Zinc finger nuclease in-embryo editing of the RELA locus generated live born domestic pigs with the warthog RELA orthologue, associated with resilience to African Swine Fever. The ability to efficiently achieve interspecies allele introgression in one generation opens unprecedented opportunities for agriculture and basic research.

Live pigs produced from genome edited zygotes.

Transcription activator-like effector nuclease (TALEN) and zinc finger nuclease (ZFN) genome editing technology enables site directed engineering of the genome. Here we demonstrate for the first time that both TALEN and ZFN injected directly into pig zygotes can produce live genome edited pigs. Monoallelic as well as heterozygous and homozygous biallelic events were identified, significantly broadening the use of genome editor technology in livestock by enabling gene knockout in zygotes from any chosen mating.

Identification of amplified restriction fragment length polymorphism markers linked to genes controlling boar sperm viability following cryopreservation.

This study investigated two hypotheses: 1) that consistent between-boar variation in frozen semen quality exists and is genetically determined, and 2) molecular markers linked to genes controlling semen freezability can be identified using amplified restriction fragment length polymorphism (AFLP) technology. Five ejaculates were collected from each of 129 boars. Semen was diluted into a commercial freezing buffer (700 mOsm/kg, 3% glycerol) and five straws (0.5 ml) per ejaculate were cryopreserved (to -5 degrees C at 6 degrees C/min, then -5 degrees C to -80 degrees C at 40 degrees C/min). Semen was assessed for percentage of motile cells, motility characteristics (computer-aided semen analysis; CASA), plasma membrane integrity (SYBR-14 positive), and acrosome integrity (positive for fluorescein-labeled peanut agglutinin; PNA). Consistent between-boar variability was detected for postthaw sperm motility (P < 0.01), membrane integrity (P < 0.01), acrosome integrity (P < 0.01), and all CASA characteristics (P < 0.05). There was no significant difference between ejaculates (P > 0.05) or straws (P > 0.05) for any viability assessment. Multivariate pattern analysis of the viability data set highlighted three groups of boars producing spermatozoa with poor, average, and good postthaw recovery (42, 63, and 24 boars, respectively). DNA from Large White boars (n = 22) previously classified as good and poor freezers was screened for AFLP markers. Twenty-eight polymerase chain reaction primer combinations generated 2182 restriction fragment bands, of which 421 were polymorphic. Sixteen candidate genetic markers (P < 0.005) were identified by comparing the AFLP profile with semen freezability using logistic regression analysis. These findings support the hypothesis that there is a genetic basis for variation in postthaw semen quality between individuals, and that AFLP technology may be able to identify molecular markers linked to genes influencing this variation.

A form of albinism in cattle is caused by a tyrosinase frameshift mutation.

We used PCR amplification of cDNA prepared from skin biopsies to determine the full-length protein-coding sequence of tyrosinase ( TYR) in cattle of several coat colors. An insertion of a cytosine was detected in an albino Braunvieh calf, which resulted in a frameshift which caused a premature stop codon at residue 316. This insertion was found in the homozygous state in this calf and the genomic DNA of two related albino calves. All six parents of these calves were heterozygous for this insertion. However, an albino Holstein calf did not have this insertion, nor was any other mutation detected in the partial TYR sequence obtained from the genomic DNA available. Diagnostic genotyping tests were developed to detect this mutation in Braunvieh cattle.

A mapping and evolutionary study of porcine sex chromosome genes.

A combination of FISH and RH mapping was used to study the evolution of sex chromosome genes in the pig. In total, 19 genes were identified, including 3 PAR genes (STS, KAL, PRK). The gene order of the porcine X Chromosome (Chr) closely resembled the human X Chr (PRK/STS/KAL-AMELX-EIF2s3X/ZFX-USP9X-DBX-SMCX), suggesting that the porcine X has undergone very little rearrangement during evolution. For the porcine Y Chr, two linkage groups of 10 NRY genes were found, and the following order was established: Ypter-(AMELY-EIF2S3Y/ZFY-USP9Y-DBY/UTY)-(TSPY-SMCY-UBE1Y-SRY)-CEN. This gene order showed greater conservation with the murine Y than with the human Y Chr. In addition, all porcine Y Chr genes mapped to Yp, which is similar to the mouse and included EIF2s3Y and UBE1Y, which are not present in humans. Interestingly, complete conservation of X/Y homologous gene order was found between the pig X and Y Chrs, indicating that the porcine Y Chr has not undergone extensive reorganisation with respect to the X. This suggests that the order of the X/Y homologous genes of the porcine X and Y Chrs may closely resemble the ancestral gene order of the eutherian sex chromosomes.