Whole blood transcriptome analysis in dairy calves experimentally challenged with bovine herpesvirus 1 (BoHV-1) and comparison to a bovine respiratory syncytial virus (BRSV) challenge.

Bovine herpesvirus 1 (BoHV-1), is associated with several clinical syndromes in cattle, among which bovine respiratory disease (BRD) is of particular significance. Despite the importance of the disease, there is a lack of information on the molecular response to infection via experimental challenge with BoHV-1. The objective of this study was to investigate the whole-blood transcriptome of dairy calves experimentally challenged with BoHV-1. A secondary objective was to compare the gene expression results between two separate BRD pathogens using data from a similar challenge study with BRSV. Holstein-Friesian calves (mean age (SD) = 149.2 (23.8) days; mean weight (SD) = 174.6 (21.3) kg) were either administered BoHV-1 inoculate (1 × 107/mL × 8.5 mL) (n = 12) or were mock challenged with sterile phosphate buffered saline (n = 6). Clinical signs were recorded daily from day (d) -1 to d 6 (post-challenge), and whole blood was collected in Tempus RNA tubes on d six post-challenge for RNA-sequencing. There were 488 differentially expressed (DE) genes (p < 0.05, False Discovery rate (FDR) < 0.10, fold change ≥2) between the two treatments. Enriched KEGG pathways (p < 0.05, FDR <0.05); included Influenza A, Cytokine-cytokine receptor interaction and NOD-like receptor signalling. Significant gene ontology terms (p < 0.05, FDR <0.05) included defence response to virus and inflammatory response. Genes that are highly DE in key pathways are potential therapeutic targets for the treatment of BoHV-1 infection. A comparison to data from a similar study with BRSV identified both similarities and differences in the immune response to differing BRD pathogens.

Messenger RNA biomarkers of Bovine Respiratory Syncytial Virus infection in the whole blood of dairy calves.

Bovine Respiratory Syncytial Virus (BRSV) is a primary viral cause of Bovine Respiratory Disease (BRD) in young calves, which is responsible for substantial morbidity and mortality. Infection with BRSV induces global gene expression changes in respiratory tissues. If these changes are observed in tissues which are more accessible in live animals, such as whole blood, they may be used as biomarkers for diagnosis of the disease. Therefore, the objective of the current study was to elucidate the whole blood transcriptomic response of dairy calves to an experimental challenge with BRSV. Calves (Holstein-Friesian) were either administered BRSV inoculate (103.5 TCID50/ml × 15 ml) (n = 12) or sterile phosphate buffered saline (n = 6). Clinical signs were scored daily and whole blood was collected in Tempus RNA tubes immediately prior to euthanasia, at day 7 post-challenge. RNA was extracted from blood and sequenced (150 bp paired-end). The sequence reads were aligned to the bovine reference genome (UMD3.1) and EdgeR was subsequently employed for differential gene expression analysis. Multidimensional scaling showed that samples from BRSV challenged and control calves segregated based on whole blood gene expression changes, despite the BRSV challenged calves only displaying mild clinical symptoms of the disease. There were 281 differentially expressed (DE) genes (p < 0.05, FDR < 0.1, fold change > 2) between the BRSV challenged and control calves. The top enriched KEGG pathways and gene ontology terms were associated with viral infection and included "Influenza A", "defense response to virus", "regulation of viral life cycle" and "innate immune response". Highly DE genes involved in these pathways may be beneficial for the diagnosis of subclinical BRD from blood samples.

A mutation in the Warburg syndrome gene, RAB3GAP1, causes a similar syndrome with polyneuropathy and neuronal vacuolation in Black Russian Terrier dogs.

An autosomal recessive disease of Black Russian Terriers was previously described as a juvenile-onset, laryngeal paralysis and polyneuropathy similar to Charcot Marie Tooth disease in humans. We found that in addition to an axonal neuropathy, affected dogs exhibit microphthalmia, cataracts, and miotic pupils. On histopathology, affected dogs exhibit a spongiform encephalopathy characterized by accumulations of abnormal, membrane-bound vacuoles of various sizes in neuronal cell bodies, axons and adrenal cells. DNA from an individual dog with this polyneuropathy with ocular abnormalities and neuronal vacuolation (POANV) was used to generate a whole genome sequence which contained a homozygous RAB3GAP1:c.743delC mutation that was absent from 73 control canine whole genome sequences. An additional 12 Black Russian Terriers with POANV were RAB3GAP1:c.743delC homozygotes. DNA samples from 249 Black Russian Terriers with no known signs of POANV were either heterozygotes or homozygous for the reference allele. Mutations in human RAB3GAP1 cause Warburg micro syndrome (WARBM), a severe developmental disorder characterized by abnormalities of the eye, genitals and nervous system including a predominantly axonal peripheral neuropathy. RAB3GAP1 encodes the catalytic subunit of a GTPase activator protein and guanine exchange factor for Rab3 and Rab18 respectively. Rab proteins are involved in membrane trafficking in the endoplasmic reticulum, axonal transport, autophagy and synaptic transmission. The neuronal vacuolation and membranous inclusions and vacuoles in axons seen in this canine disorder likely reflect alterations of these processes. Thus, this canine disease could serve as a model for WARBM and provide insight into its pathogenesis and treatment.

Characterization of the rat developmental liver transcriptome.

Gene regulation and transcriptome studies have been enabled by the development of RNA-Seq applications for high-throughput sequencing platforms. Next generation sequencing is remarkably efficient and avoids many issues inherent in hybridization-based microarray methodologies including the exon-specific dependence of probe design. Biologically relevant transcripts including messenger and regulatory RNAs may now be quantified and annotated regardless of whether they have previously been observed. We used RNA-Seq to investigate global patterns of gene expression in early developing rat liver. Liver samples from timed-pregnant Lewis rats were collected at six fetal and neonatal stages [embryonic day (E)14, E16, E18, E20, postnatal day (P)1, P7], transcripts were sequenced using an Illumina HiSeq 2000, and data analysis was performed with the Tuxedo software suite. Genes and isoforms differing in abundance were queried for enrichment within functionally related gene groups using the Functional Annotation Tool of the DAVID Bioinformatics Database. While hematopoietic gene expression is initiated by E14, hepatocyte maturation is a gradual process involving clusters of genes responsible for response to nutrients and enzymes responsible for glycolysis and fatty acid catabolism. Following birth, a large cluster of differentially abundant genes was enriched for mitochondrial gene expression and cholesterol synthesis indicating that by 1 wk of age, the liver is engaged in lipid sensing and bile production. Clustering results for differentially abundant genes and isoforms were similar with the greatest difference for the E14/E16 comparison. Finally, a bioinformatic approach was used to annotate 1,307 novel liver transcripts assembled from sequences that aligned to intergenic regions of the rat genome.

Evolution of the bovine TLR gene family and member associations with Mycobacterium avium subspecies paratuberculosis infection.

Members of the Toll-like receptor (TLR) gene family occupy key roles in the mammalian innate immune system by functioning as sentries for the detection of invading pathogens, thereafter provoking host innate immune responses. We utilized a custom next-generation sequencing approach and allele-specific genotyping assays to detect and validate 280 biallelic variants across all 10 bovine TLR genes, including 71 nonsynonymous single nucleotide polymorphisms (SNPs) and one putative nonsense SNP. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and specialized beef and dairy breeds could not be differentiated despite an average polymorphism density of 1 marker/158 bp. Collectively, 160 tagSNPs and two tag insertion-deletion mutations (indels) were sufficient to predict 100% of the variation at 280 variable sites for both Bos subspecies and their hybrids, whereas 118 tagSNPs and 1 tagIndel predictively captured 100% of the variation at 235 variable sites for B. t. taurus. Polyphen and SIFT analyses of amino acid (AA) replacements encoded by bovine TLR SNPs indicated that up to 32% of the AA substitutions were expected to impact protein function. Classical and newly developed tests of diversity provide strong support for balancing selection operating on TLR3 and TLR8, and purifying selection acting on TLR10. An investigation of the persistence and continuity of linkage disequilibrium (r2≥0.50) between adjacent variable sites also supported the presence of selection acting on TLR3 and TLR8. A case-control study employing validated variants from bovine TLR genes recognizing bacterial ligands revealed six SNPs potentially eliciting small effects on susceptibility to Mycobacterium avium spp paratuberculosis infection in dairy cattle. The results of this study will broadly impact domestic cattle research by providing the necessary foundation to explore several avenues of bovine translational genomics, and the potential for marker-assisted vaccination.

A truncating mutation in ATP13A2 is responsible for adult-onset neuronal ceroid lipofuscinosis in Tibetan terriers.

A recessive, adult-onset neuronal ceroid-lipofuscinosis (NCL) occurs in Tibetan terriers. A genome-wide association study restricted this NCL locus to a 1.3Mb region of canine chromosome 2 which contains canine ATP13A2. NCL-affected dogs were homozygous for a single-base deletion in ATP13A2, predicted to produce a frameshift and premature termination codon. Homozygous truncating mutations in human ATP13A2 have been shown by others to cause Kufor-Rakeb syndrome (KRS), a rare neurodegenerative disease. These findings suggest that KRS is also an NCL, although analysis of KRS brain tissue will be needed to confirm this prediction. Generalized brain atrophy, behavioral changes, and cognitive decline occur in both people and dogs with ATP13A2 mutations; however, other clinical features differ between the species. For example, Tibetan terriers with NCL develop cerebellar ataxia not reported in KRS patients and KRS patients exhibit parkinsonism and pyramidal dysfunction not observed in affected Tibetan terriers. To see if ATP13A2 mutations could be responsible for some cases of human adult-onset NCL (Kufs disease), we resequenced ATP13A2 from 28 Kufs disease patients. None of these patients had ATP13A2 sequence variants likely to be causal for their disease, suggesting that mutations in this gene are not common causes of Kufs disease.

Diversity and evolution of 11 innate immune genes in Bos taurus taurus and Bos taurus indicus cattle.

The Toll-like receptor (TLR) and peptidoglycan recognition protein 1 (PGLYRP1) genes play key roles in the innate immune systems of mammals. While the TLRs recognize a variety of invading pathogens and induce innate immune responses, PGLYRP1 is directly microbicidal. We used custom allele-specific assays to genotype and validate 220 diallelic variants, including 54 nonsynonymous SNPs in 11 bovine innate immune genes (TLR1-TLR10, PGLYRP1) for 37 cattle breeds. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and we were unable to differentiate between the specialized B. t. taurus beef and dairy breeds, despite an average polymorphism density of one locus per 219 bp. Ninety-nine tagSNPs and one tag insertion-deletion polymorphism were sufficient to predict 100% of the variation at all 11 innate immune loci in both subspecies and their hybrids, whereas 58 tagSNPs captured 100% of the variation at 172 loci in B. t. taurus. PolyPhen and SIFT analyses of nonsynonymous SNPs encoding amino acid replacements indicated that the majority of these substitutions were benign, but up to 31% were expected to potentially impact protein function. Several diversity-based tests provided support for strong purifying selection acting on TLR10 in B. t. taurus cattle. These results will broadly impact efforts related to bovine translational genomics.

A neonatal encephalopathy with seizures in standard poodle dogs with a missense mutation in the canine ortholog of ATF2.

Neonatal encephalopathy with seizures (NEWS) is a previously undescribed autosomal recessive disease of standard poodle puppies. Affected puppies are small and weak at birth. Many die in their first week of life. Those surviving past 1 week develop ataxia, a whole-body tremor, and, by 4 to 6 weeks of age, severe generalized clonic-tonic seizures. None have survived to 7 weeks of age. Cerebella from affected puppies were reduced in size and often contained dysplastic foci consisting of clusters of intermixed granule and Purkinje neurons. We used deoxyribonucleic acid samples from related standard poodles to map the NEWS locus to a 2.87-Mb segment of CFA36, which contains the canine ortholog of ATF2. This gene encodes activating transcription factor 2 (ATF-2), which participates in the cellular responses to a wide variety of stimuli. We amplified and sequenced all coding regions of canine ATF2 from a NEWS-affected puppy and identified a T > G transversion that predicts a methionine-to-arginine missense mutation at amino acid position 51. Methionine-51 lies within a hydrophobic docking site for mitogen-activated protein kinases that activate ATF-2 so the arginine substitution is likely to interfere with ATF-2 activation. All 20 NEWS-affected puppies in the standard poodle family were homozygous for the mutant G allele. The 58 clinically normal family members were either G/T heterozygotes or homozygous for the ancestral T allele. There are no previous reports of spontaneous ATF2 mutations in people or animals; however, atf2-knockout mice have cerebellar lesions that are similar to those in puppies with NEWS.

Genetic mapping of canine multiple system degeneration and ectodermal dysplasia loci.

We characterized a movement disorder of Chinese Crested dogs clinically and pathologically indistinguishable from canine multiple system degeneration (CMSD) previously recognized in Kerry Blue Terriers. This fatal disease segregated as an autosomal recessive in a 51-dog pedigree of both breeds and their crosses. The occurrence of affected dogs among first-generation crosses demonstrated that the mutations causing multiple system degeneration in these breeds are allelic. The CMSD locus maps to CFA1 (LOD > 18) and haplotype analysis narrowed the CFA1 target region to a 15-Mb segment that contains orthologs of genes on HSA6, including PARK2, the gene for the ubiquitin ligase parkin. Mutations in human PARK2 cause the most common form of familial Parkinson's disease, autosomal recessive juvenile parkinsonism, which has clinical and pathological similarities to canine multiple system degeneration. A second phenotype, canine ectodermal dysplasia (CED), segregated in the pedigree as an autosomal dominant with homozygous lethality. Dogs with ectodermal dysplasia have a sparse hair coat and abnormal dentition that is characteristic of the "hairless" variety of Chinese Cresteds. CED mapped to a region of CFA17 (LOD > 14) containing orthologs from HSA2. EDAR, the gene for the ectodysplasin A1 receptor, occurs on HSA2 but was excluded as the cause of canine ectodermal dysplasia.