Deletion variant near ZNF389 is associated with control of ovine lentivirus in multiple sheep flocks

A high-density genome-wide approach reveals novel genetic markers linked to small ruminant lentivirus susceptibility in sheep

Visna/Maedi virus (VMV) is lentiviral disease of sheep responsible for severe production losses. Multiple genomic regions associated with infection were reported indicating genetic complexity. In this study, a combined genome-wide approach using a high-density SNP array has been performed, comparing VMV-infected (n = 78) and non-infected (n = 66) individuals of the Valle del Belice breed. The serological tests showed a seroprevalence of 26%. The comparison among results from different approaches (GWAS, Fisher's exact test and the FST analysis) revealed two association signals: on OAR03 close to the GRIN2B gene and on OAR05 close to the TMEM232 gene. To the best of our knowledge, there has been no previous association between these genes and lentiviral infection in any species. The GRIN2B gene plays a role in pain response, synaptic transmission, and receptor clustering, while TMEM232 is involved in the development of immune-related disorders. The results highlighted new aspects of the genetic complexity related to the resistance/susceptibility to VMV in sheep, confirming that studies on different breeds can lead to different results. The ideal approach for validation of the markers identified in our study is to use samples from a population independent from the discovery population with the same phenotype used in the discovery stage.

The genetic variability of small-ruminant lentiviruses and its impact on tropism, the development of diagnostic tests and vaccines and the effectiveness of control programmes

Introduction

Maedi-visna virus and caprine arthritis encephalitis virus are two closely related lentiviruses which cause multisystemic, progressive and persistent infection in goats and sheep. Because these viruses frequently cross the species barrier, they are considered to be one genetic group called small-ruminant lentiviruses (SRLV). They have in vivo tropism mainly for monocytes and macrophages and organ tropism with unknown mechanisms. Typical clinical signs are pneumonia in sheep, arthritis in goats, and mastitis in both species. Infection with SRLV cannot currently be treated or prevented, and control programmes are the only approaches to avoiding its spread. These programmes rely mainly on annual serological testing and elimination of positive animals. However, the high genetic and antigenic variability of SRLV complicate their early and definitive diagnosis. The objective of this review is to summarise the current knowledge of SRLV genetic variation and its implications for tropism, the development of diagnostic tests and vaccines and the effectiveness of control and eradication programmes.

Material and Methods

Subject literature was selected from the PubMed and the Google Scholar databases.

Results

The high genetic diversity of SRLV affects the performance of diagnostic tools and therefore control programmes. For the early and definitive diagnosis of SRLV infection, a combination of serological and molecular tests is suggested. Testing by PCR can also be considered for sub-yearling animals. There are still significant gaps in our knowledge of the epidemiology, immunology and biology of SRLV and their impact on animal production and welfare.

Conclusion

This information may aid selection of the most effective SRLV spread reduction measures.

Field demonstration analyzing the implementation of individual animal electronic identification and genetic testing in western range sheep flocks

Adoption of electronic identification ear tags (EID) and DNA testing by commercial range sheep producers in the Western United States has been low, despite the availability of these technologies for over a decade. Jointly, these technologies offer an approach to provide individual animal performance data to improve flock health, genetic and reproductive management. This project involved a collaboration with five California sheep producers representing a broad geographic range, varying levels of pre-project EID adoption, and diverse operational practices. Tissue samples were collected from, and ear EIDs were placed in, a total of 2,936 rams and their potential lambs. We partnered with a commercial packing company, Superior Farms, to genotype the animals. Superior Farms used a targeted genotyping panel to assign parentage, and link individual animal identification (ID) to camera-graded carcass measurements. This enabled the collection of individual progeny carcass data and provided insight into sire performance, providing for the within-flock identification of prolific sires that were producing lambs with significantly more saleable meat as compared to their flock mates. Overall, almost 91% of lambs were successfully matched to their sire, and prolificacy ranging from 0-135 lambs per ram. There was as much as an $80 difference in the average edible product from camera-graded carcasses derived from lamb groups sired by different rams. A partial budget analysis modeling investment in an EID system coupled with an autodrafter and scale to collect individual weights and improve labor efficiency during processing, and a sheep flip chute to improve worker safety during foot trimmings, yielded a greater than 7:1 return on investment over a five-year time frame. Ideally, the data collection enabled by EIDs and DNA testing would feed into data-driven genetic evaluation programs to enable selection for more productive and profitable animals, and allow the US sheep industry to accelerate the rate of genetic improvement.

Meta-analysis of heritability estimates and genome-wide association for tick-borne haemoparasites in African cattle

The control of tick-borne haemoparasites in cattle largely relies on the use of acaricide drugs against the tick vectors, with some vaccination also being used against selected pathogens. These interventions can be difficult in Africa, where accessibility and cost of vaccines can be issues, and the increasing resistance of tick vectors to the widely used acaricides is a complication to disease control. A potential complementary control strategy could be the exploitation of any natural host genetic resistance to the pathogens. However, there are currently very few estimates of the extent of host resistance to tick-borne haemoparasites, and a significant contributing factor to this knowledge gap is likely to be the difficulty of collecting appropriate samples and data in the smallholder systems that predominate livestock production in low- and middle-income countries, particularly at scale. In this study, we have estimated the heritability for the presence/absence of several important haemoparasite species (including Anaplasma marginale, Babesia bigemina, Babesia bovis, and Ehrlichia ruminantium), as well as for relevant traits such as body weight and body condition score (BCS), in 1,694 cattle from four African countries (Burkina Faso, Ghana, Nigeria, and Tanzania). Heritability estimates within countries were mostly not significant, ranging from 0.05 to 0.84 across traits and countries, with standard errors between 0.07 and 0.91. However, the weighted mean of heritability estimates was moderate and significant for body weight and BCS (0.40 and 0.49, respectively), with significant heritabilities also observed for the presence of A. marginale (0.16) and E. ruminantium (0.19). In a meta-analysis of genome-wide association studies (GWAS) for these traits, two peaks were identified as reaching the suggestive significance threshold (p < 1.91 × 10-7 and p < 1.89 × 10-7, respectively): one on chromosome 24 for BCS and one on chromosome 8 for the E. ruminantium infection status. These findings indicate that there is likely to be a genetic basis that contributes to pathogen presence/absence for tick-borne haemoparasite species, which could potentially be exploited to improve cattle resistance in Africa to the economically important diseases caused by these pathogens.

Lentivirus Susceptibility in Brazilian and US Sheep with TMEM154 Mutations

Small ruminant lentiviruses (SRLVs) affect sheep and goats worldwide. The major gene related to SRLV infections is the Transmembrane Protein Gene 154 (TMEM154). We estimated the haplotype frequencies of TMEM154 in the USA (USDA-ARS) and Brazil (Embrapa) Gene Banks by using two different SNP genotyping methodologies, Fluidigm^TM and KASP^TM. We also genotyped the ZNF389_ss748775100 deletion variant in Brazilian flocks. A total of 1040 blood samples and 112 semen samples from 15 Brazilian breeds were genotyped with Fluidigm for the SNP ZNF389_ss748775100 and 12 TMEM154 SNPs. A total of 484 blood samples from the Santa Inês breed and 188 semen samples from 14 North American sheep breeds were genotyped with KASP for 6 TMEM154 SNPs. All the Brazilian samples had the "I/I" genotype for the ZNF389_ss748775100 mutation. There were 25 TMEM154 haplotypes distributed across the Brazilian breeds, and 4 haplotypes in the US breeds. Haplotypes associated with susceptibility were present in almost all breeds, which suggests that genetic testing can help to improve herd health and productivity by selecting non-susceptible animals as founders of the next generations. Fluidigm and KASP are reliable assays when compared with Beadchip arrays. Further studies are necessary to understand the unknown role of TMEM154 mutations, host-pathogen interaction and new genes associated with the clinical condition.

A locus conferring tolerance to Theileria infection in African cattle

East Coast fever, a tick-borne cattle disease caused by the Theileria parva parasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Bos indicus (Boran) cattle demonstrating heritable tolerance to infection with T. parva (h² = 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region on bovine chromosome 15 that is significantly associated with survival outcome following T. parva exposure. Testing this locus in an independent cohort of animals replicates this association with survival following T. parva infection. A stop gained variant in a paralogue of the FAF1 gene in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa’s most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection by T. parva.

More than a million cattle die of East Coast fever in Africa each year, the impact of which disproportionately falls onto low-income, smallholder farmers. The lack of a widely accessible vaccine, heavy reliance on chemicals to control the tick vector and inadequate drug treatments means that new approaches for controlling the disease are urgently required. Through a genetic study of an extended pedigree of Boran cattle that are more than three times less likely to succumb to the disease than matched controls, we identify a region on chromosome 15 of the cattle genome associated with a high level of tolerance to the disease. We show that a nonsense variant in a predicted paralogue of FAS-associated factor 1 (FAF1) in this region is also associated with survival in an independent cohort, and is linked to rates of cell expansion during infection. This genetic variant can therefore support marker-assisted selection, allowing farmers to breed tolerant cattle and offers a route to introduce this beneficial DNA to non-native breeds, enabling reduced disease incidence and increased productivity, which would be of benefit to millions of rural smallholder farmers across Africa.

Transcriptome Analysis for Genes Associated with Small Ruminant Lentiviruses Infection in Goats of Carpathian Breed

Small ruminant lentiviruses (SRLV) are economically important viral pathogens of sheep and goats. SRLV infection may interfere in the innate and adaptive immunity of the host, and genes associated with resistance or susceptibility to infection with SRLV have not been fully recognized. The presence of animals with relatively high and low proviral load suggests that some host factors are involved in the control of virus replication. To better understand the role of the genes involved in the host response to SRLV infection, RNA sequencing (RNA-seq) method was used to compare whole gene expression profiles in goats carrying both a high (HPL) and low (LPL) proviral load of SRLV and uninfected animals. Data enabled the identification of 1130 significant differentially expressed genes (DEGs) between control and LPL groups: 411 between control and HPL groups and 1434 DEGs between HPL and LPL groups. DEGs detected between the control group and groups with a proviral load were found to be significantly enriched in several gene ontology (GO) terms, including an integral component of membrane, extracellular region, response to growth factor, inflammatory and innate immune response, transmembrane signaling receptor activity, myeloid differentiation primary response gene 88 (MyD88)-dependent toll-like receptor signaling pathway as well as regulation of cytokine secretion. Our results also demonstrated significant deregulation of selected pathways in response to viral infection. The presence of SRLV proviral load in blood resulted in the modification of gene expression belonging to the toll-like receptor signaling pathway, the tumor necrosis factor (TNF) signaling pathway, the cytokine-cytokine receptor interaction, the phagosome, the Ras signaling pathway, the phosphatidylinositol 3-kinase (PI3K)/protein kinase B (AKT) (PI3K-Akt) signaling pathway and rheumatoid arthritis. It is worth mentioning that the most predominant in all pathways were genes represented by toll-like receptors, tubulins, growth factors as well as interferon gamma receptors. DEGs detected between LPL and HPL groups were found to have significantly enriched regulation of signaling receptor activity, the response to toxic substances, nicotinamide adenine dinucleotide (NADH) dehydrogenase complex assembly, cytokine production, vesicle, and vacuole organization. In turn, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway tool classified DEGs that enrich molecular processes such as B and T-cell receptor signaling pathways, natural killer cell-mediated cytotoxicity, Fc gamma R-mediated phagocytosis, toll-like receptor signaling pathways, TNF, mammalian target of rapamycin (mTOR) signaling and forkhead box O (Foxo) signaling pathways, etc. Our data indicate that changes in SRLV proviral load induced altered expression of genes related to different biological processes such as immune response, inflammation, cell locomotion, and cytokine production. These findings provide significant insights into defense mechanisms against SRLV infection. Furthermore, these data can be useful to develop strategies against SRLV infection by selection of animals with reduced SRLV proviral concentration that may lead to a reduction in the spread of the virus.

A novel 2 bp deletion variant in Ovine-DRB1 gene is associated with increased Visna/maedi susceptibility in Turkish sheep

Visna/maedi (VM) is a multisystemic lentivirus infection of sheep that affecting sheep industry across the globe. TMEM154 gene has been identified to be a major VM-associated host gene, nevertheless, a recent study showed that the frequency of the VM-resistant TMEM154 haplotypes was very low or absent in indigenous sheep. Thus, the present study was designed to determine other possible co-receptors associated with VM. For this purpose, DRB1 gene, which is renowned for its role in host immune response against various diseases was targeted. A total number of 151 case–control matched pairs were constructed from 2266 serologically tested sheep. A broad range of DRB1 haplotype diversity was detected by sequence-based genotyping. Moreover, a novel 2 bp deletion (del) in the DRB1 intron 1 was identified. For the final statistic, the sheep carrying VM-resistant TMEM154 diplotypes were removed and a McNemar’s test with a matched pairs experimental design was conducted. Consequently, it was identified for the first time that the 2 bp del variant is a genetic risk factor for VM (p value 0.002; chi-square 8.31; odds ratio 2.9; statistical power 0.90) in the dominant model. Thus, negative selection for 2 bp del variant could decrease VM infection risk in Turkish sheep.

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 log₁₀ 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.

A DNA Regulatory Element Haplotype at Zinc Finger Genes Is Associated with Host Resilience to Small Ruminant Lentivirus in Two Sheep Populations

Small ruminant lentivirus (SRLV) causes Maedi-Visna or Ovine Progressive Pneumonia in sheep and creates insidious livestock production losses. This retrovirus is closely related to human immunodeficiency virus and currently has no vaccines or cure. Genetic marker assisted selection for sheep disease resiliency presents an attractive management solution. Previously, we identified a region containing a cluster of zinc finger genes that had association with ovine SRLV proviral concentration. Trait-association analysis validated a small insertion/deletion variant near ZNF389 (rs397514112) in multiple sheep breeds. In the current study, 543 sheep from two distinct populations were genotyped at 34 additional variants for fine mapping of the regulatory elements within this locus. Variants were selected based on ChIP-seq annotation data from sheep alveolar macrophages that defined active cis-regulatory elements predicted to influence zinc finger gene expression. We present a haplotype block of variants within regulatory elements that have improved associations and larger effect sizes (up to 4.7-fold genotypic difference in proviral concentration) than the previously validated ZNF389 deletion marker. Hypotheses for the underlying causal mutation or mutations are presented based on changes to in silico transcription factor binding sites. These variants offer alternative markers for selective breeding and are targets for future functional mutation assays.

Achievements of an eradication programme against caprine arthritis encephalitis virus in South Tyrol, Italy

Small ruminant lentivirus infections in goats affect both production and animal welfare. This represents a threat to the qualitative and quantitative growth of goat farming, recently observed in mountainous regions such as the Autonomous Province of Bolzano – South Tyrol (Italy). To monitor and eradicate the caprine arthritis encephalitis virus in this goat population, a compulsory eradication campaign was launched, based on a strict census of small ruminants and yearly serological testing of all animals, followed by the consequent culling of seropositive individuals. The campaign succeeded in completely eliminating cases of clinical disease in goats, while drastically reducing the seroprevalence at the herd as well as individual animal level. The serological outcome of the introduced control measures was determined using commercially available ELISA kits, demonstrating their suitability for use in this type of campaign, aimed at reducing seroprevalence as well as clinical manifestations of these infections. However, this clear success is diminished by the failure to achieve a complete eradication of these viruses. The reasons leading to the observed tailing phenomenon and the occurrence of new infections in already sanitised flocks are discussed and implementation of further measures are proposed.

Generation of a molecular clone of an attenuated lentivirus, a first step in understanding cytopathogenicity and virulence

Small ruminant lentiviruses infect goats and sheep, inducing clinical disease in a minority of infected animals. Following an eradication campaign, clinical cases may disappear in a population. The complete elimination of these lentiviruses is however difficult to achieve and the spreading of less virulent strains often parallels the elimination of their virulent counterparts. Here, we characterized three such strains isolated from a flock in the post-eradication phase. We completely sequenced their genomes, showing that one of the isolates was most probably the product of a recombination event between the other two viruses. By comparing the sequences of these isolates with those of virulent strains, we found evidence that particular LTR mutations may explain their attenuated phenotype. Finally, we constructed an infectious molecular clone representative of these viruses, analyzing its replication characteristics in different target cells. This clone will permit us to explore the molecular correlates of cytopathogenicity and virulence.

Genome-Wide Association Identifies SLC2A9 and NLN Gene Regions as Associated with Entropion in Domestic Sheep

Entropion is an inward rolling of the eyelid allowing contact between the eyelashes and cornea that may lead to blindness if not corrected. Although many mammalian species, including humans and dogs, are afflicted by congenital entropion, no specific genes or gene regions related to development of entropion have been reported in any mammalian species to date. Entropion in domestic sheep is known to have a genetic component therefore, we used domestic sheep as a model system to identify genomic regions containing genes associated with entropion. A genome-wide association was conducted with congenital entropion in 998 Columbia, Polypay, and Rambouillet sheep genotyped with 50,000 SNP markers. Prevalence of entropion was 6.01%, with all breeds represented. Logistic regression was performed in PLINK with additive allelic, recessive, dominant, and genotypic inheritance models. Two genome-wide significant (empirical P<0.05) SNP were identified, specifically markers in SLC2A9 (empirical P = 0.007; genotypic model) and near NLN (empirical P = 0.026; dominance model). Six additional genome-wide suggestive SNP (nominal P<1x10^-5) were identified including markers in or near PIK3CB (P = 2.22x10^-6; additive model), KCNB1 (P = 2.93x10^-6; dominance model), ZC3H12C (P = 3.25x10^-6; genotypic model), JPH1 (P = 4.68x20^-6; genotypic model), and MYO3B (P = 5.74x10^-6; recessive model). This is the first report of specific gene regions associated with congenital entropion in any mammalian species, to our knowledge. Further, none of these genes have previously been associated with any eyelid traits. These results represent the first genome-wide analysis of gene regions associated with entropion and provide target regions for the development of sheep genetic markers for marker-assisted selection.