An in-depth investigation of the microbiota and its virulence factors associated with severe udder cleft dermatitis lesions.

Udder cleft dermatitis (UCD) is a skin condition affecting the anterior parts of the udder in dairy cattle. In the present study, we aimed to shed light on the microbiota in severe UCD lesions versus healthy udder skin by putting forward a taxonomic and functional profile based on a virulence factor analysis. Through shotgun metagenomic sequencing, we found a high proportion of bacteria in addition to a low abundance of archaea. A distinct clustering of healthy udder skin versus UCD lesion samples was shown by applying principal component analysis and (sparse) partial least squares analysis on the metagenomic data. Proteobacteria, Bacillota, and Actinomycetota were among the most abundant phyla in healthy udder skin samples. In UCD samples, Bacteroidota was the most abundant phylum. At genus level, Bifidobacterium spp. had the highest relative abundance in healthy skin samples, whereas Porphyromonas spp. and Corynebacterium spp. had the highest relative abundance in UCD samples. In the differential abundance analysis, Porphyromonas spp. and Bacteroides spp. were significantly differentially abundant in UCD samples, whereas Bifidobacterium spp., Staphylococcus sp. AntiMn-1, and Staphylococcus equorum were more commonly found in healthy samples. Moreover, the abundance of several treponeme phylotypes was significantly higher in lesion samples. The streptococcal cysteine protease speB was among the most abundant virulence factors present in severe UCD lesions, while a plethora of virulence factors such as the antitoxin relB were downregulated, possibly contributing to creating the ideal wound climate for the dysbiotic community. Network analysis showed healthy lesion samples had a large network ofpositive, correlations between the abundances of beneficial species such as Aerococcus urinaeequi and Bifidobacterium angulatum, indicating that the healthy skin microbiome forms an active protective bacterial network, which is disrupted in case of UCD. In UCD samples, a smaller microbial network mainly consisting of positive correlations between the abundances of Bacteroides fragilis and anaerobic Bacteroidota was exposed. Moreover, a high correlation between the taxonomic data and virulence factors was revealed, concurrently with 2 separate networks of microbes and virulence factors. One network, matching with the taxonomic findings in the healthy udder skin samples, showcased a community of harmless or beneficial bacteria, such as Bifidobacterium spp. and Butyrivibrio proteoclasticus, associated with hcnB, hcnC, relB, glyoxalase, and cupin 2. The other network, corresponding to UCD samples, consisted of pathogenic or facultative pathogenic and mainly anaerobic bacteria such as Treponema spp., Mycoplasmopsis spp., and bovine gammaherpesvirus 4, that correlated with virulence factors SpvB, fhaB, and haemagglutination activity domain-associated factor. Our results point toward a dysbiotic community with a notable decrease in diversity and evenness, with a loss of normal skin inhabitants and innocuous or useful species making way for predominantly anaerobic, facultative pathogens. The shift in the abundance of virulence factors such as fhaB and SpvB could play a role in the manifestation of a local micro-environment favorable to the microbiome associated with udder skin lesions. Lastly, the presence of specific networks between microbial species, and between microbes and virulence factors was shown.

Severe udder cleft dermatitis lesion transcriptomics points to an impaired skin barrier, defective wound repair and a dysregulated inflammatory response as key elements in the pathogenesis.

This study is the first to investigate the transcriptomic changes occurring in severe udder cleft dermatitis lesions (UCD) in Holstein-Friesian cows. An examination of the gene expression levels in natural UCD lesions and healthy udder skin through RNA Seq-Technology provided a deeper insight into the inflammatory pathways associated with this disease. A clear distinction between the gene expression patterns of UCD lesions and healthy skin was shown in the principal component analysis. Genes coding for inflammatory molecules were upregulated such as the chemokines C-X-C motif ligand 2 (CXCL2), 5 (CXCL5) and 8 (CXCL8), and C-C motif ligand 11 (CCL11). Moreover, the genes coding for the multifunctional molecules ADAM12 and SLPI were amongst the highest upregulated ones, whereas the most downregulated genes included the ones coding for keratins and keratin-associated molecules. Predominantly inflammatory pathways such as the chemokine signaling, cytokine receptor interaction and IL-17 signaling pathway were significantly upregulated in the pathway analysis. These results point towards a fulminant, dysregulated inflammatory response concomitant with a disruption of the skin barrier integrity and a hampered wound repair mechanism in severe UCD lesions.