The forgotten variable? Does the euthanasia method and sample storage condition influence an organisms transcriptome - a gene expression analysis on multiple tissues in pigs.

BACKGROUND: Transcriptomic studies often require collection of fresh tissues post euthanasia. The chosen euthanasia method might have the potential to induce variations in gene expressions that are unlinked with the experimental design. The present study compared the suitability of 'nitrogen gas in foam' (ANOXIA) in comparison to a non-barbiturate anaesthetic, T-61® (T61), for euthanizing piglets used in transcriptome research. Further, the effect of common tissue storage conditions, RNAlater™ (RL) and snap freezing in liquid nitrogen (LN2), on gene expression profiles were also analysed. RESULTS: On comparison of the 3'mRNA-Seq data generated from pituitary, hypothalamus, liver and lung tissues, no significant differential expression in the protein coding genes were detected between the euthanasia methods. This implies that the nitrogen anoxia method could be a suitable alternative for euthanasia of piglets used in transcriptomic research. However, small nuclear RNAs (snRNAs) that constitute the eukaryotic spliceosomal machinery were found to be significantly higher (log2fold change ≥ 2.0, and adjusted p value ≤ 0.1) in pituitary samples collected using ANOXIA. Non-protein coding genes like snRNAs that play an important role in pre-mRNA splicing can subsequently modify gene expression. Storage in RL was found to be superior in preserving RNA compared to LN2 storage, as evidenced by the significantly higher RIN values in representative samples. However, storage in RL as opposed to LN2, also influenced differential gene expression in multiple tissues, perhaps as a result of its inability to inhibit biological activity during storage. Hence such external sources of variations should be carefully considered before arriving at research conclusions. CONCLUSIONS: Source of biological variations like euthanasia method and storage condition can confound research findings. Even if we are unable to prevent the effect of these external factors, it will be useful to identify the impact of these variables on the parameter under observation and thereby prevent misinterpretation of our results.

Risk of contamination of semen, vaginal secretions, follicular fluid and ovarian medulla with SARS-CoV-2 in patients undergoing ART.

STUDY QUESTION: Can severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mRNA be detected in the reproductive tract of asymptomatic patients undergoing ART? SUMMARY ANSWER: SARS-CoV-2 mRNA is not detectable in semen, follicular fluid, vaginal secretions or residual medulla from ovarian tissue cryopreservation procedures in asymptomatic patients who undergo ART, irrespective of the results of a triage questionnaire and a nasopharyngeal SARS-CoV-2 RNA detection test. WHAT IS KNOWN ALREADY: The SARS-CoV-2 pandemic had a huge impact on the activities of fertility clinics. Although some studies reported the presence of SARS-CoV-2 mRNA in the reproductive system during or after acute COVID-19 symptomatic infections, uncertainties remain regarding the presence of viral mRNA in the reproductive material and follicular fluid of asymptomatic patients undergoing ART. STUDY DESIGN, SIZE, DURATION: An observational cohort trial of residual material samples including semen, follicular fluid, vaginal secretions and ovarian medulla was conducted during the second pandemic wave in Brussels from September 2020 to April 2021. PARTICIPANTS/MATERIALS, SETTING, METHODS: All patients who underwent ART (IUI, IVF/ICSI, oocyte and ovarian tissue cryopreservation) responded to a triage questionnaire at the beginning and end of the cycle and underwent nasopharyngeal swab collection for SARS-CoV-2 RNA detection by RT-PCR before the procedure according to standard recommendations. For semen analysis, only the questionnaire was requested the day before the sample collection. The ART cycles of patients with positive nasopharyngeal SARS-CoV-2 RNA detection tests and/or questionnaires were cancelled except for those that could not be postponed. After providing informed consent, swabs on residual materials were collected the day of the oocyte, ovarian tissue or semen collection and were processed for RT-qPCR. MAIN RESULTS AND THE ROLE OF CHANCE: A total of 394 samples from 291 patients were analysed. Amongst them, 20 samples were obtained from patients with a positive questionnaire but negative nasopharyngeal SARS-CoV-2 test and 20 others were from patients with a positive nasopharyngeal SARS-CoV-2 test. The remaining samples were collected from patients with a negative or unknown nasopharyngeal SARS-CoV-2 test and/or a negative or unknown triage questionnaire. Viral RNA for SARS-CoV-2 was undetectable in all of the samples. LIMITATIONS, REASONS FOR CAUTION: Considering the cancellation policy, only a limited number of samples from patients with positive triage questionnaires or nasopharyngeal SARS-CoV-2 tests were included in the analysis. WIDER IMPLICATIONS OF THE FINDINGS: The study suggested that there was no risk of reproductive tract contamination by SARS-CoV-2 in asymptomatic patients, irrespective of the results from a triage questionnaire or nasopharyngeal SARS-CoV-2 test. The results suggested that no additional measures to prevent staff or cross-patient contamination need to be implemented in the IVF and andrology laboratories. STUDY FUNDING/COMPETING INTEREST(S): This study was funded by the Université Libre de Bruxelles and by a grant from Ferring. A.D. and I.D. received a grant from Ferring for the study. The authors have no other conflict of interest to declare related to this study. TRIAL REGISTRATION NUMBER: N/A.

How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species.

BACKGROUND: PLINK is probably the most used program for analyzing SNP genotypes and runs of homozygosity (ROH), both in human and in animal populations. The last decade, ROH analyses have become the state-of-the-art method for inbreeding assessment. In PLINK, the --homozyg function is used to perform ROH analyses and relies on several input settings. These settings can have a large impact on the outcome and default values are not always appropriate for medium density SNP array data. Guidelines for a robust and uniform ROH analysis in PLINK using medium density data are lacking, albeit these guidelines are vital for comparing different ROH studies. In this study, 8 populations of different livestock and pet species are used to demonstrate the importance of PLINK input settings. Moreover, the effects of pruning SNPs for low minor allele frequencies and linkage disequilibrium on ROH detection are shown. RESULTS: We introduce the genome coverage parameter to appropriately estimate FROH and to check the validity of ROH analyses. The effect of pruning for linkage disequilibrium and low minor allele frequencies on ROH analyses is highly population dependent and such pruning may result in missed ROH. PLINK's minimal density requirement is crucial for medium density genotypes and if set too low, genome coverage of the ROH analysis is limited. Finally, we provide recommendations for the maximal gap, scanning window length and threshold settings. CONCLUSIONS: In this study, we present guidelines for an adequate and robust ROH analysis in PLINK on medium density SNP data. Furthermore, we advise to report parameter settings in publications, and to validate them prior to analysis. Moreover, we encourage authors to report genome coverage to reflect the ROH analysis' validity. Implementing these guidelines will substantially improve the overall quality and uniformity of ROH analyses.

SNP genotypes reveal breed substructure, selection signatures and highly inbred regions in Piétrain pigs.

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed's genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and FST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population-specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34-126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.

Unraveling the genetic diversity of Belgian Milk Sheep using medium-density SNP genotypes.

The present study focuses on the Belgian Milk Sheep in Flanders (Belgium) and compares its genetic diversity and relationship with the Flemish Sheep, the Friesian Milk Sheep, the French Lacaune dairy sheep and other Northern European breeds. For this study, 94 Belgian Milk Sheep, 23 Flemish Sheep and 22 Friesian Milk Sheep were genotyped with the OvineSNP50 array. In addition, 29 unregistered animals phenotypically similar to Belgian Milk Sheep were genotyped using the 15K ISGC chip. Both Belgian and Friesian Milk Sheep as well as the East Friesian Sheep were found to be less diverse than the other seven breeds included in this study. Genomic inbreeding coefficients based on runs of homozygosity (ROH) were estimated at 14.5, 12.4 and 10.2% for Belgian Milk Sheep, Flemish Sheep and Friesian Milk Sheep respectively. Out of 29 unregistered Belgian Milk Sheep, 28 mapped in the registered Belgian Milk Sheep population. Ancestry analysis, PCA and FST calculations showed that Belgian Milk Sheep are more related to Friesian Milk Sheep than to Flemish Sheep, which was contrary to the breeders' expectations. Consequently, breeders may prefer to crossbreed Belgian Milk Sheep with Friesian sheep populations (Friesian Milk Sheep or East Friesian Sheep) in order to increase diversity. This research underlines the usefulness of SNP chip genotyping and ROH analyses for monitoring genetic diversity and studying genetic links in small livestock populations, profiting from internationally available genotypes. As assessment of genetic diversity is vital for long-term breed survival, these results will aid flockbooks to preserve genetic diversity.

Nectar traits differ between pollination syndromes in Balsaminaceae.

BACKGROUND AND AIMS: The attractiveness of nectar rewards depends both on the quantity of nectar produced and on its chemical composition. It is known that nectar quantity and chemical composition can differ in plant species depending on the main pollinator associated with the species. The main aims of this study were to test formally whether nectar traits are adapted to pollination syndromes in the speciose Balsaminaceae and, if so, whether a combination of nectar traits mirrors pollination syndromes. METHODS: Comparative methods based on Ornstein-Uhlenbeck models were used to test whether nectar volume, nectar sucrose proportion, sugar and amino acid concentration and amino acid composition had evolved as a function of pollination syndromes in 57 species of Balsaminaceae. Cluster analysis and ordination were performed to derive clusters of species resembling each other in nectar composition. KEY RESULTS: Evolutionary models for nectar volume and nectar sucrose proportion performed best when including information on pollination syndrome, while including such information improve model fit neither for sugar and amino acid concentration nor for amino acid composition. A significant relationship emerged between pollination syndrome and the combined nectar traits. CONCLUSIONS: Our results show that nectar volume and nectar sucrose proportion evolve rapidly towards optimal values associated with different pollination syndromes. The detection of a signal indicating that nectar traits in combination are to a certain extent able to predict pollination syndromes in Balsaminaceae suggests that a holistic approach including the whole set of nectar traits helps us to better understand evolution of nectar composition in response to pollinators.

Expanding the clinical spectrum of biallelic ZNF335 variants.

ZNF335 plays an essential role in neurogenesis and biallelic variants in ZNF335 have been identified as the cause of severe primary autosomal recessive microcephaly in 2 unrelated families. We describe, herein, 2 additional affected individuals with biallelic ZNF335 variants, 1 individual with a homozygous c.1399 T > C, p.(Cys467Arg) variant, and a second individual with compound heterozygous c.2171_2173delTCT, p.(Phe724del) and c.3998A > G, p.(Glu1333Gly) variants with the latter variant predicted to affect splicing. Whereas the first case presented with early death and a severe phenotype characterized by anterior agyria with prominent extra-axial spaces, absent basal ganglia, and hypoplasia of the brainstem and cerebellum, the second case had a milder clinical presentation with hypomyelination and otherwise preserved brain structures on MRI. Our findings expand the clinical spectrum of ZNF335-associated microcephaly.

Refractive index gratings in electro-optic polymer thin films.

Refractive index gratings have been inscribed in polymer thin films by permanently photobleaching the organic chromophore PYR-3 dopant. The grating inscription process was investigated in detail for the purpose of improving the diffraction efficiency (η) of the PYR-3 doped polymer gratings. Three processes were identified that contributed to the η of the first diffracted order: a periodic change in the refractive index due to photobleaching of the PYR-3, formation of the surface relief grating as a consequence of free volume change during bleaching, and the introduction of periodic, strain-induced changes in the refractive index.

Half of 23 Belgian dog breeds has a compromised genetic diversity, as revealed by genealogical and molecular data analysis.

The genetic diversity in 23 dog breeds raised in Belgium was investigated using both genealogical analysis and microsatellite markers. Some of these breeds are native breeds, with only small populations maintained. Pedigree and molecular data, obtained from the Belgian kennel club, were used to calculate the inbreeding coefficients, realised effective population size as well as probabilities of gene origin and average observed heterozygosity. Inbreeding coefficients ranged from 0.8 to 44.7% and realised effective population size varied between 3.2 and 829.1, according to the used method and breed. Mean observed heterozygosity ranged from 0.47 to 0.73. Both pedigree and molecular methods reveal low genetic diversity and presence of bottlenecks, especially in native Belgian breeds with small population sizes. Furthermore, principal component analysis on the set of investigated diversity parameters revealed no groups of breeds that could be identified in which similar breeding strategies could be applied to maintain genetic diversity.

Developing 3D SEM in a broad biological context.

When electron microscopy (EM) was introduced in the 1930s it gave scientists their first look into the nanoworld of cells. Over the last 80 years EM has vastly increased our understanding of the complex cellular structures that underlie the diverse functions that cells need to maintain life. One drawback that has been difficult to overcome was the inherent lack of volume information, mainly due to the limit on the thickness of sections that could be viewed in a transmission electron microscope (TEM). For many years scientists struggled to achieve three-dimensional (3D) EM using serial section reconstructions, TEM tomography, and scanning EM (SEM) techniques such as freeze-fracture. Although each technique yielded some special information, they required a significant amount of time and specialist expertise to obtain even a very small 3D EM dataset. Almost 20 years ago scientists began to exploit SEMs to image blocks of embedded tissues and perform serial sectioning of these tissues inside the SEM chamber. Using first focused ion beams (FIB) and subsequently robotic ultramicrotomes (serial block-face, SBF-SEM) microscopists were able to collect large volumes of 3D EM information at resolutions that could address many important biological questions, and do so in an efficient manner. We present here some examples of 3D EM taken from the many diverse specimens that have been imaged in our core facility. We propose that the next major step forward will be to efficiently correlate functional information obtained using light microscopy (LM) with 3D EM datasets to more completely investigate the important links between cell structures and their functions.

Pick-up and drop transfer of diamond nanosheets.

Nanocrystalline diamond (NCD) is a promising material for electronic and mechanical micro- and nanodevices. Here we introduce a versatile pick-up and drop technique that makes it possible to investigate the electrical, optical and mechanical properties of as-grown NCD films. Using this technique, NCD nanosheets, as thin as 55 nm, can be picked-up from a growth substrate and positioned on another substrate. As a proof of concept, electronic devices and mechanical resonators are fabricated and their properties are characterized. In addition, the versatility of the method is further explored by transferring NCD nanosheets onto an optical fiber, which allows measuring its optical absorption. Finally, we show that NCD nanosheets can also be transferred onto two-dimensional crystals, such as MoS2, to fabricate heterostructures. Pick-up and drop transfer enables the fabrication of a variety of NCD-based devices without requiring lithography or wet processing.

Ecological and evolutionary drivers of range size in Coenagrion damselflies.

Geographic range size is a key ecological and evolutionary characteristic of a species, yet the causal basis of variation in range size among species remains largely unresolved. One major reason for this is that several ecological and evolutionary traits may jointly shape species' differences in range size. We here present an integrated study of the contribution of ecological (dispersal capacity, body size and latitudinal position) and macroevolutionary (species' age) traits in shaping variation in species' range size in Coenagrion damselflies. We reconstructed the phylogenetic tree of this genus to account for evolutionary history when assessing the contribution of the ecological traits and to evaluate the role of the macroevolutionary trait (species' age). The genus invaded the Nearctic twice independently from the Palearctic, yet this was not associated with the evolution of larger range sizes or dispersal capacity. Body size and species' age did not explain variation in range size. There is higher flight ability (as measured by wing aspect ratio) at higher latitudes. Species with a larger wing aspect ratio had a larger range size, also after correcting for phylogeny, suggesting a role for dispersal capacity in shaping the species' ranges. More northern species had a larger species' range, consistent with Rapoport's rule, possibly related to niche width. Our results underscore the importance of integrating macroecology and macroevolution when explaining range size variation among species.

A vestibular schwannoma in a patient with Birt-Hogg-Dube syndrome.

Birt-Hogg Dubé syndrome is an autosomal dominant disease with variable clinical expression. It is characterized by cutaneous manifestations, renal tumors and lung cysts. Other tumors, such as adrenal tumors and tumors originating from the neural crest cells such as meningioma and neurothekeoma have also been described. This syndrome is caused by germline mutations in the folliculin (FLCN) gene located on chromosome 17p. We report, for the first time, a patient with BHDS and a history of a vestibular schwannoma in adolescence. The diagnosis of BHDS was confirmed, by identifying a nonsense mutation in exon 10 of the FLCN gene. A vestibular schwannoma also originates from neural crest cells, just as other neural tumors, previously encountered in patients with BHDS. The reported mutations cause a truncation of the protein, folliculin. The exact role of folliculin is still undetermined. Two different theories suggest the effect of tumorigenesis. One is that folliculin plays an important role in the AMPK-mTOR pathway which leads to proliferation of cells when activated. The other is that the folliculin acts as a possible tumor suppressor gene, since there is a high frequency of second hits in the FLCN-gene. In order to confirm a possible relation of BHDS and neural crest tumors, further research is necessary in the tumorigenesis of the folliculin gene.

Genetic parameters for chronic progressive lymphedema in Belgian Draught Horses.

Genetic parameters for chronic progressive lymphedema (CPL)-associated traits in Belgian Draught Horses were estimated, using a multitrait animal model. Clinical scores of CPL in the four limbs/horse (CPLclin ), skinfold thickness and hair samples (hair diameter) were studied. Due to CPLclin uncertainty in younger horses (progressive CPL character), a restricted data set (D_3+) was formed, excluding records from horses under 3 years from the complete data set (D_full). Age, gender, coat colour and limb hair pigmentation were included as fixed, permanent environment and date of recording as random effects. Higher CPLclin certainty (D_3+) increased heritability coefficients of, and genetic correlations between traits, with CPLclin heritabilities (SE) for the respective data sets: 0.11 (0.06) and 0.26 (0.05). A large proportion of the CPLclin variance was attributed to the permanent environmental effect in D_full, but less in D_3+. Date of recording explained a proportion of variance from 0.09 ± 0.03 to 0.61 ± 0.08. Additive genetic correlations between CPLclin and both skinfold thickness and hair diameter showed the latter two traits cannot be used as a direct diagnostic aid for CPL. Due to the relatively low heritability of CPLclin , selection should focus on estimated breeding values (from repeated clinical examinations) to reduce CPL occurrence in the Belgian Draught Horse.

Effect of sire type and a by-product based diet on performance and meat quality in growing-finishing pigs.

For many years, pig production has focused on maximizing performance by selecting for maximal muscle growth and feeding diets that allow the animals to express their genetic potential. However, it is unclear whether this selection for muscle deposition has affected the capacity of pigs to cope with by-product-based diets, which rely on fat as the primary energy source instead of starches and sugars. Therefore, an experiment was set up to investigate if different types of boars affect how their progeny cope with alternative ingredients in the diet, with a possible need for adapted breeding schemes. Two types of boars within the Piétrain sire line were used based on either a high or low estimated breeding value for daily feed intake (HFI: high feed intake, low feed intake). When their progeny reached 14 weeks of age, two dietary strategies were compared: a control (CON) vs a by-product-based diet high in fat and fiber (HFF). The CON diet was mainly based on cereals (corn, wheat, barley) and soybean meal. The HFF diet was formulated to contain the same net energy, CP and digestible amino acid levels without any cereals or soybean meal. In total 192 animals were included in the experiment (48 animals/type of boar/diet) and performance, digestibility, carcass and meat quality were compared. None of the parameters showed a significant interaction (P < 0.05) between the type of boar and diet, suggesting that shifting to diets that are less prone to feed-food competition is equally feasible in different types of pigs. Type of boar did affect performance, carcass quality and intramuscular fat content. HFI pigs showed higher daily feed intake (DFI) and daily gain (P < 0.001), with no significant difference in feed conversion ratio (P = 0.205), lower carcass quality (P < 0.001) and higher intramuscular fat content (P = 0.030). For both boar types, pigs fed the CON diet performed better, with a higher daily gain (P = 0.028), DFI (P = 0.011) and dressing yield (P = 0.009) and better digestibility (P < 0.001), but without differences in feed conversion ratio or meat quality. In conclusion, there was no indication that pigs differing in feed intake capacity cope differently with a high-fat, high-fiber diet based on by-products. Different types of pigs may cope well with diets that are less prone to feed-food competition.

Hydrogen termination of CVD diamond films by high-temperature annealing at atmospheric pressure.

A high-temperature procedure to hydrogenate diamond films using molecular hydrogen at atmospheric pressure was explored. Undoped and doped chemical vapour deposited (CVD) polycrystalline diamond films were treated according to our annealing method using a H2 gas flow down to ~50 ml∕min (STP) at ~850 °C. The films were extensively evaluated by surface wettability, electron affinity, elemental composition, photoconductivity, and redox studies. In addition, electrografting experiments were performed. The surface characteristics as well as the optoelectronic and redox properties of the annealed films were found to be very similar to hydrogen plasma-treated films. Moreover, the presented method is compatible with atmospheric pressure and provides a low-cost solution to hydrogenate CVD diamond, which makes it interesting for industrial applications. The plausible mechanism for the hydrogen termination of CVD diamond films is based on the formation of surface carbon dangling bonds and carbon-carbon unsaturated bonds at the applied tempera-ture, which react with molecular hydrogen to produce a hydrogen-terminated surface.

On local bifurcations in neural field models with transmission delays.

Neural field models with transmission delays may be cast as abstract delay differential equations (DDE). The theory of dual semigroups (also called sun-star calculus) provides a natural framework for the analysis of a broad class of delay equations, among which DDE. In particular, it may be used advantageously for the investigation of stability and bifurcation of steady states. After introducing the neural field model in its basic functional analytic setting and discussing its spectral properties, we elaborate extensively an example and derive a characteristic equation. Under certain conditions the associated equilibrium may destabilise in a Hopf bifurcation. Furthermore, two Hopf curves may intersect in a double Hopf point in a two-dimensional parameter space. We provide general formulas for the corresponding critical normal form coefficients, evaluate these numerically and interpret the results.

An integration of external information for foreign stallions into the Belgian genetic evaluation for jumping horses.

The aim of this study was to test the integration of external information, i.e. foreign estimated breeding values (EBV) and the associated reliabilities (REL), for stallions into the Belgian genetic evaluation for jumping horses. The Belgian model is a bivariate repeatability Best Linear Unbiased Prediction animal model only based on Belgian performances, while Belgian breeders import horses from neighbouring countries. Hence, use of external information is needed as prior to achieve more accurate EBV. Pedigree and performance data contained 101382 horses and 712212 performances, respectively. After conversion to the Belgian trait, external information of 98 French and 67 Dutch stallions was integrated into the Belgian evaluation. Resulting Belgian rankings of the foreign stallions were more similar to foreign rankings according to the increase of the rank correlations of at least 12%. REL of their EBV were improved of at least 2% on average. External information was partially to totally equivalent to 4 years of contemporary horses' performances or to all the stallions' own performances. All these results showed the interest to integrate external information into the Belgian evaluation.

Myostatin mutation causing double muscling could affect increased psoroptic mange sensitivity in dual purpose Belgian Blue cattle.

Belgian Blue cattle are known for their high degree of muscling and good carcass qualities. This high degree of muscling is mainly caused by a mutation in the myostatin gene (MSTN). Although the MSTN mutation is considered as fixed in the Belgian Blue breed, segregation is occurring in a sub-population bred for dual purpose. In the latter population, we observed an association between the mutation in MSTN and susceptibility to psoroptic mange, a skin disease caused by Psoroptes ovis mites that heavily plagues Belgian Blue cattle. In total, 291 animals were sampled and screened for their susceptibility for mange lesions and their MSTN genotype. Via linear mixed modelling, we observed that homozygous mutant animals had a significant increase in the size of mange lesions (+2.51% lesion extent) compared to homozygous wild type. These findings were confirmed with zero-inflated modelling, an animal model and odds analysis. Risk ratios for developing severe mange lesions were 5.9 times as high for homozygous mutant animals. All analyses confirmed an association between the MSTN genotype and psoroptic mange lesion size.