Influence of the geographic distribution of prion protein gene sequence variation on patterns of chronic wasting disease spread in white-tailed deer (Odocoileus virginianus)

ABSTRACT

Managing and controlling the spread of diseases in wild animal populations is challenging, especially for social and mobile species. Effective management benefits from information about disease susceptibility, allowing limited resources to be focused on areas or populations with a higher risk of infection. Chronic wasting disease (CWD), a transmissible spongiform encephalopathy that affects cervids, was detected in Colorado in the late 1960s. CWD was detected in Illinois and Wisconsin in 2002 and has since spread through many counties. Specific nucleotide variations in the prion protein gene (PRNP) sequence have been associated with reduced susceptibility to CWD in white-tailed deer. Though genetic resistance is incomplete, the frequency of deer possessing these mutations in a population is an important factor in disease spread (i.e. herd immunity). In this study we sequenced 625 bp of the PRNP gene from a sampling of 2433 deer from Illinois and Wisconsin. In north-central Illinois where CWD was first detected, counties had a low frequency of protective haplotypes (frequency <0.20); whereas in northwestern Illinois counties, where CWD cases have only more recently been detected, the frequency of protective haplotypes (frequency >0.30) was much higher (p < 0.05). Protective haplotype frequencies varied significantly among infected and uninfected geographic areas. The frequency of protective PRNP haplotypes may contribute to population level susceptibility and may shape the way CWD has spread through Illinois. Analysis of PRNP haplotype distribution could be a useful tool to assess CWD risk and allocate resources to contain and reduce the spread of infection.     

Prion protein gene sequence and chronic wasting disease susceptibility in white-tailed deer (Odocoileus virginianus)

ABSTRACT

The sequence of the prion protein gene (PRNP) affects susceptibility to spongiform encephalopathies, or prion diseases in many species. In white-tailed deer, both coding and non-coding single nucleotide polymorphisms have been identified in this gene that correlate to chronic wasting disease (CWD) susceptibility. Previous studies examined individual nucleotide or amino acid mutations; here we examine all nucleotide polymorphisms and their combined effects on CWD. A 626 bp region of PRNP was examined from 703 free-ranging white-tailed deer. Deer were sampled between 2002 and 2010 by hunter harvest or government culling in Illinois and Wisconsin. Fourteen variable nucleotide positions were identified (4 new and 10 previously reported). We identified 68 diplotypes comprised of 24 predicted haplotypes, with the most common diplotype occurring in 123 individuals. Diplotypes that were found exclusively among positive or negative animals were rare, each occurring in less than 1% of the deer studied. Only one haplotype (C, odds ratio 0.240) and 2 diplotypes (AC and BC, odds ratios of 0.161 and 0.108 respectively) has significant associations with CWD resistance. Each contains mutations (one synonymous nucleotide 555C/T and one nonsynonymous nucleotide 286G/A) at positions reported to be significantly associated with reduced CWD susceptibility. Results suggest that deer populations with higher frequencies of haplotype C or diplotypes AC and BC might have a reduced risk for CWD infection – while populations with lower frequencies may have higher risk for infection. Understanding the genetic basis of CWD has improved our ability to assess herd susceptibility and direct management efforts within CWD infected areas.     

Geographically widespread honeybee‐gut symbiont subgroups show locally distinct antibiotic‐resistant patterns

How long‐term antibiotic treatment affects host bacterial associations is still largely unknown. The honeybee‐gut microbiota has a simple composition, so we used this gut community to investigate how long‐term antibiotic treatment affects host‐associated microbiota. We investigated the phylogenetic relatedness, genomic content (GC percentage, genome size, number of genes and CRISPR) and antibiotic‐resistant genes (ARG) for strains from two abundant members of the honeybee core gut microbiota (Gilliamella apicola and Snodgrassella alvi). Domesticated honeybees are subjected to geographically different management policies, so we used two research apiaries, representing different antibiotic treatment regimens in their apiculture: low antibiotic usage (Norway) and high antibiotic usage (Arizona, USA). We applied whole‐genome shotgun sequencing on 48 G. apicola and 22 S. alvi. We identified three predominating subgroups of G. apicola in honeybees from both Norway and Arizona. For G. apicola, genetic content substantially varied between subgroups and distance similarity calculations showed similarity discrepancy between subgroups. Functional differences between subgroups, such as pectin‐degrading enzymes (G. apicola), were also identified. In addition, we identified horizontal gene transfer (HGT) of transposon (Tn10)‐associated tetracycline resistance (Tet B) across the G. apicola subgroups in the Arizonan honeybees, using interspace polymorphisms in the Tet B determinant. Our results support that honeybee‐gut symbiont subgroups can resist long‐term antibiotic treatment and maintain functionality through acquisition of geographically distinct antibiotic‐resistant genes by HGT.     

Gene expression alterations in Rocky Mountain elk infected with chronic wasting disease

Chronic wasting disease (CWD) is an invariably fatal neurologic disease that naturally infects mule deer, white tailed deer and elk. The understanding of CWD neurodegeneration at a molecular level is very limited. In this study, microarray analysis was performed to determine changes in the gene expression profiles in six different tissues including brain, midbrain, thalamus, spleen, RPLN and tonsil of CWD-infected elk in comparison to non-infected healthy elk, using 24,000 bovine specific oligo probes. In total, 329 genes were found to be differentially expressed (> 2.0-fold) between CWD negative and positive brain tissues, with 132 genes upregulated and 197 genes downregulated. There were 249 DE genes in the spleen (168 up- and 81 downregulated), 30 DE genes in the retropharyngeal lymph node (RPLN) (18 up- and 12 downregulated), and 55 DE genes in the tonsil (21 up- and 34 downregulated). Using Gene Ontology (GO), the DE genes were assigned to functional groups associated with cellular process, biological regulation, metabolic process, and regulation of biological process. For all brain tissues, the highest ranking networks for DE genes identified by Ingenuity Pathway Analysis (IPA) were associated with neurological disease, cell morphology, cellular assembly and organization. Quantitative real-time PCR (qRT-PCR) validated the expression of DE genes primarily involved in different regulatory pathways, including neuronal signaling and synapse function, calcium signaling, apoptosis and cell death and immune cell trafficking and inflammatory response. This is the first study to evaluate altered gene expression in multiple organs including brain from orally infected elk and the results will improve our understanding of CWD neurodegeneration at the molecular level.     

Gene expression alterations in Rocky Mountain elk infected with chronic wasting disease

Chronic wasting disease (CWD) is an invariably fatal neurologic disease that naturally infects mule deer, white tailed deer and elk. The understanding of CWD neurodegeneration at a molecular level is very limited. In this study, microarray analysis was performed to determine changes in the gene expression profiles in six different tissues including brain, midbrain, thalamus, spleen, RPLN and tonsil of CWD-infected elk in comparison to non-infected healthy elk, using 24,000 bovine specific oligo probes. In total, 329 genes were found to be differentially expressed (> 2.0-fold) between CWD negative and positive brain tissues, with 132 genes upregulated and 197 genes downregulated. There were 249 DE genes in the spleen (168 up- and 81 downregulated), 30 DE genes in the retropharyngeal lymph node (RPLN) (18 up- and 12 downregulated), and 55 DE genes in the tonsil (21 up- and 34 downregulated). Using Gene Ontology (GO), the DE genes were assigned to functional groups associated with cellular process, biological regulation, metabolic process, and regulation of biological process. For all brain tissues, the highest ranking networks for DE genes identified by Ingenuity Pathway Analysis (IPA) were associated with neurological disease, cell morphology, cellular assembly and organization. Quantitative real-time PCR (qRT-PCR) validated the expression of DE genes primarily involved in different regulatory pathways, including neuronal signaling and synapse function, calcium signaling, apoptosis and cell death and immune cell trafficking and inflammatory response. This is the first study to evaluate altered gene expression in multiple organs including brain from orally infected elk and the results will improve our understanding of CWD neurodegeneration at the molecular level.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.