Genomic regions associated with Mycoplasma ovipneumoniae presence in nasal secretions of domestic sheep.

Mycoplasma ovipneumoniae contributes to polymicrobial pneumonia in domestic sheep. Elucidation of host genetic components to M. ovipneumoniae nasal detection would have the potential to reduce the incidence of pneumonia. 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 of the production cycle (February, April, and September/October) over four years (2015 to 2018). The presence and quantity of M. ovipneumoniae was determined using a species-specific qPCR with 3 to 10 sampling times per sheep. 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 principle 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.