Skip to main content

Selective breeding can contribute to bovine tuberculosis control and eradication

Abstract

Bovine tuberculosis (bTB) persists in many countries having a significant impact on public health and livestock industry finances. The incidence and prevalence of new cases in parts of the UK and elsewhere over the past decades warrant intensified efforts towards achieving Officially Tuberculosis Free (OTF) status in the respective regions. Genetic selection aiming to identify and remove inherently susceptible animals from breeding has been proposed as an additional measure in ongoing programmes towards controlling the disease. The presence of genetic variation among individual animals in their capacity to respond to Mycobacterium bovis exposure has been documented and heritability estimates of 0.06-0.18 have been reported. Despite their moderate magnitude, these estimates suggest that host resistance to bTB is amenable to improvement with selective breeding. Although relatively slow, genetic progress can be constant, cumulative and permanent, thereby complementing ongoing disease control measures. Importantly, mostly no antagonistic genetic correlations have been found between bTB resistance and other animal traits suggesting that carefully incorporating the former in breeding decisions should not adversely affect bovine productivity. Simulation studies have demonstrated the potential impact of genetic selection on reducing the probability of a breakdown to occur or the duration and severity of a breakdown that has already been declared. Furthermore, research on the bovine genome has identified multiple genomic markers and genes associated with bTB resistance. Nevertheless, the combined outcomes of these studies suggest that host resistance to bTB is a complex, polygenic trait, with no single gene alone explaining the inherent differences between resistant and susceptible animals. Such results support the development of accurate genomic breeding values that duly capture the collective effect of multiple genes to underpin selective breeding programmes. In addition to improving host resistance to bTB, scientists and practitioners have considered the possibility of reducing host infectivity. Ongoing studies have suggested the presence of genetic variation for infectivity and confirmed that bTB eradication would be accelerated if selective breeding considered both host resistance and infectivity traits. In conclusion, research activity on bTB genetics has generated knowledge and insights to support selective breeding as an additional measure towards controlling and eradicating the disease.

Introduction

Tuberculosis is a key disease with multiple pathogen and host species involved in transmission and consistently ranks in the top five most significant diseases worldwide [1]. According to the World Health Organisation’s “End TB strategy” new tuberculosis cases and deaths worldwide must be reduced by 90% and 95%, respectively, by the end of the 2016-2035 period [2]. The ambitiousness of these goals combined with the fact that about one third of the time has already elapsed set the stage for a considerable amount of intensified activity and attention, at multiple levels, to all Mycobacterium infections, including Mycobacterium bovis (M. bovis), the key pathogen agent responsible for bovine tuberculosis.

Bovine tuberculosis (bTB) presents a global problem to the industry. As a zoonotic disease, bTB receives due attention as a public health, financial, and political issue. As such, monitoring, control, and eradication programmes have been established in multiple afflicted countries and regions.

In the United Kingdom (UK), bTB is endemic and costly to the tune of £175-200 million annually [3, 4]. The disease has returned and consolidated its presence in parts of Great Britain and Northern Ireland over the past 30-40 years. The key policy in place at the moment is national surveillance based mostly on animal reaction to the single intradermal comparative cervical tuberculin (SICCT) test, commonly known as the “skin test”, culling of positive animals, post-mortem examination and confirmation, and declaration of a bTB breakdown in the herd where the positive animal(s) was found combined with movement restrictions. Scotland has been the only Officially Tuberculosis Free (OTF) nation in the UK since 2009, while the goal is for England to achieve this status by 2038 followed by the other regions. To achieve OTF status, the percentage of new breakdowns must not exceed 0.1% per year for six consecutive years [5]. According to 2009 data, 5-6% of herds in the non-OTF parts of the UK were affected [6] and the figure did not improve until recently [7]. However, the number of new breakdowns in England in the 12-month period to March 2022 was 2,892, down 8.80% from the previous 12-month period ending in March 2021 [7]. A similar reduction was observed in the number of cattle slaughtered due to a bTB incident in the corresponding time period. An average annual decreasing trend of about 8% in new herd incidents was also observed in the previous four 12-month periods. The challenge now is not only to maintain but, importantly, accelerate this downwards trend over the next decades to ensure the OTF target is met.

Selective breeding has been proposed as an additional possible means to control bTB, thereby contributing to the ongoing monitoring and eradication programmes. The main premise is that individual animals exposed to the same pathogenic conditions responsible for a disease react and respond differently: some animals remain healthy, others become a little sick, others more so, and some actually die. This variation is partly attributed to genetic factors manifested in protein-coding genes and/or regulatory sequence elements that render certain animals inherently more resistant or susceptible to the disease than others. In the case of bTB, identifying animals that carry genetic variants contributing to infection resistance and selecting them for breeding can promote the overall population resistance and contribute to bTB control. In a similar but more conservative scenario, individual animals bearing genetic variants that render them susceptible to the disease are identified and removed from the breeding process, thereby removing said variants from the subsequent generations. The potential of selective breeding in bTB control is reflected in a relatively recent strategy review commissioned by the UK government, led by Sir Charles Godfray of the University of Oxford [8] where a relevant excerpt reads as follows: ”There is evidence that cattle vary genetically in their susceptibility to bovine TB and, once infected, there may be genetic variation in their propensity to spread the infection to other animals. It is thus possible to select for better disease susceptibility and transmission traits as part of cattle breeding programmes.” [8].

The present review discusses the genetic component of variation in host response to bTB infection and demonstrates how this variation may be applied in selective breeding programmes towards enhancing cattle resistance to bTB. Results from relevant studies are summarised showing and quantifying the potential impact on disease control and eradication. Further to disease resistance, the role of genetics in reducing animal infectivity regarding bTB is briefly discussed.

Genetic variation among animals for resistance to bTB

The genetic basis of host resistance and susceptibility to bTB infection was addressed in a 2010 review [9] and potential future benefits were already acknowledged at the time. Multiple studies were subsequently conducted and quantified the genetic variation for bTB in different cattle populations and countries, using a variety of statistical models and trait definitions, and reported heritability estimates that range mostly between 0.06 and 0.18 [10,11,12,13,14,15]. The majority of studies were based on large amounts of field data combined with extensive animal pedigree information. Although different datasets and trait definitions were considered in these studies, all phenotypes related to animal response to bTB infection and reflected their susceptibility or resistance to the disease. The range of reported heritability estimates means that 6-18% of the differences observed among individual animals in their response to the pathogen are attributed to heritable genetic factors. In two studies based on the UK national bTB surveillance data, heritability of bTB resistance was 0.09-0.12 in dairy cattle [14] and 0.09-0.11 in beef cattle (author’s unpublished results, 2020).

More recent studies included animal genomic data based on genotyping with genome-wide DNA arrays in the calculations and reported higher heritability estimates of 0.21-0.27 [16,17,18,19]. The advantage of using such data is that they more accurately capture finer details of the genetic background dictating animal resistance to the disease. The caveat is that the amount of genomic data is usually smaller compared to the size of pedigree data on which conventional genetic studies are based; heritability estimates also tend to be sometimes inflated. Fortunately, animal genotyping costs are decreasing rapidly, and the number of individual animals genotyped has been notably increasing in many breeding programmes worldwide, leading to improved accuracy of the genomic estimates. More discussion on genome-wide data in bTB studies is provided in a later section of this review.

The overall magnitude of heritability reported so far may seem moderately low but for most animal traits of interest estimates normally do not exceed 0.50. Heritability of health-related traits is known to be generally modest to low, usually less than 0.15. Nevertheless, the genetic variation manifested in the above heritability estimates for bTB resistance is enough to achieve progress with selective breeding. The latter has been decisively demonstrated for other cattle health traits with similar or even lower heritability values such as clinical mastitis, lameness, and metabolic diseases, which now feature prominently in selective breeding programmes worldwide [20].

In conclusion, there is confidence that genetic variation in cattle for bTB resistance is indeed sufficient to identify animals that are inherently more resistant than others to M. Bovis infection. Selective breeding programmes would then need appropriate genetic evaluation systems and tools in place to properly identify the inherently resistant and susceptible animals in order to inform breeding decisions that will underpin enhanced resistance in the subsequent generations.

Genetic evaluations for bTB resistance

Research conducted on bTB surveillance data in the UK led to the development of a national genetic evaluation of dairy cattle for bTB resistance [14]. This system estimates the genetic merit or breeding value of each individual dairy bovine reared in the country for resistance to the infection. Supported by the Agriculture and Horticulture Development Board, the national farm levy body, the genetic evaluation project collated relevant data from different national databases, including: bTB records from surveillance in Great Britain and Northern Ireland, animal location and movement data, production history, and extensive animal pedigree information.

The study first combined SICCT test and post-mortem data, namely records on bTB lesions and M. Bovis culture results, to define the health status of each animal, which would constitute their individual phenotype for bTB infection. Furthermore, probabilities of animal infection at different stages of a breakdown were assigned to each animal that eventually tested positive to the SICCT test [14]. These probabilities were based on the sensitivity estimates of the test as reported in the literature. Animals that remained healthy throughout the epidemic were included in the data as controls.

A genetic model was developed to analyse these bTB phenotypes [14]. The model was subsequently expanded and, in its final form, accounts for multiple non-genetic factors that are known to affect the trait, such as the location and duration of the breakdown, the calendar year and month the breakdown started, the herd size, the age and lactation number of the cow, and factors such as heterosis and recombination, since the genetic evaluation combines data from multiple dairy breeds and must account for breed composition and potential crossbreeding.

Further to non-genetic effects, the model includes an animal additive genetic effect that, after data analysis, provides an estimate of the genetic merit of each individual animal for resistance to bTB that accounts for all other sources of systematic variation fitted in the model. Here, although the phenotypic data pertain only to milking cows, breeding values are also estimated for their sires, through the pedigree connections that are included in the data analysis. Actually, the sire breeding values are more accurate than those of individual cows on the account of the former being associated with multiple progeny with records in the system. In a validation exercise, the correlation of the sire estimated breeding values with the proportion of infected milking daughters was -0.68, suggesting that sires with highly positive values of genetic merit are associated with reduced infection rate in their progeny [14]. Similarly, the daughters of the 20 genetically most resistant sires are 21% less likely to become infected that those of the 20 most susceptible sires [14].

Importantly, the correlation between genetic evaluations for bTB resistance and other animal traits is mostly zero, suggesting that selective breeding to enhance animal resistance to the disease would not compromise improvement in the other traits [14]. In fact, moderate positive desirable correlations have been reported of bTB with lifespan and the Profitable Lifetime Index (£PLI), the overall dairy selection index in the UK that combines the economically important animal traits in one single value [14]. This result consolidates the expectation that enhancing the inherent cattle resistance to bTB would have favourable effects on animal longevity, with the associated benefits of enriched animal welfare, reduced waste, and improved environmental footprint.

A subsequent study expanded the bTB genetic evaluations to include genome-wide marker data from animal genotyping using DNA arrays [21]. Different models of genomic analysis were investigated and the single-step Best Linear Unbiased Prediction (BLUP) method [22, 23], which allows the simultaneous combination of all phenotypic, genomic and pedigree data available, resulted in the most accurate estimates of animal genetic merit for bTB resistance [21]. Results from this type of analysis would facilitate the identification and selection of young breeding animals early in life, thereby delivering faster benefits of selective breeding in bTB control and eradication programmes.

Research on the genetics of bTB resistance in dairy cattle has led to the development of the official UK national genetic evaluation, run since 2016 at the Edinburgh Genetic Evaluation Services on behalf of the Agriculture and Horticulture Development Board. The genetic evaluation outcome, the actual selection index, has been termed “TB Advantage” [24] and is now available for dairy sires of all breeds with progeny (milking daughters) in the national bTB surveillance programme. TB Advantage is expressed on a scale normally ranging from -4 to +4, where, for every +1 point in the index, 1% fewer daughters are expected to become infected during a bTB breakdown. This information allows dairy farmers to focus on the most positive or remove the most negative TB Advantage bulls from the list of potential sires selected to artificially inseminate their cows.

Another research project investigated the feasibility of a genetic evaluation for bTB resistance in UK beef cattle. Data analyses were conducted on 437,235 animals, progeny of 36,790 sires from 53 beef breeds. These data pertained to 8,648 bTB breakdowns where 7.3% of the animals were infected (author’s unpublished results, 2020). Significant genetic variation was detected for bTB resistance and was manifested in variation in the resulting estimated breeding values of the sires, thereby demonstrating the feasibility to distinguish breeding animals expected to produce bTB resistant from susceptible progeny. As was the case with dairy breeds, correlations of these breeding values with beef traits such as animal growth, carcass characteristics and lifespan were either negligible or moderately favourable, suggesting that genetic selection to enhance bTB resistance in beef cattle would not compromise the current breeding goal. Incorporation of bTB in the official national genetic evaluation process of beef cattle is imminent.

Genetic evaluations for bTB resistance are also available in the Republic of Ireland [25] and are offered as tools for genetic selection to complement the on-going national bTB control programme towards accelerating the disease eradication process. These evaluations are based on relevant population genetic research conducted on Irish bTB surveillance data [13, 15] that demonstrated the presence of sufficient genetic variation and the feasibility of accurate genetic evaluations for the trait. The actual evaluations are run since 2019 by the Irish Cattle Breeding Federation for both dairy and beef breeds in the country [25]. Estimated breeding values are expressed as the predicted proportion of infected progeny of a sire. In this regard, contrary to the UK genetic evaluation, lower values are desirable. For example, when two sires with estimated breeding values of 5% and 10% are assessed, 5% fewer progeny of the first sire are expected to be diagnosed as infected compared to the second.

Despite being expressed on distinctly different scales in Ireland and the UK, estimated breeding values for resistance to bTB offer essentially the same information, namely, they identify and separate genetically resistant and susceptible individuals. Thus, they have the same potential utility in supporting selective breeding programmes in the respective countries. The key is in the interpretation and proper use.

These genetic evaluations are now well tested and validated nationally, so that combination and a joint analysis of data from the two countries to produce across-country estimated breeding values for bTB become a realistic option. Experience with the international genetic evaluations provided by the International Bull Evaluation Service (Interbull) [26] since 1994 for an array of other bovine traits has been very positive. It is clearly possible to combine data from different countries and populations with good genetic links, account for differences in the national definitions and systems, analyse jointly, and then express the estimated breeding values of all selection candidates available internationally on the national scale used in each country. In addition to enhancing the accuracy of the genetic evaluation and broadening the selection pool for the farmers, outcomes facilitate trade across the border, streamline the use of bovine genetics available globally, and promote collaboration and harmonisation of practices in the participating countries.

Genomic architecture of cattle resistance to bTB

The advent and broad implementation of bovine genome-wide DNA arrays (chips) based on tens or hundreds of thousands of Single Nucleotide Polymorphism (SNP) markers have enabled research that cast closer insights on genomic regions and genes associated with animal traits of interest. The key principle is that these SNPs are within or in high linkage disequilibrium with genes and genomic regulatory regions that control important animal traits. This approach has found widespread use in livestock breeding programmes [27].

Multiple genomic association studies using data from the above mentioned DNA arrays have been conducted on bTB resistance in various bovine populations in different countries [15, 16, 18, 19, 28,29,30,31]. The majority of the studies reported single markers with relatively modest effects on the trait. Furthermore, mainly different markers were found associated with cattle resistance to bTB across these studies, partly reflecting differences among populations, trait definitions, DNA array specifications, and models of statistical analysis. Nevertheless, the overall picture is that of a trait that is not being controlled by only a single or a few genes.

In addition, a few regional heritability mapping studies have been conducted using genome-wide SNP data [16, 18, 32]. Some of these studies revealed genomic regions associated with bTB resistance that included significant SNPs identified in the genomic association, but others were not necessarily consistent with single marker results on the same data. In the latter case, these would be genomic regions that harbour multiple markers and genes, each with a small undetectable effect, but combined with a more substantial and discernible effect on bTB resistance.

Using publicly available bovine genome assemblies, several annotated genes were identified neighbouring SNPs or being located within genomic regions associated with bTB traits in the aforementioned studies. Further studies more closely examined specific candidate genes affecting the host’s resistance or susceptibility to the disease [33,34,35,36,37,38]. These studies offered novel insights into the function and role of the individual genes in networks and pathways that are linked with the animal’s capacity to fight off infection.

Table 1 summarises the genes that have been reportedly associated with cattle resistance to bTB in different genomic association, regional heritability mapping, and candidate gene studies. The collective outcome of these studies is that cattle resistance to bTB is a largely polygenic trait dictated by multiple genes across the bovine genome. Therefore, selective breeding based on genetic and genomic evaluations as previously discussed would be the best approach towards genetic improvement of the trait. Furthermore, identification of single markers and genomic regions with noticeable effects may inform the genetic evaluation process, with due emphasis being placed on these markers and regions resulting in an increase in the accuracy of the estimate of the genetic merit of individual animals for resistance to bTB.

Table 1 Genes associated with host resistance to bovine tuberculosis in different research studies in the past 10 years demonstrate the polygenic nature of the trait

Potential impact of selective breeding on the disease dynamics

Genetic improvement with selective breeding is a long-term process and benefits become visible in future generations. From the moment the genetically most resistant breeding animals are identified and mated, it would take at least one generation to see the effect in their progeny. The process will continue in the subsequent generation and then in the following one and so on, and benefits gradually accumulate. The process is slow, but the effects are permanent and cumulative. The benefits of selective breeding have been observed in multiple livestock traits not least those associated with animal production, health, fertility, and longevity over the past century and have been duly documented [20].

Presence of heritable genetic variation and the feasibility of harnessing it in the estimation of accurate genetic evaluations of animals, as discussed above, may underpin selective breeding practices aiming to enhance resistance to bTB and contribute to the overall control and eventual eradication of the disease. In a generational context, these selection tools have become available relatively recently and the full effects are expected to materialise in due course in the subsequent generations.

In the meantime, simulation studies have been conducted to predict and quantify the long-term effect of selective breeding for enhanced animal resistance or reduced susceptibility to bTB. In one such study, animal data were simulated reflecting the current bTB dynamics in the UK regarding disease prevalence, and size and duration of bTB breakdowns [39]. A model of disease transmission was implemented based on a sequence of stages, where the animal: (i) is at first susceptible but not infected, (ii) is then exposed to the pathogen, (iii) becomes infectious and (iv) is finally detected as infected with the SICCT test and, therefore, removed from the herd. This choice of model was based on empirical data showing that infected cattle may become infectious before entering the detectable stage [40]. The same level of genetic variation among individual animals in bTB resistance found in the real data analyses [14] was introduced in the simulations. Subsequently, various levels of selection emphasis on bTB resistance were examined: in the most intense selection scenario, the 10% genetically most resistant sires were selected to breed the next generation, whereas in the least intense case 70% were selected. The latter reflects a conservative approach where the 30% genetically most susceptible animals are removed from breeding. A baseline scenario of no selection was also examined for comparison. Twenty generations of selection were simulated, roughly corresponding to 80-100 years in a dairy cattle conventional breeding programme.

Results showed that any level of selective breeding would be associated with a faster ending of the epidemic compared to the scenario of no selection [39]. As expected, higher selection intensities had a stronger impact on the epidemic profile. The probability of a breakdown to occur would half after about 5 generations of intense selection of the top genetic merit sires. A more conservative approach based on avoiding the genetically most susceptible sires in breeding decisions would take about 12 generations to reduce the probability of a breakdown to one half. When a breakdown actually occurs, selective breeding would decrease its duration and severity manifested in the number of infected animals after the breakdown had been declared. After one generation of selection, duration and severity could decrease by 5-10% and 11-17%, respectively, dependent on the intensity of selection [39]. In practice, selection programmes may be adapted to the geographical risk levels, with actively intense selection applied in high-risk areas and moderate conservative practices elsewhere.

The above research was based on a model reflecting the “worst-case” scenario where an animal becomes infectious and starts transmitting the disease before it can be detected as infected and removed from the herd. Although this assertion is consistent with the suggestion that cattle with bTB lesions in the respiratory tract should be considered potential transmitters of the pathogen and sources of infection for other animals in the herd [41], previous non-genetic epidemiological studies on bTB typically assumed that infected animals usually react positively to the SICCT test and are removed from the herd before they become infectious [42, 43]. The two models are discussed in detail in [39]. Further genetic epidemiological studies reversing the order of the two stages so that infected animals were detectable before they became infectious demonstrated that the benefits of selective breeding were proportionally still the same compared to the no-selection option [39].

Key conclusions drawn for this research were that genetic selection could accelerate bTB eradication, as a complementary measure to already existing ones. Thus, selective breeding for enhanced animal resistance to bTB can contribute to herds and regions achieving their OTF goals faster.

Genetic background of bTB infectivity

Further to improving animal resistance to bTB, a disease eradication programme would also benefit if the rate of transmission from one animal to another decreased. There is increasing evidence, albeit mostly empirical, that individual animals differ in their capacity to infect their herd-mates. The presence of inherently super-spreader individuals is now widely accepted [44] and identifying them early and removing from the breeding process would yield additional benefits. Selective breeding practices could then be expanded to harness this novel host genetic variation and contribute to decreasing infection rates.

Relevant genetic-epidemiological simulation studies have suggested that selective breeding for reduced infectivity together with enhanced resistance to bTB would lead to faster disease eradication than selecting only for enhanced resistance or, of course, not selecting at all [45, 46]. These simulation studies showed that the magnitude of the benefit would vary depending on selection intensities, the heritability of infectivity, and the genetic correlation between infectivity and resistance. Several selective breeding scenarios were assessed in these studies, and in all cases the epidemic risk decreased dramatically, when infectivity and resistance were considered together. For example, in one scenario assuming moderate selection intensity, trait heritability and prediction accuracy, the epidemic risk after five generations of selection decreased by 50% when selective breeding decisions considered infectivity [45]. In another scenario with similar assumed parameters, it took seven generations for the bTB epidemic to die out with selective breeding for reduced susceptibility and three generations when selective breeding aimed at reducing both animal susceptibility and infectivity [46].

However, the genetic basis of infectivity is still largely unknown. Defining and characterising the proper animal phenotype is a serious challenge. Work is currently being conducted on this topic and a recent study on dairy cattle data reported significant genetic variation among index cases regarding the number of secondary cases attributed to them [47]. Here an index case would be the first animal tested and confirmed bTB positive, thereby instigating the breakdown, and secondary cases would be the animals that tested positive in the first stages of the breakdown and assumed to have been infected by the index case. Thus, the trait may be considered as a proxy for animal infectiousness. Different models of analysis were considered and preliminary heritability estimates of 0.09-0.21 derived were of similar magnitude as the heritability estimates for bTB resistance discussed previously. Continuing work addresses refinements in the definition of animal infectivity and modelling improvements. Although more research is definitely warranted to properly define and dissect the genetic basis of infectivity and understand the potential for genetic improvement, these first results are encouraging.

Novel bTB phenotypes

Availability of informative individual animal phenotypes for the traits of interest is the key and most crucial component of modern selective breeding programmes. This entails data of high quality that reflect the true biology of the animal as closely as possible, available in sufficient quantities to allow accurate calculations. Without said phenotypes, even the most advanced genomic technologies and the most elaborate statistical analysis models would have very little opportunity to demonstrate their utility.

In the case of bTB, SICCT test records combined with post-mortem visible lesions and M. bovis culture results constitute the phenotypes in current national genetic evaluations [24, 25]. Thanks to the ongoing national surveillance programmes, such data span multiple decades and cover all affected geographic regions in a country. The depth and amount of data in combination with the high specificity of the SICCT test render these phenotypes and resulting genetic evaluations valuable for selective breeding purposes.

Other diagnostics that may provide new animal phenotypes with improved sensitivity are worth consideration. Gamma-interferon assays have been long used for bTB diagnosis [48] often in combination with the SICCT test [49], especially in high-risk geographic regions. Nevertheless, the development of individual animal phenotypes for bTB based on gamma-interferon to complement the existing ones based on SICCT tests and post-mortem examination results and use in bovine selective breeding programmes has not been explored.

Alternative phenotyping approaches have been proposed in livestock improvement based on molecular traces in the milk of dairy cows that emanate from various animal physiological processes [50]. Mid-infrared (MIR) spectroscopy of milk samples may capture these traces in an effective, non-invasive, and non-interfering process. This has led to the accumulation of large amounts of MIR data in routine milk recording systems in many countries. Analysis of MIR data already provides useful proxy phenotypes for different animal traits associated with milk composition [51], cow body energy [52], pregnancy diagnosis [53], greenhouse gas emissions [54], and feed intake [55].

In a recent study, extensive cow milk MIR data were combined with bTB national surveillance records and machine learning was deployed to jointly analyse the data and deliver a prediction pipeline of the bTB status of individual animals from their MIR spectral profiles [56]. After optimisation of the training process, including data and network, the specificity, sensitivity and accuracy of predictions reached 0.94, 0.96 and 0.95, respectively [56]. These results lend confidence in the method to predict future bTB status of individual animals, thereby offering farmers a useful tool for early management decisions on cows likely to become infected. Furthermore, these predictions may constitute additional novel animal phenotypes for bTB, warranting genetic studies to determine their suitability for genetic evaluations and improvement with selective breeding.

Closing remarks and outlook

Bovine tuberculosis is a complex issue, and an array of factors may determine whether an animal becomes infected, and a breakdown is declared. Complex problems require complex strategies and solutions, and indeed this is the case with ongoing bTB control and eradication programmes in different countries. Genetic selection of inherently resistant individuals that will breed enhanced natural resistance into the subsequent generations is proposed as an additional component within these programmes.

Further improvements in genomic evaluations are in scope with more targeted genotyping both in terms of animals (bulls and cows) and DNA arrays, and incorporation of identified genomic regions with significant effects on the trait in the calculation of breeding values. Decreasing DNA array genotyping costs, genotype imputation improvements, and new approaches such as genotyping by low pass sequencing combined with imputation provide useful opportunities in this regard.

As previously mentioned, the polygenicity of host resistance to M. bovis infection is now well accepted. Future genomic studies of DNA array and whole genome sequence data combined with post-genomic (transcriptional and translational) research will cast further insights into the complex architecture of the trait, including focus on the location and function of both protein-coding genes and regulatory elements [57]. Improvements in the functional annotation of the bovine genome [58] certainly facilitate such research activities, leading to enhanced precision and biological relevance of the genomic evaluation and selective breeding.

Beyond host resistance to M. bovis infection, understanding the genetics of infectivity at the individual animal level is only now starting and we have some way to go before being able to usefully inform genetic selection programmes aiming to decrease the rate of animal-to-animal transmission. Nevertheless, this area of scientific research and development is topical and warrants continuing attention and effort. Future bTB control programmes should be able to combine selective breeding for increased host resistance with reduced animal infectivity.

Furthermore, as has been experienced with the genetics of just about any important animal trait, multi-actor collaboration is a prerequisite to the successful development and implementation of selective breeding programmes. There is scope for continuing interdisciplinary research to understand why and how individual genes affect host traits and how host genetics interact with pathogen variation. Moreover, there are opportunities for across-country cooperation to facilitate international trade that is of key importance to the bovine industry. There are national genetic evaluation systems for bTB in place in at least the UK and Ireland and it is possible to pull data together and produce international genetic evaluations, as is the case with several dozens of other animal traits at the moment [26].

As already mentioned, selective breeding is a long-term process, it takes time to see the benefits, but it can work. Of course, this requires faith and good judgement, and everybody to be on board, not least farmers, practitioners, breeder associations, breeding companies, scientists, governments, and policy makers. Only then would the benefits discussed here truly materialise, contributing decisively to the control and eventual eradication of the disease.

Availability of data and materials

Not applicable.

Abbreviations

bTB:

Bovine Tuberculosis

OTF:

Officially Tuberculosis Free

M. bovis:

Mycobacterium bovis

UK:

United Kingdom

SICCT:

Single Intradermal Comparative Cervical Tuberculin

£PLI:

Profitable Lifetime Index

BLUP:

Best Linear Unbiased Prediction

Interbull:

International Bull Evaluation Service

SNP:

Single Nucleotide Polymorphism

MIR:

Mid-infrared

References

  1. IHME Global Burden of Disease; https://www.healthdata.org/gbd/2019, Last Access 17/10/2022.

  2. World Health Organisation; https://www.who.int/teams/global-tuberculosis-programme/the-end-tb-strategy, Last Access 17/10/2022.

  3. Abernethy DA, Upton P, Higgins IM, McGrath G, Goodchild AV, Rolfe SJ, Broughan JM, Downs SH, Clifton-Hadley R, Menzies FD, de la Rua-Domenech R, Blissitt MJ, Duignan A, More SJ. Bovine tuberculosis trends in the UK and the Republic of Ireland, 1995–2010. Vet Rec. 2013;172:312. https://doi.org/10.1136/vr.100969.

    Article  CAS  PubMed  Google Scholar 

  4. Allen AR, Skuce RA, Byrne AW. Bovine tuberculosis in Britain and Ireland – A perfect storm? the confluence of potential ecological and epidemiological impediments to controlling a chronic infectious disease. Front Vet Sci. 2018;5:109. https://doi.org/10.3389/fvets.2018.00109.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Reviriego Gordejo FJ, Vermeersch JP. Towards eradication of bovine tuberculosis in the European Union. Vet Microbiol. 2006;112(2–4):101–9. https://doi.org/10.1016/j.vetmic.2005.11.034.

    Article  CAS  PubMed  Google Scholar 

  6. European Food Safety Authority and European Centre for Disease Prevention and Control. Scientific Report of EFSA and CDC - The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2009. EFSA J. 2011;9(3):2090.

    Google Scholar 

  7. Department for Environment, Food and Rural Affairs; https://www.gov.uk/government/collections/bovine-tb, Last Access 11/10/2022.

  8. Department for Environment, Food and Rural Affairs; https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/756942/tb-review-final-report-corrected.pdf, Last Access 10/10/2022.

  9. Allen AR, Minozzi G, Glass EJ, Skuce RA, McDowell SWJ, Woolliams JA, Bishop SC. Bovine tuberculosis: the genetic basis of host susceptibility. Proc R Soc B. 2010;277:2737–45. https://doi.org/10.1098/rspb.2010.0830.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Petukhov V, Kochnev N, Panov B, Korotkevich O, Kulikova S, Marenkov V. Genetics of cattle resistance to tuberculosis. Proc 6th World Congress on Genetics Applied to Livestock Production. 1998.

    Google Scholar 

  11. Bermingham ML, More SJ, Good M, Cromie AR, Higgins IM, Brotherstone S, Berry DP. Genetics of tuberculosis in Irish Holstein-Friesian dairy herds. J Dairy Sci. 2009;92(7):3447–56.

    Article  CAS  PubMed  Google Scholar 

  12. Brotherstone S, White I, Coffey M, Downs S, Mitchell A, Clifton-Hadley R, More S, Good M, Woolliams J. Evidence of genetic resistance of cattle to infection with Mycobacterium bovis. J Dairy Sci. 2010;93(3):1234–42.

    Article  CAS  PubMed  Google Scholar 

  13. Richardson IW, Bradley DG, Higgins IM, More SJ, McClure J, Berry DP. Variance components for susceptibility to Mycobacterium bovis infection in dairy and beef cattle. Genet Sel Evol. 2014;46:77. https://doi.org/10.1186/s12711-014-0077-1.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Banos G, Winters M, Mrode R, Mitchell A, Bishop SC, Woolliams JA, Coffey MP. Genetic evaluation for bovine tuberculosis resistance in dairy cattle. J Dairy Sci. 2016;100(2):1272–81.

    Article  PubMed  Google Scholar 

  15. Ring SC, Purfield DC, Good M, Breslin P, Ryan E, Blom A, Evans RD, Doherty ML, Bradley DG, Berry DP. Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data. PLoS One. 2019. https://doi.org/10.1371/journal.pone.0212067.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Bermingham ML, Bishop SC, Woolliams JA, Pong-Wong R, Allen AR, McBride SH, Ryder JJ, Wright DM, Skuce RA McDowell SWJ, Glass EJ. Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis. Heredity. 2014;112(5):543–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Tsairidou S, Woolliams JA, Allen AR, Skuce RA, McBride SH, Wright DM, Birmingham ML, Pong-Wong R, Matika O, McDowell SWJ, Glass EJ, Bishop SC. Genomic prediction for tuberculosis resistance in dairy cattle. PLoS One. 2014;9(5):e96728.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Raphaka K, Matika O, Sánchez-Molano E, Mrode R, Coffey MP, Riggio V, Glass EJ, Woolliams JA, Bishop SC, Banos G. Genomic regions underlying susceptibility to bovine tuberculosis in Holstein-Friesian cattle. BMC Genet. 2017;18:27. https://doi.org/10.1186/s12863-017-0493-7.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Callaby R, Kelly R, Mazeri S, Egbe F, Benedictus L, Clark E, Doeschl-Wilson A, Bronsvoort B, Salavati M, Muwonge A. Genetic diversity of Cameroon cattle and a putative genomic map for resistance to bovine tuberculosis. Front Genet. 2020; 11 https://doi.org/10.3389/fgene.2020.550215.

  20. Miglior F, Fleming A, Malchiodi F, Brito LF, Martin P, Baes CF. A 100-Year Review: Identification and genetic selection of economically important traits in dairy cattle. J Dairy Sci. 2017;100(2):10251–71. https://doi.org/10.3168/jds.2017-12968.

    Article  CAS  PubMed  Google Scholar 

  21. Mrode R, Banos G, Winters M, Coffey MP. Application of various models for the genomic evaluation of bovine tuberculosis in dairy cattle. Interbull Bulletin 2016; 50 https://journal.interbull.org/index.php/ib/article/view/1634/1636.

  22. Aguilar I, Misztal I, Johnson DL, Legarra A, Tsuruta S, Lawlor TJ. Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J Dairy Sci. 2010;93:743–52.

    Article  CAS  PubMed  Google Scholar 

  23. Christensen OF, Lund MS. Genomic prediction when some animals are not genotyped. Genet Sel Evol. 2010;42:2.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Agriculture and Horticulture Development Board https://ahdb.org.uk/knowledge-library/tb-advantage-genetic-index-factsheet, Last Access 18/10/2022.

  25. Irish Cattle Breeding Federation https://www.icbf.com/?page_id=12693, Last Access 19/10/2022.

  26. International Bull Evaluation Service https://interbull.org/index, Last Access 14/10/2022.

  27. Meuwissen THE, Hayes BJ, Goddard ME. Prediction of Total Genetic Value Using Genome-Wide Dense Marker Maps. Genetics. 2001;15(4):1819–29. https://doi.org/10.1093/genetics/157.4.1819.

    Article  Google Scholar 

  28. Finlay EK, Berry DP, Wickham B, Gormley EP, Bradley DG. A genome wide association scan of bovine tuberculosis susceptibility in Holstein-Friesian dairy cattle. PLoS One. 2012;7(2):e30545. https://doi.org/10.1371/journal.pone.0030545.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kassahun Y, Mattiangeli V, Ameni G, Hailu E, Aseffa A, Young DB, Hewinson RG, Vordermeier HM, Bradley DG. Admixture mapping of tuberculosis and pigmentation-related traits in an African-European hybrid cattle population. Front Genet. 2015;6:210. https://doi.org/10.3389/fgene.2015.00210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Richardson IW, Berry DP, Wiencko HL, Higgins IM, More SJ, McClure J, Lynn DL, Bradley DG. A genome-wide association study for genetic susceptibility to Mycobacterium bovis infection in dairy cattle identifies a susceptibility QTL on chromosome 23. Genet Sel Evol. 2016;48:19. https://doi.org/10.1186/s12711-016-0197-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. González-Ruiz S, Strillacci MG, Durán-Aguilar M, Cantó-Alarcón GJ, Herrera-Rodríguez SE, Bagnato A, Guzmán LF, Milián-Suazo F, Román-Ponce SI. Genome-wide association study in Mexican Holstein cattle reveals novel quantitative trait loci regions and confirms mapped loci for resistance to bovine tuberculosis. Animals. 2019;9:636. https://doi.org/10.3390/ani9090636.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Wilkinson S, Bishop SC, Allen AR, McBride SH, Skuce RA, Bermingham ML, Woolliams JA, Glass EJ. Fine-mapping host genetic variation underlying outcomes to Mycobacterium bovis infection in dairy cows. BMC Genomics. 2017;18:477. https://doi.org/10.1186/s12864-017-3836-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Cheng Y, Huang C, Tsai HJ. Relationship of bovine SLC11A1 (Formerly NRAMP1) polymorphisms to the risk of bovine tuberculosis in Holstein cattle. J Vet Sci Technol. 2015;6:5. https://doi.org/10.4172/2157-7579.1000247.

    Article  CAS  Google Scholar 

  34. Cheng Y, Huang C, Tsai HJ. Relationship of bovine TNF-α gene polymorphisms with the risk of bovine tuberculosis in Holstein cattle. J Vet Med Sci. 2016;78(5):727–32. https://doi.org/10.1292/jvms.15-0506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Liu K, Zhang B, Teng Z, Wang Y, Dong G, Xu C, Qin B, Song C, Chai J, Li Y, Shi X, Shu X, Zhang Y. Association between SLC11A1 (NRAMP1) polymorphisms and susceptibility to tuberculosis in Chinese Holstein cattle. Tuberculosis (Edinb). 2017;103:10–5. https://doi.org/10.1016/j.tube.2016.11.003.

    Article  CAS  PubMed  Google Scholar 

  36. Bhaladhare A, Chauhan A, Sonwane A, Kumar A, Singh RV, Prakash C, Kumar S, Kumar P, Kumar S, Bhushan B. Single nucleotide polymorphism in TLR1 and TNFα genes and their association with susceptibility to bovine tuberculosis. Indian J Anim Res. 2020;54(7):805–12. https://doi.org/10.18805/ijar.B-3831.

    Article  Google Scholar 

  37. Holder A, Garty R, Elder C, Mesnard P, Laquerbe C, Bartens MC, Salavati M, Shabbir MZ, Tzelos T, Connelly T, Villarreal-Ramos B. Werling D (2020) Analysis of genetic variation in the bovine SLC11A1 gene, its influence on the Expression of NRAMP1 and potential association with resistance to bovine tuberculosis. Front Microbiol. 2020;11:1420. https://doi.org/10.3389/fmicb.2020.01420.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Chai J, Wang Q, Qin B, Wang S, Wang Y, Shahid M, Liu K, Zhang Y, Qu W. Association of NOS2A gene polymorphisms with susceptibility to bovine tuberculosis in Chinese Holstein cattle. PLoS One. 2021;16(6):e0253339. https://doi.org/10.1371/journal.pone.0253339.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Raphaka K, Sánchez-Molano E, Tsairidou S, Anacleto O, Glass EJ, Woolliams JA, Doeschl-Wilson A, Banos G. Impact of genetic selection for increased cattle resistance to bovine tuberculosis on disease transmission dynamics. Front Vet Sci. 2018;5:237. https://doi.org/10.3389/fvets.2018.00237.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Conlan AJ, McKinley TJ, Karolemeas K, Pollock EB, Goodchild AV, Mitchell AP, Birch CPD, Clifton-Hadley RS, Wood JL. Estimating the hidden burden of bovine tuberculosis in Great Britain. PLoS Comput Biol. 2012;8:e1002730. https://doi.org/10.1371/journal.pcbi.1002730.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. McIlroy S, Neill S, McCracken R. Pulmonary lesions and Mycobacterium bovis excretion from the respiratory tract of tuberculin reacting cattle. Vet Rec. 1986;118:718–21. https://doi.org/10.1136/vr.118.26.718.

    Article  CAS  PubMed  Google Scholar 

  42. Fischer E, Van Roermund H, Hemerik L, Van Asseldonk M, De Jong M. Evaluation of surveillance strategies for bovine tuberculosis (Mycobacterium bovis) using an individual based epidemiological model. Prev Vet Med. 2005;67:283–301. https://doi.org/10.1016/j.prevetmed.2004.12.002.

    Article  CAS  PubMed  Google Scholar 

  43. Brooks-Pollock E, Roberts GO, Keeling MJ. A dynamic model of bovine tuberculosis spread and control in Great Britain. Nature. 2014;511:228–31. https://doi.org/10.1038/nature13529.

    Article  CAS  PubMed  Google Scholar 

  44. Fielding HR, McKinley TJ, Delahay RJ, Silk MJ, McDonald RA. Characterization of potential superspreader farms for bovine tuberculosis: A review. Vet Med Sci. 2021;7:310–21. https://doi.org/10.1002/vms3.358.

    Article  CAS  PubMed  Google Scholar 

  45. Tsairidou S, Allen A, Banos G, Coffey MP, Anacleto O, Byrne AW, Skuce RA, Glass EJ, Woolliams JA, Doeschl-Wilson A. Can We Breed Cattle for Lower bovine TB Infectivity? Front Vet Sci. 2018;5:310. https://doi.org/10.3389/fvets.2018.00310.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Tsairidou S, Anacleto O, Woolliams JA, Doeschl-Wilson A. Enhancing genetic disease control by selecting for lower host infectivity and susceptibility. Heredity. 2019;122:742–58. https://doi.org/10.1038/s41437-018-0176-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Sánchez-Molano E, Madenci D, Tsairidou S, Winters M, Mitchell AP, Banos G, Doeschl-Wilson A. Detection of genetic variability in cattle infectivity for bovine tuberculosis. Proc 12th World Congress on Genetics Applied to Livestock Production 2022; https://www.wageningenacademic.com/pb-assets/wagen/WCGALP2022/57_014.pdf, Last Access 14/10/2022.

  48. Rothel JS, Jones SL, Corner LA, Cox JC, Wood PR. The gamma-interferon assay for diagnosis of bovine tuberculosis in cattle: conditions affecting the production of gamma-interferon in whole blood culture. Australian Vet J. 1992;69:1–4. https://doi.org/10.1111/j.1751-0813.1992.tb09848.x.

    Article  CAS  Google Scholar 

  49. Ryan TJ, Buddle BM, de Lisle GW. An evaluation of the gamma interferon test for detecting bovine tuberculosis in cattle 8 to 28 days after tuberculin skin testing. Research in Veterinary Sci. 2000;69(1):57–61. https://doi.org/10.1053/rvsc.2000.0386.

    Article  CAS  Google Scholar 

  50. Soyeurt H, Dardenne P, Dehareng F, Lognay G, Veselko D, Marlier M, Bertozzi C, Mayeres P, Gengler N. Estimating fatty acid content in cow milk using mid-infrared spectrometry. J Dairy Sci. 2006;89(9):3690–5. https://doi.org/10.3168/jds.S0022-0302(06)72409-2.

    Article  CAS  PubMed  Google Scholar 

  51. Soyeurt H, Dehareng F, Gengler N, McParland S, Wall E, Berry DP, Coffey MP, Dardenne P. Mid-infrared prediction of bovine milk fatty acids across multiple breeds, production systems, and countries. J Dairy Sci. 2011;94(4):1657–67. https://doi.org/10.3168/jds.2010-3408.

    Article  CAS  PubMed  Google Scholar 

  52. McParland S, Banos G, Wall E, Coffey MP, Soyeurt H, Veerkamp RF, Berry DP. The use of mid-infrared spectrometry to predict body energy status of Holstein cows. J Dairy Sci. 2011;94(7):3651–61. https://doi.org/10.3168/jds.2010-3965.

    Article  CAS  PubMed  Google Scholar 

  53. Brand W, Wells AT, Smith SL, Denholm SJ, Wall E, Coffey MP. Predicting pregnancy status from mid-infrared spectroscopy in dairy cow milk using deep learning. J Dairy Sci. 2021;104(4):4980–90. https://doi.org/10.3168/jds.2020-18367.

    Article  CAS  PubMed  Google Scholar 

  54. Dehareng F, Delfosse C, Froidmont E, Soyeurt H, Martin C, Gengler N, Vanlierde A, Dardenne P. Potential use of milk mid-infrared spectra to predict individual methane emission of dairy cows. Animal. 2012;6(10):1694–701. https://doi.org/10.1017/S1751731112000456.

    Article  CAS  PubMed  Google Scholar 

  55. Wallén SE, Prestløkken E, Meuwissen THE, McParland S, Berry DP. Milk mid-infrared spectral data as a tool to predict feed intake in lactating Norwegian Red dairy cows. J Dairy Sci. 2018;101(7):6232–43. https://doi.org/10.3168/jds.2017-13874.

    Article  CAS  PubMed  Google Scholar 

  56. Denholm SJ, Brand W, Mitchell AP, Wells AT, Krzyzelewski T, Smith SL, Wall E, Coffey MP. Predicting bovine tuberculosis status of dairy cows from mid-infrared spectral data of milk using deep learning. J Dairy Sci. 2020;103(10):9355–67. https://doi.org/10.3168/jds.2020-18328.

    Article  CAS  PubMed  Google Scholar 

  57. Goddard ME, Kemper KE, MacLeod IM, Chamberlain AJ, Hayes BJ. Genetics of complex traits: prediction of phenotype, identification of causal polymorphisms and genetic architecture. Proc Biol Sci. 1835;2016(283):20160569. https://doi.org/10.1098/rspb.2016.0569.

    Article  CAS  Google Scholar 

  58. Clark EL, Archibald AL, Daetwyler HD, Groenen MAM, Harrison PW, Houston RD, Kühn C, Lien S, Macqueen DJ, Reecy JM, Robledo D, Watson M, Tuggle CK, Giuffra E. From FAANG to fork: application of highly annotated genomes to improve farmed animal production. Genome Biol. 2020;21:285. https://doi.org/10.1186/s13059-020-02197-8.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The journal editor and reviewers for useful comments and suggestions. The organisers of the Seventh International Conference on Mycobacterium bovis for their invitation that led to this publication. Colleagues and co-authors at Scotland’s Rural College and the Roslin Institute for collaboration on bTB research. The Agriculture and Horticulture Development Board, UK, and the Animal and Plant Health Agency, UK, for continuous support of bTB research.

About this supplement

This article has been published as part of Irish Veterinary Journal Volume 76 Supplement 1, 2023: Seventh International Conference on Mycobacterium bovis. The full contents of the supplement are available online at https://irishvetjournal.biomedcentral.com/articles/supplements/volume-76-supplement-1.

Funding

The author acknowledges the financial support from Scotland’s Rural College for conducting this review.

Author information

Authors and Affiliations

Authors

Contributions

The author conducted the review and wrote and revised the manuscript.

Corresponding author

Correspondence to Georgios Banos.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The author declares no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Banos, G. Selective breeding can contribute to bovine tuberculosis control and eradication. Ir Vet J 76 (Suppl 1), 19 (2023). https://doi.org/10.1186/s13620-023-00250-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13620-023-00250-z

Keywords