Proteome-wide analysis of Coxiella burnetii for conserved T-cell epitopes with presentation across multiple host species.

BACKGROUND: Coxiella burnetii is the Gram-negative bacterium responsible for Q fever in humans and coxiellosis in domesticated agricultural animals. Previous vaccination efforts with whole cell inactivated bacteria or surface isolated proteins confer protection but can produce a reactogenic immune responses. Thereby a protective vaccine that does not cause aberrant immune reactions is required. The critical role of T-cell immunity in control of C. burnetii has been made clear, since either CD8+ or CD4+ T cells can empower clearance. The purpose of this study was to identify C. burnetii proteins bearing epitopes that interact with major histocompatibility complexes (MHC) from multiple host species (human, mouse, and cattle). RESULTS: Of the annotated 1815 proteins from the Nine Mile Phase I (RSA 493) assembly, 402 proteins were removed from analysis due to a lack of inter-isolate conservation. An additional 391 proteins were eliminated from assessment to avoid potential autoimmune responses due to the presence of host homology. We analyzed the remaining 1022 proteins for their ability to produce peptides that bind MHCI or MHCII. MHCI and MHCII predicted epitopes were filtered and compared between species yielding 777 MHCI epitopes and 453 MHCII epitopes. These epitopes were further examined for presentation by both MHCI and MHCII, and for proteins that contained multiple epitopes. There were 31 epitopes that overlapped positionally between MHCI and MHCII across host species. Of these, there were 9 epitopes represented within proteins containing ≥ 5 total epitopes, where an additional 24 proteins were also epitope dense. In all, 55 proteins were found to contain high scoring T-cell epitopes. Besides the well-studied protein Com1, most identified proteins were novel when compared to previously studied vaccine candidates. CONCLUSION: These data represent the first proteome-wide evaluation of C. burnetii peptide epitopes. Furthermore, the inclusion of human, mouse, and bovine data capture a range of hosts for this zoonotic pathogen plus an important model organism. This work provides new vaccine targets for future vaccination efforts and enhances opportunities for selecting multiple T-cell epitope types to include within a vaccine.

GO-FAANG meeting: a Gathering On Functional Annotation of Animal Genomes.

The Functional Annotation of Animal Genomes (FAANG) Consortium recently held a Gathering On FAANG (GO-FAANG) Workshop in Washington, DC on October 7-8, 2015. This consortium is a grass-roots organization formed to advance the annotation of newly assembled genomes of domesticated and non-model organisms (www.faang.org). The workshop gathered together from around the world a group of 100+ genome scientists, administrators, representatives of funding agencies and commodity groups to discuss the latest advancements of the consortium, new perspectives, next steps and implementation plans. The workshop was streamed live and recorded, and all talks, along with speaker slide presentations, are available at www.faang.org. In this report, we describe the major activities and outcomes of this meeting. We also provide updates on ongoing efforts to implement discussions and decisions taken at GO-FAANG to guide future FAANG activities. In summary, reference datasets are being established under pilot projects; plans for tissue sets, morphological classification and methods of sample collection for different tissues were organized; and core assays and data and meta-data analysis standards were established.

Reduced lentivirus susceptibility in sheep with TMEM154 mutations.

Visna/Maedi, or ovine progressive pneumonia (OPP) as it is known in the United States, is an incurable slow-acting disease of sheep caused by persistent lentivirus infection. This disease affects multiple tissues, including those of the respiratory and central nervous systems. Our aim was to identify ovine genetic risk factors for lentivirus infection. Sixty-nine matched pairs of infected cases and uninfected controls were identified among 736 naturally exposed sheep older than five years of age. These pairs were used in a genome-wide association study with 50,614 markers. A single SNP was identified in the ovine transmembrane protein (TMEM154) that exceeded genome-wide significance (unadjusted p-value 3×10(-9)). Sanger sequencing of the ovine TMEM154 coding region identified six missense and two frameshift deletion mutations in the predicted signal peptide and extracellular domain. Two TMEM154 haplotypes encoding glutamate (E) at position 35 were associated with infection while a third haplotype with lysine (K) at position 35 was not. Haplotypes encoding full-length E35 isoforms were analyzed together as genetic risk factors in a multi-breed, matched case-control design, with 61 pairs of 4-year-old ewes. The odds of infection for ewes with one copy of a full-length TMEM154 E35 allele were 28 times greater than the odds for those without (p-value<0.0001, 95% CI 5-1,100). In a combined analysis of nine cohorts with 2,705 sheep from Nebraska, Idaho, and Iowa, the relative risk of infection was 2.85 times greater for sheep with a full-length TMEM154 E35 allele (p-value<0.0001, 95% CI 2.36-3.43). Although rare, some sheep were homozygous for TMEM154 deletion mutations and remained uninfected despite a lifetime of significant exposure. Together, these findings indicate that TMEM154 may play a central role in ovine lentivirus infection and removing sheep with the most susceptible genotypes may help eradicate OPP and protect flocks from reinfection.

A single codon insertion in PICALM is associated with development of familial subvalvular aortic stenosis in Newfoundland dogs.

Familial subvalvular aortic stenosis (SAS) is one of the most common congenital heart defects in dogs and is an inherited defect of Newfoundlands, golden retrievers and human children. Although SAS is known to be inherited, specific genes involved in Newfoundlands with SAS have not been defined. We hypothesized that SAS in Newfoundlands is inherited in an autosomal dominant pattern and caused by a single genetic variant. We studied 93 prospectively recruited Newfoundland dogs, and 180 control dogs of 30 breeds. By providing cardiac screening evaluations for Newfoundlands we conducted a pedigree evaluation, genome-wide association study and RNA sequence analysis to identify a proposed pattern of inheritance and genetic loci associated with the development of SAS. We identified a three-nucleotide exonic insertion in phosphatidylinositol-binding clathrin assembly protein (PICALM) that is associated with the development of SAS in Newfoundlands. Pedigree evaluation best supported an autosomal dominant pattern of inheritance and provided evidence that equivocally affected individuals may pass on SAS in their progeny. Immunohistochemistry demonstrated the presence of PICALM in the canine myocardium and area of the subvalvular ridge. Additionally, small molecule inhibition of clathrin-mediated endocytosis resulted in developmental abnormalities within the outflow tract (OFT) of Xenopus laevis embryos. The ability to test for presence of this PICALM insertion may impact dog-breeding decisions and facilitate reduction of SAS disease prevalence in Newfoundland dogs. Understanding the role of PICALM in OFT development may aid in future molecular and genetic investigations into other congenital heart defects of various species.

Effects of Low and High Maternal Protein Intake on Fetal Skeletal Muscle miRNAome in Sheep.

Prenatal maternal feeding plays an important role in fetal development and has the potential to induce long-lasting epigenetic modifications. MicroRNAs (miRNAs) are non-coding, single-stranded RNAs that serve as one epigenetic mechanism. Though miRNAs have crucial roles in fetal programming, growth, and development, there is limited data regarding the maternal diet and miRNA expression in sheep. Therefore, we analyzed high and low maternal dietary protein for miRNA expression in fetal longissimus dorsi. Pregnant ewes were fed an isoenergetic high-protein (HP, 160-270 g/day), low-protein (LP, 73-112 g/day), or standard-protein diet (SP, 119-198 g/day) during pregnancy. miRNA expression profiles were evaluated using the Affymetrix GeneChip miRNA 4.0 Array. Twelve up-regulated, differentially expressed miRNAs (DE miRNAs) were identified which are targeting 65 genes. The oar-3957-5p miRNA was highly up-regulated in the LP and SP compared to the HP. Previous transcriptome analysis identified that integrin and non-receptor protein tyrosine phosphatase genes targeted by miRNAs were detected in the current experiment. A total of 28 GO terms and 10 pathway-based gene sets were significantly (padj < 0.05) enriched in the target genes. Most genes targeted by the identified miRNAs are involved in immune and muscle disease pathways. Our study demonstrated that dietary protein intake during pregnancy affected fetal skeletal muscle epigenetics via miRNA expression.

Gene editing of pigs to control influenza A virus infections.

Proteolytic activation of the hemagglutinin (HA) glycoprotein by host cellular proteases is pivotal for influenza A virus (IAV) infectivity. Highly pathogenic avian influenza viruses possess the multibasic cleavage site of the HA which is cleaved by ubiquitous proteases, such as furin; in contrast, the monobasic HA motif is recognized and activated by trypsin-like proteases, such as the transmembrane serine protease 2 (TMPRSS2). Here, we aimed to determine the effects of TMPRSS2 on the replication of pandemic H1N1 and H3N2 subtype IAVs in the natural host, the pig. The use of the CRISPR/Cas 9 system led to the establishment of homozygous gene edited (GE) TMPRSS2 knockout (KO) pigs. Delayed IAV replication was demonstrated in primary respiratory cells of KO pigs in vitro. IAV infection in vivo resulted in significant reduction of virus shedding in the upper respiratory tract, and lower virus titers and pathological lesions in the lower respiratory tract of TMPRSS2 KO pigs as compared to WT pigs. Our findings could support the commercial use of GE pigs to minimize (i) the economic losses caused by IAV infection in pigs, and (ii) the emergence of novel IAVs with pandemic potential through genetic reassortment in the "mixing vessel", the pig.

Extent of linkage disequilibrium in large-breed dogs: chromosomal and breed variation.

The aim of this study was to better define the extent of linkage disequilibrium (LD) in populations of large-breed dogs and its variation by breed and chromosomal region. Understanding the extent of LD is a crucial component for successful utilization of genome-wide association studies and allows researchers to better define regions of interest and target candidate genes. Twenty-four Golden Retriever dogs, 28 Rottweiler dogs, and 24 Newfoundland dogs were genotyped for single-nucleotide polymorphism (SNP) data using a high-density SNP array. LD was calculated for all autosomes using Haploview. Decay of the squared correlation coefficient (r (2)) was plotted on a per-breed and per-chromosome basis as well as in a genome-wide fashion. The point of 50 % decay of r (2) was used to estimate the difference in extent of LD between breeds. Extent of LD was significantly shorter for Newfoundland dogs based upon 50 % decay of r (2) data at a mean of 344 kb compared to Golden Retriever and Rottweiler dogs at 715 and 834 kb, respectively (P < 0.0001). Notable differences in LD by chromosome were present within each breed and not strictly related to the length of the corresponding chromosome. Extent of LD is breed and chromosome dependent. To our knowledge, this is the first report of SNP-based LD for Newfoundland dogs, the first report based on genome-wide SNPs for Rottweilers, and an almost tenfold improvement in marker density over previous genome-wide studies of LD in Golden Retrievers.

Genome-Wide Scan of Wool Production Traits in Akkaraman Sheep.

The objective of this study was to uncover the genetic background of wool quality, a production trait, by estimating genomic heritability and implementing GWAS in Akkaraman sheep. The wool characteristics measured included fibre diameter (FD) and staple length (SL) at the age of 8 months and yearling fibre diameter (YFD), yearling staple length (YSL) and yearling greasy fleece weight (YGFW) at 18 months of age. Animals were genotyped using the Axiom 50 K Ovine Genotyping Array. Maximum likelihood estimations of a linear mixed model (LMM) were used to estimate genomic heritability, where GWAS was conducted following a score test of each trait. Genomic heritability estimates for the traits ranged between 0.22 and 0.63, indicating that phenotypes have a moderate range of heritability. One genome- and six chromosome-wide significant SNPs were associated with the wool traits in Akkaraman lambs. Accordingly, TRIM2, MND1, TLR2, RNF175, CEP290, TMTC3, RERE, SLC45A1, SOX2, MORN1, SKI, FAAP20, PRKCZ, GABRD, CFAP74, CALML6 and TMEM52 genes as well as nine uncharacterized regions (LOC101118971, LOC105609137, LOC105603067, LOC101122892, LOC106991694, LOC106991467, LOC106991455, LOC105616534 and LOC105609719) were defined as plausible candidates. The findings of this study shed light on the genetics of wool quality and yield for the Akkaraman breed and suggests targets for breeders during systematic breeding programmes.

Genome-wide association study of early liveweight traits in fat-tailed Akkaraman lambs.

Small ruminants, especially sheep, are essential for sustainable agricultural production systems, future food/nutrition security, and poverty reduction in developing countries. Within developed countries, the ability of sheep to survive on low-quality forage intake could act as buffer against climate change. Besides sheep's importance in sustainable agricultural production, there has been less ongoing work in terms of sheep genetics in Near East, Middle East and in Africa. For lamb meat production, body weight and average daily gain (ADG) until weaning are critical economic traits that affects the profitability of the industry. The current study aims to identify single nucleotide polymorphisms (SNPs) that are significantly associated with pre-weaning growth traits in fat tail Akkaraman lambs using a genome-wide association study (GWAS). A total of 196 Akkaraman lambs were selected for analysis. After quality control, a total of 31,936 SNPs and 146 lambs were used for subsequent analyses. PLINK 1.9 beta software was used for the analyses. Based on Bonferroni-adjusted p-values, one SNP (rs427117280) on chromosome 2 (OAR2) had significant associations with weaning weight at day 90 and ADG from day 0 to day 90, which jointly explains a 0.8% and 0.9% of total genetic variation respectively. The Ovis aries natriuretic peptide C (NPPC) could be considered as a candidate gene for the defined significant associations. The results of the current study will help to increase understanding of the variation in weaning weight and ADG until weaning of Akkaraman lambs and help enhance selection for lambs with improved weaning weight and ADG. However, further investigations are required for the identification of causal variants within the identified genomic regions.

Genome-wide association with footrot in hair and wool sheep.

Ovine footrot is an infectious disease with important contributions from Dichelobacter nodosus and Fusobacterium necrophorum. Footrot is characterized by separation of the hoof from underlying tissue, and this causes severe lameness that negatively impacts animal wellbeing, growth, and profitability. Large economic losses result from lost production as well as treatment costs, and improved genetic tools to address footrot are a valuable long-term goal. Prior genetic studies had examined European wool sheep, but hair sheep breeds such as Katahdin and Blackbelly have been reported to have increased resistance to footrot, as well as to intestinal parasites. Thus, footrot condition scores were collected from 251 U.S. sheep including Katahdin, Blackbelly, and European-influenced crossbred sheep with direct and imputed genotypes at OvineHD array (>500,000 single nucleotide polymorphism) density. Genome-wide association was performed with a mixed model accounting for farm and principal components derived from animal genotypes, as well as a random term for the genomic relationship matrix. We identified three genome-wide significant associations, including SNPs in or near GBP6 and TCHH. We also identified 33 additional associated SNPs with genome-wide suggestive evidence, including a cluster of 6 SNPs in a peak near the genome-wide significance threshold located near the glutamine transporter gene SLC38A1. These findings suggest genetic susceptibility to footrot may be influenced by genes involved in divergent biological processes such as immune responses, nutrient availability, and hoof growth and integrity. This is the first genome-wide study to investigate susceptibility to footrot by including hair sheep and also the first study of any kind to identify multiple genome-wide significant associations with ovine footrot. These results provide a foundation for developing genetic tests for marker-assisted selection to improve resistance to ovine footrot once additional steps like fine mapping and validation are complete.

A splice site mutation in a gene encoding for PDK4, a mitochondrial protein, is associated with the development of dilated cardiomyopathy in the Doberman pinscher.

Familial dilated cardiomyopathy is a primary myocardial disease that can result in the development of congestive heart failure and sudden cardiac death. Spontaneous animal models of familial dilated cardiomyopathy exist and the Doberman pinscher dog is one of the most commonly reported canine breeds. The objective of this study was to evaluate familial dilated cardiomyopathy in the Doberman pinscher dog using a genome-wide association study for a genetic alteration(s) associated with the development of this disease in this canine model. Genome-wide association analysis identified an area of statistical significance on canine chromosome 14 (p(raw) = 9.999e-05 corrected for genome-wide significance), fine-mapping of additional SNPs flanking this region localized a signal to 23,774,190-23,781,919 (p = 0.001) and DNA sequencing identified a 16-base pair deletion in the 5' donor splice site of intron 10 of the pyruvate dehydrogenase kinase 4 gene in affected dogs (p < 0.0001). Electron microscopy of myocardium from affected dogs demonstrated disorganization of the Z line, mild to moderate T tubule and sarcoplasmic reticulum dilation, marked pleomorphic mitochondrial alterations with megamitochondria, scattered mitochondria with whorling and vacuolization and mild aggregates of lipofuscin granules. In conclusion, we report the identification of a splice site deletion in the PDK4 gene that is associated with the development of familial dilated cardiomyopathy in the Doberman pinscher dog.

Estimates of genomic heritability and genome-wide association studies for blood parameters in Akkaraman sheep.

The aim of this study was to estimate genomic heritability and the impact that genetic backgrounds have on blood parameters in Akkaraman sheep by conducting genome-wide association studies and regional heritability mapping analysis. Genomic heritability estimates for blood parameters ranged from 0.00 to 0.55, indicating that measured phenotypes have a low to moderate heritability. A total of 7 genome- and 13 chromosome-wide significant SNPs were associated with phenotypic changes in 15 blood parameters tested. Accordingly, SCN7A, SCN9A, MYADM-like, CCDC67, ITGA9, MGAT5, SLC19A1, AMPH, NTRK2, MSRA, SLC35F3, SIRT6, CREB3L3, and NAV3 genes as well as three undefined regions (LOC101117887, LOC106991526 and LOC105608461) were suggested as candidates. Most of the identified genes were involved in basic biological processes that are essential to immune system function and cellular growth; specific functions include cellular transport, histone deacetylation, cell differentiation, erythropoiesis, and endocytosis. The top significant SNP for HCT, MCH, and MCHC was found within a genomic region mainly populated by the MYADM-like gene family. This region was previously suggested to be under historical selection pressure in many sheep breeds from various parts of the world. These results have implications on animal breeding program studies due to the effect that the genetic background has on blood parameters, which underlying many productive and wellness related traits.

Genetic Parameters and Genomic Regions Underlying Growth and Linear Type Traits in Akkaraman Sheep.

In the current study, the genetic architecture of growth and linear type traits were investigated in Akkaraman sheep. Estimations of genomic heritability, genetic correlations, and phenotypic correlations were implemented for 17 growth and linear type traits of 473 Akkaraman lambs by the univariate and multivariate analysis of animal mixed models. Correspondingly, moderate heritability estimates, as well as high and positive genetic/phenotypic correlations were found between growth and type traits. On the other hand, 2 genome-wide and 19 chromosome-wide significant single nucleotide polymorphisms were found to be associated with the traits as a result of animal mixed model-based genome-wide association analyses. Accordingly, we propose several genes located on different chromosomes (e.g., PRDM2, PTGDR, PTPRG, KCND2, ZNF260, CPE, GRID2, SCD5, SPIDR, ZNF407, HCN3, TMEM50A, FKBP1A, TLE4, SP1, SLC44A1, and MYOM3) as putative quantitative trait loci for the 22 growth and linear type traits studied. In our study, specific genes (e.g., TLE4, PTGDR, and SCD5) were found common between the traits studied, suggesting an interplay between the genetic backgrounds of these traits. The fact that four of the proposed genes (TLE4, MYOM3, SLC44A1, and TMEM50A) are located on sheep chromosome 2 confirms the importance of these genomic regions for growth and morphological structure in sheep. The results of our study are therefore of great importance for the development of efficient selection indices and marker-assisted selection programs, as well as for the understanding of the genetic architecture of growth and linear traits in sheep.

Genomic Analysis of Gastrointestinal Parasite Resistance in Akkaraman Sheep.

Genome-wide association studies (GWAS) have been used as an effective tool to understand the genetics of complex traits such as gastrointestinal parasite (GIP) resistance. The aim of this study was to understand the genetics of gastrointestinal parasite (nematodes, Moniezia spp., Eimeria spp.) resistance in Akkaraman sheep by performing genomic heritability estimations and conducting GWAS to uncover responsible genomic regions. This is one of the first studies to examine the genetic resistance of Akkaraman sheep to the tapeworm parasite. The samples from 475 animals were genotyped using the Axiom 50K Ovine Genotyping Array. Genomic heritability estimates ranged from 0.00 to 0.34 for parasite resistance traits. This indicates that measured phenotypes have low to moderate heritability estimates. A total of two genome-wide significant SNP associated with TNEM3 and ATRNL1 genes and 10 chromosome-wide significant SNPs related with 10 genes namely NELL1, ST6GALNAC3, HIPK1, SYT1, ALK, ZNF596, TMCO5A, PTH2R, LARGE1, and SCG2 were suggested as candidates for parasite resistance traits. The majority of these candidate genes were involved in several basic biological processes that are essential and important for immune system functions and cellular growth; specifically, inflammatory responses, cellular transport, cell apoptosis, cell differentiation, histone de-acetylation, and endocytosis. These results have implications for animal breeding program studies due to the effect that the genetic background has on parasite resistance, which underlies many productive, health, and wellness-related traits.

A high-density genome-wide association with absolute blood monocyte count in domestic sheep identifies novel loci.

Monocytes are a core component of the immune system that arise from bone marrow and differentiate into cells responsible for phagocytosis and antigen presentation. Their derivatives are often responsible for the initiation of the adaptive immune response. Monocytes and macrophages are central in both controlling and propagating infectious diseases such as infection by Coxiella burnetii and small ruminant lentivirus in sheep. Genotypes from 513 Rambouillet, Polypay, and Columbia sheep (Ovis aries) were generated using the Ovine SNP50 BeadChip. Of these sheep, 222 animals were subsequently genotyped with the Ovine Infinium® HD SNP BeadChip to increase SNP coverage. Data from the 222 HD genotyped sheep were combined with the data from an additional 258 unique sheep to form a 480-sheep reference panel; this panel was used to impute the low-density genotypes to the HD genotyping density. Then, a genome-wide association analysis was conducted to identify loci associated with absolute monocyte counts from blood. The analysis used a single-locus mixed linear model implementing EMMAX with age and ten principal components as fixed effects. Two genome-wide significant peaks (p < 5x10-7) were identified on chromosomes 9 and 1, and ten genome-wide suggestive peaks (p < 1x10-5) were identified on chromosomes 1, 2, 3, 4, 9, 10, 15, and 16. The identified loci were within or near genes including KCNK9, involved into cytokine production, LY6D, a member of a superfamily of genes, some of which subset monocyte lineages, and HMGN1, which encodes a chromatin regulator associated with myeloid cell differentiation. Further investigation of these loci is being conducted to understand their contributions to monocyte counts. Investigating the genetic basis of monocyte lineages and numbers may in turn provide information about pathogens of veterinary importance and elucidate fundamental immunology.

OvHV-2 Glycoprotein B Delivered by a Recombinant BoHV-4 Is Immunogenic and Induces Partial Protection against Sheep-Associated Malignant Catarrhal Fever in a Rabbit Model.

An efficacious vaccine for sheep-associated malignant catarrhal fever (SA-MCF) is important for the livestock industry. Research towards SA-MCF vaccine development is hindered by the absence of culture systems to propagate the causative agent, ovine herpesvirus-2 (OvHV-2), which means its genome cannot be experimentally modified to generate an attenuated vaccine strain. Alternative approaches for vaccine development are needed to deliver OvHV-2 antigens. Bovine herpesvirus 4 (BoHV-4) has been evaluated as a vaccine vector for several viral antigens with promising results. In this study, we genetically engineered BoHV-4 to express OvHV-2 glycoprotein B (gB) and evaluated its efficacy as an SA-MCF vaccine using a rabbit model. The construction of a viable recombinant virus (BoHV-4-AΔTK-OvHV-2-gB) and confirmation of OvHV-2 gB expression were performed in vitro. The immunization of rabbits with BoHV-4-AΔTK-OvHV-2-gB elicited strong humoral responses to OvHV-2 gB, including neutralizing antibodies. Following intra-nasal challenge with a lethal dose of OvHV-2, 42.9% of the OvHV-2 gB vaccinated rabbits were protected against SA-MCF, while all rabbits in the mock-vaccinated group succumbed to SA-MCF. Overall, OvHV-2 gB delivered by the recombinant BoHV-4 was immunogenic and partly protective against SA-MCF in rabbits. These are promising results towards an SA-MCF vaccine; however, improvements are needed to increase protection rates.

Natural genetic variation in Drosophila melanogaster reveals genes associated with Coxiella burnetii infection.

The gram-negative bacterium Coxiella burnetii is the causative agent of Query (Q) fever in humans and coxiellosis in livestock. Host genetics are associated with C. burnetii pathogenesis both in humans and animals; however, it remains unknown if specific genes are associated with severity of infection. We employed the Drosophila Genetics Reference Panel to perform a genome-wide association study to identify host genetic variants that affect host survival to C. burnetii infection. The genome-wide association study identified 64 unique variants (P < 10-5) associated with 25 candidate genes. We examined the role each candidate gene contributes to host survival during C. burnetii infection using flies carrying a null mutation or RNAi knockdown of each candidate. We validated 15 of the 25 candidate genes using at least one method. This is the first report establishing involvement of many of these genes or their homologs with C. burnetii susceptibility in any system. Among the validated genes, FER and tara play roles in the JAK/STAT, JNK, and decapentaplegic/TGF-ß signaling pathways which are components of known innate immune responses to C. burnetii infection. CG42673 and DIP-ε play roles in bacterial infection and synaptic signaling but have no previous association with C. burnetii pathogenesis. Furthermore, since the mammalian ortholog of CG13404 (PLGRKT) is an important regulator of macrophage function, CG13404 could play a role in host susceptibility to C. burnetii through hemocyte regulation. These insights provide a foundation for further investigation regarding the genetics of C. burnetii susceptibility across a wide variety of hosts.

Genes involved in immune, gene translation and chromatin organization pathways associated with Mycoplasma ovipneumoniae presence in nasal secretions of domestic sheep.

Mycoplasma ovipneumoniae contributes to polymicrobial pneumonia in domestic sheep. Elucidation of host genetic influences of M. ovipneumoniae nasal detection has the potential to reduce the incidence of polymicrobial pneumonia in sheep through implementation of selective breeding strategies. Nasal mucosal secretions were collected from 647 sheep from a large US sheep flock. Ewes of three breeds (Polypay n = 222, Rambouillet n = 321, and Suffolk n = 104) ranging in age from one to seven years, were sampled at three different times in the production cycle (February, April, and September/October) over four years (2015 to 2018). The presence and DNA copy number of M. ovipneumoniae was determined using a newly developed species-specific qPCR. Breed (P<0.001), age (P<0.024), sampling time (P<0.001), and year (P<0.001) of collection affected log10 transformed M. ovipneumoniae DNA copy number, where Rambouillet had the lowest (P<0.0001) compared with both Polypay and Suffolk demonstrating a possible genetic component to detection. Samples from yearlings, April, and 2018 had the highest (P<0.046) detected DNA copy number mean. Sheep genomic DNA was genotyped with the Illumina OvineHD BeadChip. Principal component analysis identified most of the variation in the dataset was associated with breed. Therefore, genome wide association analysis was conducted with a mixed model (EMMAX), with principal components 1 to 6 as fixed and a kinship matrix as random effects. Genome-wide significant (P<9x10-8) SNPs were identified on chromosomes 6 and 7 in the all-breed analysis. Individual breed analysis had genome-wide significant (P<9x10-8) SNPs on chromosomes 3, 4, 7, 9, 10, 15, 17, and 22. Annotated genes near these SNPs are part of immune (ANAPC7, CUL5, TMEM229B, PTPN13), gene translation (PIWIL4), and chromatin organization (KDM2B) pathways. Immune genes are expected to have increased expression when leukocytes encounter M. ovipneumoniae which would lead to chromatin reorganization. Work is underway to narrow the range of these associated regions to identify the underlying causal mutations.