Genome-wide association study of footrot in Texel sheep

Genome-wide association studies of parasite resistance, productivity and immunology traits in Scottish Blackface sheep

Gastrointestinal parasitism represents a global problem for grazing ruminants, which can be addressed sustainably by breeding animals to be more resistant against infection by parasites. The aim of this study was to assess the genetic architecture underlying traits associated with gastrointestinal parasite resistance, immunological profile and production in meat sheep, and identify and characterise candidate genes affecting these traits. Data on gastrointestinal parasite infection (faecal egg counts for Strongyles (FECS) and Nematodirus (FECN) and faecal oocyst counts for Coccidia, along with faecal soiling scores (DAG), characterised by the accumulation of faeces around the perineum) and production (live weight (LWT)) were gathered from a flock Scottish Blackface lambs at three and four months of age. Data on the immune profile were also collected from a subset of these lambs at two and five months of age. Immune traits included the production of Interferon-γ (IFN-γ), Interleukin (IL)-4 and IL-10 following stimulation of whole blood with pokeweed mitogen (PWM) or antigen from the gastric parasite Teladorsagia circumcincta (T-ci), and serum levels of T. circumcincta-specific immunoglobulin A (IgA). Animals were genotyped with genome-wide DNA arrays, and a total of 1 766 animals and 45 827 Single Nucleotide Polymorphisms (SNPs) were retained following quality control and imputation. Genome-wide association studies were performed for 24 traits. The effects of individual markers with significant effects were estimated, and the genotypic effect solutions were used to estimate additive and dominance effects, and the proportion of additive genetic variance attributed to each SNP locus. A total of 15 SNPs were associated at least at a suggestive level with FECS, FECN, DAG, IgA, PWM-induced IFN-γ and IL-4, and T-ci-induced IL-10. This study uncovered 52 genes closely related to immune function in proximity to these SNPs. A number of genes encoding C-type lectins and killer cell lectin-like family members were close to a SNP associated with FECN, while several genes encoding IL-1 cytokine family members were found to be associated with IgA. Potential candidate genes belonging to or in close proximity with the Major Histocompatibility Complex (MHC) were revealed, including Homeostatic Iron Regulator and butyrophilin coding genes associated with IFN-γ(PWM), and IL-17 coding genes associated with IgA. Due to the importance of the MHC in the control of immune responses, these genes may play an important role in resistance to parasitic infections. Our results reveal a largely complex and polygenic genetic profile of the studied traits in this Scottish Blackface sheep population.

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.

Genome-wide association study of health and production traits in meat sheep

Genotypes are currently widely used in animal breeding programmes to enhance the speed of genetic progress. With sufficient data, a Genome-Wide Association Study (GWAS) can be performed to identify informative markers. The aim of this study was to investigate the genetic background of health (footrot and mastitis) and production (birth weight, weaning weight, scan weight, and fat and muscle depth) traits using the available phenotypic and Single Nucleotide Polymorphism (SNP) data collected on the UK Texel sheep population. Initially, 10 193 genotypes were subject to quality control, leaving 9 505 genotypes for further analysis. Selected genotypes, recorded on four different Illumina chip types from low density (15 k SNPs) to high density (606 006 SNPs), were imputed to a subset of 45 686 markers from 50 k array, distributed on 27 chromosomes. Phenotypes collected on 32 farms across the UK for footrot and mastitis and extracted from the UK National database (iTexel) for the production traits were used along with pre-estimated variance components to obtain de-regressed breeding values and used to perform GWAS. Results showed three SNPs being significant on the genome-wise level ('OAR8_62240378.1' on chromosome 8 for birth weight, 's14444.1' on chromosome 19 for weaning weight and 's65197.1' on chromosome 23 for scan weight). Fourteen subsequent SNPs were found to be significant at the chromosome-wise level. These SNPs are located within or close to previously reported QTLs impacting on animal health (such as faecal egg count or somatic cell count) and production (such as body or carcass weight and fat amount). These results indicate that the studied traits are highly polygenic with complex genetic architecture.

Experimental infection with Neospora caninum in Texel ewes at different stages of gestation

In this study Texel sheep, at different stages of pregnancy, were experimentally infected with Neospora caninum. Eleven ewes, seronegative for N. caninum and Toxoplasma gondii, were inoculated 30 days before breeding (Group A), or at 65, 100, and 120 days of gestation (Groups B, C, and D). The group E (control) was inoculated with PBS. Blood samples were collected at -2, 2, 5, and 7 days post-infection (dpi), and weekly up to 42 dpi, for hematology, parasitemia (PCR), and serology (RIFI) assessments. Blood and tissue samples were collected from the lambs for molecular and histological analyses. All animals in Groups B, C, and D were seroconverted, whilst those in groups A and E remained seronegative. Parasitic DNA was detected in the blood of two ewes (groups B and D) and a lamb (group D), and in the brain of a lamb (group B). The parasitemia-positive ewe in group B delivered a weak and seropositive lamb, and had parasitic DNA in its placenta. These results confirm the vertical transmission of N. caninum in ewes inoculated at the beginning and end of pregnancy. The absence of abortions and other clinical signs suggest that Texel sheep may potentially have resistance to N. caninum.

Risk factors associated with the infection of sheep with Dichelobacter nodosus

Ovine footrot is a highly contagious foot disease caused by the gram-negative bacterium Dichelobacter nodosus (D. nodosus). In a recent report, we showed a prevalence of 42.9% D. nodosus positive swabs across Germany. In this follow-up study, we used real-time PCR results for D. nodosus and footrot scores of 9297 sheep from 208 flocks and collated these data with survey data on herd and animal characteristics and herd management. The aims of the present study were to investigate herd and animal factors associated with D. nodosus infection and footrot scores in individual sheep. Multivariable analyses with generalized mixed models showed that month of recording, breed, herdbook membership, use of antibiotics, and footbaths in the past 3–10 years, signs of footrot in the past 12 months and flock environment of the sheep, modelled as a random farm effect within region, were significant risk factors. Among the 21 different breeds, Romney had the lowest risk of D. nodosus infection, while Swifter had the highest risk and German Merino and German White Heath were the next breeds at highest risk of D. nodosus infection. The variance between farms in the prevalence of D. nodosus was large and accounted for 84% of the total variance in the mixed model analysis. We conclude that specific and as yet unknown effects influencing D. nodosus infections in flocks, as well as breed and weather, are the most important effects on D. nodosus infection in sheep, pointing towards the need to establish adequate infection control at farm level.

Screening of stable resistant accessions and identification of resistance loci to Barley yellow mosaic virus disease

Background

The disease caused by Barley yellow mosaic virus (BaYMV) infection is a serious threat to autumn-sown barley (Hordeum vulgare L.) production in Europe, East Asia and Iran. Due to the rapid diversification of BaYMV strains, it is urgent to discover novel germplasm and genes to assist breeding new varieties with resistance to different BaYMV strains, thus minimizing the effect of BaYMV disease on barley cropping.

Methods

A natural population consisting of 181 barley accessions with different levels of resistance to BaYMV disease was selected for field resistance identification in two separate locations (Yangzhou and Yancheng, Jiangsu Province, China). Additive main effects and multiplicative interaction (AMMI) analysis was used to identify accessions with stable resistance. Genome-wide association study (GWAS) of BaYMV disease resistance was broadly performed by combining both single nucleotide polymorphisms (SNPs) and specific molecular markers associated with the reported BaYMV disease resistance genes. Furthermore, the viral protein genome linked (VPg) sequences of the virus were amplified and analyzed to assess the differences between the BaYMV strains sourced from the different experimental sites.

Results

Seven barley accessions with lower standardized Area Under the Disease Progress Steps (sAUDPS) index in every environment were identified and shown to have stable resistance to BaYMV disease in each assessed location. Apart from the reported BaYMV disease resistance genes rym4 and rym5, one novel resistance locus explaining 24.21% of the phenotypic variation was identified at the Yangzhou testing site, while two other novel resistance loci that contributed 19.23% and 19.79% of the phenotypic variation were identified at the Yancheng testing site, respectively. Further analysis regarding the difference in the VPg sequence of the predominant strain of BaYMV collected from these two testing sites may explain the difference of resistance loci differentially identified under geographically distinct regions. Our research provides novel genetic resources and resistance loci for breeding barley varieties for BaMYV disease resistance.

British Sheep Breeds as a Part of World Sheep Gene Pool Landscape: Looking into Genomic Applications

Sheep farming has been an important sector of the UK's economy and rural life for many centuries. It is the favored source of wool, meat and milk products. In the era of exponential progress in genomic technologies, we can now address the questions of what is special about UK sheep breed genotypes and how they differ genetically form one another and from other countries. We can reflect how their natural history has been determined at the level of their genetic code and what traces have been left in their genomes because of selection for phenotypic traits. These include adaptability to certain environmental conditions and management, as well as resistance to disease. Application of these advancements in genetics and genomics to study sheep breeds of British domestic selection has begun and will continue in order to facilitate conservation solutions and production improvement.

Ovine footrot: A review of current knowledge

Footrot is a contagious foot disease mainly affecting sheep. It is caused by the Gram-negative anaerobic bacterium Dichelobacter nodosus. Warm, wet environmental conditions favour development of footrot, and under perfect conditions, it takes just 2-3 weeks from infection to manifestation of clinical signs. Affected sheep show lameness of various degrees and often graze while resting on their carpi. Local clinical signs vary in severity and extent from interdigital inflammation (benign footrot) to underrunning of the complete horn shoe in advanced stages of virulent footrot. Laboratory diagnosis ideally involves collection of four-foot interdigital swab samples followed by competitive real time PCR, allowing for detection of the presence of D. nodosus and differentiation between benign and virulent strains. Laboratory-based diagnostics at the flock level based on risk-based sampling and pooling of interdigital swab samples are recommended. The list of treatment options of individual sheep includes careful removal of the loose undermined horn, local or systemic administration of antimicrobials, systemic administration of non-steroidal anti-inflammatories (NSAIDs) and disinfectant footbathing. Strategies for control at the flock level are manifold and depend on the environmental conditions and the procedures traditionally implemented by the respective country. Generally, measures consist of treatment/culling of infected sheep, vaccination and prevention of reinfection of disease-free flocks. Gaining deeper insight into the beneficial effects of NSAIDs, screening for eco-friendly footbath solutions, developing better vaccines, including the development of a robust, reproducible infection model and elucidation of protective immune responses, as well as the elaboration of effective awareness training programs for sheep farmers, are relevant research gaps.

Whole genome studies of origin, selection and adaptation of the Russian cattle breeds

Our review presents several recent studies on the genetic history and signatures of selection in genomes of the native Russian cattle breeds. Most of these works are not easily accessible for the Russian-speaking audience. We describe the origins of appearance of the Russian cattle breeds from the genetics perspective. We point to the links between most of the Russian breeds with the taurine breeds of the European origin and for some Russian breeds with the breeds of the Asian origin. We describe major phylogenetic clusters of the Russian breeds and point to those that still maintain their unique genetics, meaning that their preservation is a priority. In addition, we review the results of the search for signatures of selection in genomes of the Russian cattle breeds. Some unique signatures of selection present in the genomes of so-called “turano-mongolian” cattle (i. e. the Yakut cattle) are described which allowed the Yakut cattle to adapt to harsh environments found above the Polar Circle. Signatures of selection which could help other cattle breeds of the Russian origin to adapt to various climatic condition of the Russian Federation are reviewed. The Russian cattle genomes also contain known signatures of selection related to cattle domestication about 8–10 thousand years ago. The most profound ones include genes related to changes of the coat colour. This phenotype in many cases could be related to the distinction of the first domesticated populations and lead to the formation of so-called land races (primitive breeds). Whole-genome association studies of Russian cattle breeds pointed to a novel gene which could be related to the “white-faced” phenotype and to a gene which is related to body temperature support under the acute cold stress. The data presented in our review could be used for identification of genetic markers to focus on in future efforts on designing new highly productive cattle breeds adapted to climates of the Russian Federation and other countries with similar climates.

Aetiology, Risk Factors, Diagnosis and Control of Foot-Related Lameness in Dairy Sheep

During the last twenty years, considerable research efforts have recognized the consequences of foot-related lameness primarily in cattle, and meat and wool sheep. Despite the lack of extensive epidemiological studies, field observations and isolated research reports in dairy sheep have suggested that the problem might be more severe in semi-intensive and intensive farming systems. Footrot, contagious ovine digital dermatitis, ovine interdigital dermatitis, white line disease, and pedal joint abscess are the most common causes of foot-related lameness. Dichelobacter nodosus, Fusobacterium necrophorum, Treponema spp., and Actinomyces pyogenes are the most significant foot-related lameness-associated pathogens. Despite a documented hereditary predisposition, environmental factors are the most important in determining the occurrence of foot-related lameness. Moist and warm environment, increased parity and milk yield, inappropriate housing conditions and infrastructures, inadequate hygiene status, imbalanced nutrition, and insufficient foot care are the most critical risk factors. Furthermore, a foot-lameness control plan should include targeted implementation of claw trimming and footbathing, evidence-based planning of hygiene measures in preventive veterinary practices (i.e., antibiotic administration, vaccinations against footrot), selective breeding to footrot resistance, and, most importantly, the continuous training of farming personnel. Controlling foot-lameness in dairy sheep is critical in determining the well-being of animals, and strongly affects the farm's profitability and sustainability.

Enhancing genetic disease control by selecting for lower host infectivity and susceptibility

Infectious diseases have a huge impact on animal health, production and welfare, and human health. Understanding the role of host genetics in disease spread is important for developing disease control strategies that efficiently reduce infection incidence and risk of epidemics. While heritable variation in disease susceptibility has been targeted in livestock breeding, emerging evidence suggests that there is additional genetic variation in host infectivity, but the potential benefits of including infectivity into selection schemes are currently unknown. A Susceptible-Infected-Recovered epidemiological model incorporating polygenic genetic variation in both susceptibility and infectivity was combined with quantitative genetics selection theory to assess the non-linear impact of genetic selection on field measures of epidemic risk and severity. Response to 20 generations of selection was calculated in large simulated populations, exploring schemes differing in accuracy and intensity. Assuming moderate genetic variation in both traits, 50% selection on susceptibility required seven generations to reduce the basic reproductive number R₀ from 7.64 to the critical threshold of <1, below which epidemics die out. Adding infectivity in the selection objective accelerated the decline towards R₀ < 1, to 3 generations. Our results show that although genetic selection on susceptibility reduces disease risk and prevalence, the additional gain from selection on infectivity accelerates disease eradication and reduces more efficiently the risk of new outbreaks, while it alleviates delays generated by unfavourable correlations. In conclusion, host infectivity was found to be an important trait to target in future genetic studies and breeding schemes, to help reducing the occurrence and impact of epidemics.

Genome-wide genotyping uncovers genetic profiles and history of the Russian cattle breeds

One of the most economically important areas within the Russian agricultural sector is dairy and beef cattle farming contributing about $11 billion to the Russian economy annually. Trade connections, selection and breeding have resulted in the establishment of a number of breeds that are presumably adapted to local climatic conditions. Little however is known about the ancestry and history of Russian native cattle. To address this question, we genotyped 274 individuals from 18 breeds bred in Russia and compared them to 135 additional breeds from around the world that had been genotyped previously. Our results suggest a shared ancestry between most of the Russian cattle and European taurine breeds, apart from a few breeds that shared ancestry with the Asian taurines. The Yakut cattle, belonging to the latter group, was found to be the most diverged breed in the whole combined dataset according to structure results. Haplotype sharing further suggests that the Russian cattle can be divided into four major clusters reflecting ancestral relations with other breeds. Herein, we therefore shed light on to the history of Russian cattle and identified closely related breeds to those from Russia. Our results will facilitate future research on detecting signatures of selection in cattle genomes and eventually inform future genetics-assisted livestock breeding programs in Russia and in other countries.

A genome-wide significant association on chromosome 2 for footrot resistance/susceptibility in Swiss White Alpine sheep

Footrot is one of the most important causes of lameness in global sheep populations and is characterized by a bacterial infection of the interdigital skin. As a multifactorial disease, its clinical representation depends not only on pathogen factors and environmental components but also on the individual resistance/susceptibility of the host. A genetic component has been shown in previous studies; however, so far no causative genetic variant influencing the risk of developing footrot has been identified. In this study, we genotyped 373 Swiss White Alpine sheep, using the ovine high-density 600k SNP chip, in order to run a DNA-based comparison of individuals with known clinical footrot status. We performed a case-control genome-wide association study, which revealed a genome-wide significant association for SNP rs418747104 on ovine chromosome 2 at 81.2 Mb. The three best associated SNP markers were located at the MPDZ gene, which codes for the multiple PDZ domain crumbs cell polarity complex component protein, also known as multi-PDZ domain protein 1 (MUPP1). This protein is possibly involved in maintaining the barrier function and integrity of tight junctions. Therefore, we speculate that individuals carrying MPDZ variants may differ in their footrot resistance/susceptibility due to modified horn and interdigital skin integrity. In conclusion, our study reveals that MPDZ might represent a functional candidate gene, and further research is needed to explore its role in footrot affected sheep.

Detection of QTL for greasy fleece weight in sheep using a 50 K single nucleotide polymorphism chip

Genome-wide association studies (GWAS) have introduced an influential tool in the search for quantitative trait loci (QTL) influencing economically important traits in sheep. To identify QTL associated with greasy fleece weight, a GWAS with 50 K single nucleotide polymorphisms (SNPs) was performed in a Baluchi sheep population. Association with greasy fleece weights was tested using the software Plink. The results of our GWAS provided three novel SNP markers and candidate genes associated with greasy fleece weight. A total of three chromosome-wide significant associations were detected for SNP on chromosomes 17 and 20 affecting greasy fleece weight across the four shearing. One of the significant SNP markers was located within ovine known genes namely FAM101A. Further investigation of these identified regions in validation studies will facilitate the identification of strong candidate genes for wool production in sheep.

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

BACKGROUND

Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.

RESULTS

Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.

CONCLUSIONS

Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.

The Use of SNP Microarrays for Biodiversity Studies of Sheep – A Review

The conservation of farm animal genetic resources and their protection against genetic erosion requires knowledge of biodiversity status. Genetic variation in populations can be estimated using both traditional pedigree-based methods and molecular techniques. SNP microarrays are a new generation of molecular genetic tools, which have found application in analysis of biodiversity in populations of domestic and wild sheep, in studies of resistance to intestinal parasites and foot rot, and in searching for markers associated with meat and milk yield, or colour inheritance traits. The aim of the study is the review of recent literature on the biodiversity and the use of molecular markers for population genetics in different breeds and populations of sheep.

Genetic diversity analyses reveal first insights into breed-specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH ) with pedigree-based inbreeding coefficients (F_PED ). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.