Using triallelic SNPs for determining parentage in North American yak ( Bos grunniens) and estimating cattle ( B. taurus) introgression.

Background: Genetic testing for pedigree accuracy is critical for managing genetic diversity in North American (NA) yak ( Bos grunniens), a population expanded mostly from imported zoological park specimens.  DNA testing also enhances species conservation by identifying recent B. taurus F1 hybrid ancestors (within three generations).  Biallelic single nucleotide polymorphisms (SNPs) can accomplish either task, but increases the marker count and costs necessary to achieve both.  Our aim was to identify novel, multifunctional, triallelic yak SNPs (tySNPs), with each having two alleles for yak parentage testing, and a third allele for identifying recent cattle introgression.  Methods:  Genome sequences were aligned to the cattle UMD3.1 assembly and SNPs were screened for 1) heterozygosity in a NA and a Chinese yak, 2) a third allele at high frequency in cattle, and 3) flanking sequences conserved in both species.  Subsequently, tySNPs were filtered for unique alignment to the haplotype-resolved F1 yak assembly.  Allele frequencies were estimated in a subset of 87 tySNPs by genotyping 170 NA yak. Results:  We identified 610 autosomal tySNPs, distributed in 441 clusters with 5 Mb average genome spacing.  The average NA yak minor allele frequency was high (0.296), while average introgressed cattle alleles were low (0.004).  In simulations with tySNPs, 28 were sufficient for globally-unique animal identification (P I=5.81x10 -12), 87 were able to exclude 19 random bulls from parentage at the 99% level without using the dam's genotype (P E=5.3x10 -4), and 87 were able to detect F1 hybridization events after three generations of yak backcrosses (1/16th B. taurus germplasm). Conclusions:  Identifying animals, determining parentage and detecting recent hybridization events was efficient with as few as 87 tySNPs.  A similar triallelic approach could be used with other bottlenecked Bos species that hybridize with cattle, such as NA plains bison ( B. bison).

Evaluation of EPAS1 variants for association with bovine congestive heart failure.

Background:  Bovine congestive heart failure (BCHF) has become increasingly prevalent in feedlot cattle in the Western Great Plains of North America. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. A protein variant of hypoxia-inducible factor 2 alpha (HIF2α, encoded by the endothelial PAS domain-containing protein 1 gene, EPAS1) was previously reported to be associated with pulmonary hypertension at altitudes exceeding 2,000 m. Our aim was to evaluate EPAS1 haplotypes for association with BCHF in feedlot cattle raised at moderate altitudes (1,200 m). Methods: Paired samples of clinical cases and unaffected controls were collected at four feedlots in Nebraska and Wyoming. Each pair (n =102) was matched for source, pen, breed type, sex, arrival date, and management conditions. Cases were identified by animal caretakers, euthanized, and diagnosis was confirmed at necropsy. Cases were derived from 30 different ranch operations, with the largest source contributing 32. Animals were tested for eight EPAS1 haplotypes encoding 36 possible different diploid combinations. Results: The common, ancestral EPAS1 haplotype encoding HIF2α with alanine (A) at position 606 and glycine (G) at position 610 was equally frequent in cases and controls (0.67). The EPAS1 variant haplotype reported to be associated with disease (encoding threonine (T) at position 606 and serine (S) at position 610) was not enriched in cases compared with controls (0.21 and 0.25, respectively). Frequencies of other EPAS1 haplotypes (e.g., encoding Q270, L362, or G671) were each less than 0.05 overall. McNemar's test with 45 discordant pairs showed the linked T606/S610 variant was not associated with BCHF (OR = 0.73, CI 95 0.38 -1.4, p-value = 0.37). Conclusions: HIF2α polypeptide variants were not significantly associated with BCHF in feedlot cattle at moderate altitudes. Thus, a wider search is needed to identify genetic risk factors underlying this disease.

A bovine CD18 signal peptide variant with increased binding activity to Mannheimia hemolytica leukotoxin.

Background:Mannheimia haemolytica is the major bacterial infectious agent of bovine respiratory disease complex and causes severe morbidity and mortality during lung infections. M. haemolytica secretes a protein leukotoxin (Lkt) that binds to the CD18 receptor on leukocytes, initiates lysis, induces inflammation, and causes acute fibrinous bronchopneumonia. Lkt binds the 22-amino acid CD18 signal peptide domain, which remains uncleaved in ruminant species. Our aim was to identify missense variation in the bovine CD18 signal peptide and measure the effects on Lkt binding. Methods: Missense variants in the integrin beta 2 gene ( ITGB2) encoding CD18 were identified by whole genome sequencing of 96 cattle from 19 breeds, and targeted Sanger sequencing of 1238 cattle from 46 breeds. The ability of different CD18 signal peptide variants to bind Lkt was evaluated by preincubating the toxin with synthetic peptides and applying the mixture to susceptible bovine cell cultures in cytotoxicity-blocking assays. Results: We identified 14 missense variants encoded on 15 predicted haplotypes, including a rare signal peptide variant with a cysteine at position 5 (C 5) instead of arginine (R 5). Preincubating Lkt with synthetic signal peptides with C 5 blocked cytotoxicity significantly better than those with R 5. The most potent synthetic peptide (C 5PQLLLLAGLLA) had 30-fold more binding activity compared to that with R 5. Conclusions: The results suggest that missense variants in the CD18 signal peptide affect Lkt binding, and animals carrying the C 5 allele may be more susceptible to the effects of Lkt. The results also identify a potent class of non-antibiotic Lkt inhibitors that could potentially protect cattle from cytotoxic effects during acute lung infections.

Using diverse U.S. beef cattle genomes to identify missense mutations in EPAS1, a gene associated with pulmonary hypertension.

The availability of whole genome sequence (WGS) data has made it possible to discover protein variants in silico. However, existing bovine WGS databases do not show data in a form conducive to protein variant analysis, and tend to under represent the breadth of genetic diversity in global beef cattle. Thus, our first aim was to use 96 beef sires, sharing minimal pedigree relationships, to create a searchable and publicly viewable set of mapped genomes relevant for 19 popular breeds of U.S. cattle. Our second aim was to identify protein variants encoded by the bovine endothelial PAS domain-containing protein 1 gene ( EPAS1), a gene associated with pulmonary hypertension in Angus cattle. The identity and quality of genomic sequences were verified by comparing WGS genotypes to those derived from other methods. The average read depth, genotype scoring rate, and genotype accuracy exceeded 14, 99%, and 99%, respectively. The 96 genomes were used to discover four amino acid variants encoded by EPAS1 (E270Q, P362L, A671G, and L701F) and confirm two variants previously associated with disease (A606T and G610S). The six EPAS1 missense mutations were verified with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry assays, and their frequencies were estimated in a separate collection of 1154 U.S. cattle representing 46 breeds. A rooted phylogenetic tree of eight polypeptide sequences provided a framework for evaluating the likely order of mutations and potential impact of EPAS1 alleles on the adaptive response to chronic hypoxia in U.S. cattle. This public, whole genome resource facilitates in silico identification of protein variants in diverse types of U.S. beef cattle, and provides a means of translating WGS data into a practical biological and evolutionary context for generating and testing hypotheses.

SNPs for parentage testing and traceability in globally diverse breeds of sheep.

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.

The Red clover necrotic mosaic virus origin of assembly is delimited to the RNA-2 trans-activator.

The bipartite RNA genome of Red clover necrotic mosaic virus (RCNMV) is encapsidated into icosahedral virions that exist as two populations: i) virions that co-package both genomic RNAs and ii) virions packaging multiple copies of RNA-2. To elucidate the packaging mechanism, we sought to identify the RCNMV origin of assembly sequence (OAS). RCNMV RNA-1 cannot package in the absence of RNA-2 suggesting that it does not contain an independent packaging signal. A 209 nt RNA-2 element expressed from the Tomato bushy stunt virus CP subgenomic promoter is co-assembled with genomic RNA-1 into virions. Deletion mutagenesis delimited the previously characterized 34 nt trans-activator (TA) as the minimal RCNMV OAS. From this study we hypothesize that RNA-1 must be base-paired with RNA-2 at the TA to initiate co-packaging. The addition of viral assembly illustrates the critical importance of the multifunctional TA element as a key regulatory switch in the RCNMV life cycle.

Controlled encapsidation of gold nanoparticles by a viral protein shell.

Icosahedral virus capsids demonstrate a high degree of selectivity in packaging cognate nucleic acid components during assembly. This packaging specificity, when integrated as part of a nanotechnological protocol, has the potential to encapsidate a wide array of foreign materials for delivery of therapeutics or biosensors into target cells. Red clover necrotic mosaic virus (RCNMV) exclusively packages two genomic ssRNAs initiated by a specific protein:RNA interaction between the RCNMV coat protein (CP) and the viral RNA origin of assembly (OAS) element. In the present work, an oligonucleotide mimic of the RCNMV OAS sequences is attached to Au nanoparticles as a recognition signal to initiate the virion-like assembly by RCNMV CP. Covalent linkage of the OAS to Au functions as a trigger for specific encapsidation and demonstrates that foreign cargo can be packaged into RCNMV virions.

The genomic RNA packaging scheme of Red clover necrotic mosaic virus.

Red clover necrotic mosaic virus (RCNMV) is a small icosahedral plant virus with a bipartite RNA genome. While the RCNMV genome consists of two RNAs, it has not been definitively established whether these RNAs are co-packaged into a single virion or packaged individually into separate virions. Biochemical evidence exists to support both hypotheses. To determine the genomic RNA complement within RCNMV, virions were subjected to heat treatments and UV crosslinking. A stable RNA-1:RNA-2 heterodimer was formed with both treatments establishing that RCNMV genomic RNAs are co-packaged into a single virion. Furthermore, RNA-2 homodimer and homotrimers were also observed indicating that some virions contain multiple copies of RNA-2 exclusively. These results indicate that RCNMV virions consist of two distinct populations: (i) virions containing both genomic RNAs; and (ii) virions with multiple copies of RNA-2. This type of hybrid packaging arrangement was unexpected and appears to be unique among the multipartite RNA viruses.