A metagenomics approach to characterize the footrot microbiome in Merino sheep

Metagenomic profiling of viral and microbial communities from the pox lesions of lumpy skin disease virus and sheeppox virus-infected hosts

Introduction

It has been recognized that capripoxvirus infections have a strong cutaneous tropism with the manifestation of skin lesions in the form of nodules and scabs in the respective hosts, followed by necrosis and sloughing off. Considering that the skin microbiota is a complex community of commensal bacteria, fungi and viruses that are influenced by infections leading to pathological states, there is no evidence on how the skin microbiome is affected during capripoxvirus pathogenesis.

Methods

In this study, shotgun metagenomic sequencing was used to investigate the microbiome in pox lesions from hosts infected with lumpy skin disease virus and sheep pox virus.

Results

The analysis revealed a high degree of variability in bacterial community structures across affected skin samples, indicating the importance of specific commensal microorganisms colonizing individual hosts. The most common and abundant bacteria found in scab samples were Fusobacterium necrophorum, Streptococcus dysgalactiae, Helcococcus ovis and Trueperella pyogenes, irrespective of host. Bacterial reads belonging to the genera Moraxella, Mannheimia, Corynebacterium, Staphylococcus and Micrococcus were identified.

Discussion

This study is the first to investigate capripox virus-associated changes in the skin microbiome using whole-genome metagenomic profiling. The findings will provide a basis for further investigation into capripoxvirus pathogenesis. In addition, this study highlights the challenge of selecting an optimal bioinformatics approach for the analysis of metagenomic data in clinical and veterinary practice. For example, direct classification of reads using a kmer-based algorithm resulted in a significant number of systematic false positives, which may be attributed to the peculiarities of the algorithm and database selection. On the contrary, the process of de novo assembly requires a large number of target reads from the symbiotic microbial community. In this work, the obtained sequencing data were processed by three different approaches, including direct classification of reads based on k-mers, mapping of reads to a marker gene database, and de novo assembly and binning of metagenomic contigs. The advantages and disadvantages of these techniques and their practicality in veterinary settings are discussed in relation to the results obtained.

Genome-wide association study identifies genetic variants underlying footrot in Portuguese Merino sheep

BACKGROUND

Ovine footrot caused by Dichelobacter nodosus (D. nodosus) is a contagious disease with serious economic and welfare impacts in sheep production systems worldwide. A better understanding of the host genetic architecture regarding footrot resistance/susceptibility is crucial to develop disease control strategies that efficiently reduce infection and its severity. A genome-wide association study was performed using a customized SNP array (47,779 SNPs in total) to identify genetic variants associated to footrot resistance/susceptibility in two Portuguese native breeds, i.e. Merino Branco and Merino Preto, and a population of crossbred animals. A cohort of 1375 sheep sampled across 17 flocks, located in the Alentejo region (southern Portugal), was included in the analyses.

RESULTS

Phenotypes were scored from 0 (healthy) to 5 (severe footrot) based on visual inspection of feet lesions, following the Modified Egerton System. Using a linear mixed model approach, three SNPs located on chromosome 24 reached genome-wide significance after a Bonferroni correction (p < 0.05). Additionally, six genome-wide suggestive SNPs were identified each on chromosomes 2, 4, 7, 8, 9 and 15. The annotation and KEGG pathway analyses showed that these SNPs are located within regions of candidate genes such as the nonsense mediated mRNA decay associated PI3K related kinase (SMG1) (chromosome 24) and the RALY RNA binding protein like (RALYL) (chromosome 9), both involved in immunity, and the heparan sulfate proteoglycan 2 (HSPG2) (chromosome 2) and the Thrombospodin 1 (THBS1) (chromosome 7) implicated in tissue repair and wound healing processes.

CONCLUSION

This is the first attempt to identify molecular markers associated with footrot in Portuguese Merino sheep. These findings provide relevant information on a likely genetic association underlying footrot resistance/susceptibility and the potential candidate genes affecting this trait. Genetic selection strategies assisted on the information obtained from this study could enhance Merino sheep-breeding programs, in combination with farm management strategies, for a more effective and sustainable long-term solution for footrot control.

Shotgun Metagenomics Reveals Minor Micro"bee"omes Diversity Defining Differences between Larvae and Pupae Brood Combs

Bees represent not only a valuable asset in agriculture, but also serve as a model organism within contemporary microbiology. The metagenomic composition of the bee superorganism has been substantially characterized. Nevertheless, traditional cultural methods served as the approach to studying brood combs in the past. Indeed, the comb microbiome may contribute to determining larval caste differentiation and hive immunity. To further this understanding, we conducted a shotgun sequencing analysis of the brood comb microbiome. While we found certain similarities regarding species diversity, it exhibits significant differentiation from all previously described hive metagenomes. Many microbiome members maintain a relatively constant ratio, yet taxa with the highest abundance level tend to be ephemeral. More than 90% of classified metagenomes were Gammaproteobacteria, Bacilli and Actinobacteria genetic signatures. Jaccard dissimilarity between samples based on bacteria genus classifications hesitate from 0.63 to 0.77, which for shotgun sequencing indicates a high consistency in bacterial composition. Concurrently, we identified antagonistic relationships between certain bacterial clusters. The presence of genes related to antibiotic synthesis and antibiotic resistance suggests potential mechanisms underlying the stability of comb microbiomes. Differences between pupal and larval combs emerge in the total metagenome, while taxa with the highest abundance remained consistent. All this suggests that a key role in the functioning of the comb microbiome is played by minor biodiversity, the function of which remains to be established experimentally.