Authoritative subspecies diagnosis tool for European honey bees based on ancestry informative SNPs.

BACKGROUND: With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and FST) to select the most informative SNPs for ancestry inference. RESULTS: Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof. CONCLUSIONS: The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

Genetic structure of Iranian indigenous sheep breeds: insights for conservation.

Genetic structure and genetic diversity levels of indigenous Iranian sheep breeds are not clear, despite the interest this region has in itself as an important center for domestication of livestock. Early population genetic studies have reported high levels of diversity among Iranian sheep breeds until recently, when high admixture levels and genetic homogeneity have been detected. The rapid reduction of diversity observed in Iranian breeds might be due to an increasing trend of intensive crossbreeding practices or even total replacement of native breeds by highly specialized and productive ones. From a conservative perspective, this situation is highly concerning; thus, it might be wise to consider a conservation program in Iran to preserve the original genetic diversity in native sheep breeds. In this study, a total of 1065 animals with the purest morphological features representing 24 Iranian indigenous sheep breeds were sampled, corresponding to ancestral breed diversity. These samples were genotyped for 17 microsatellite loci in order to (1) determine the native ancestral diversity of Iranian breeds, (2) define the degree of genetic relationship among studied breeds, and (3) assess conservation priorities among defined groups. Our results showed no recent loss of diversity, but high genetic diversity levels for indigenous sheep breeds in Iran. Indeed, the analysis of conservation priorities pointed out the importance of 8 breeds for maintaining Iranian sheep breeds' maximum genetic diversity. Thus, under a genetic perspective, these 8 breeds should be the ones included into conservation programs for restocking endangered areas.

Genomic selection signatures in sheep from the Western Pyrenees.

BACKGROUND: The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades. RESULTS: Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index Fst. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program. CONCLUSIONS: Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

Clupeiformes' Egg Envelope Proteins characterization: The case of Engraulis encrasicolus as a proxy for stock assessment through a novel molecular tool.

Zona radiata proteins are essential for ensuring bactericidal resistance, oocyte nutrients uptake and functional buoyancy, sperm binding and guidance to the micropyle, and protection to the growing oocyte or embryo from the physical environment. Such glycoproteins have been characterized in terms of molecular structure, protein composition and phylogenetics in several chordate models. Nevertheless, research on teleost has not been extensive. In Clupeiformes, one of the most biologically relevant and commercially important order which accounts for over 400 species and totally contributes to more than a quarter of the world fish catch, Egg Envelope Protein (EEP) information exist only for the Clupea pallasii and Engraulis japonicus species. The European anchovy, Engraulis encrasicolus, the target of a well-consolidated fishery in the Mediterranean Sea, has been ignored until now and the interest on the Otocephala superorder has been fragmentally limited to some Cypriniformes and Gonorynchiformes, as well. The aim of the present study was to fill the ZP protein-wise gap of knowledge afflicting the understanding of the European anchovy's reproductive process and to expand the background on Clupeiformes. We cloned the five Engraulis encrasicolus' zp genes and deduced their products, determined their tissue distribution, quantified their mRNA expression throughout the reproductive cycle and provided an insight into their evolution through phylogenetic tools. Furthermore, we proposed a multivariate statistics-based method to objectively infer and/or confirm the classification of Engraulis encrasicolus' sexual maturity stages by analyzing data of zp mRNAs' relative abundance.

Genetic variants associated with rheumatoid arthritis patients and serotypes in European populations.

OBJECTIVES: To replicate the association of rheumatoid arthritis (RA) susceptibility loci in an independent European sample and to assess their specificity with anti-citrullinated protein antibodies (ACPA) status. METHODS: A selection of 64 SNP previously associated with RA have been typed in a cohort of 267 RA patients (169 ACPA-positive and 98 ACPA-negative) and 152 controls from the Rheumatology Units of the University Hospital of Pisa (Italy) and the University of Pécs Medical Center (Hungary). Regression analyses were performed first considering overall RA patients and secondly, taking both serotype subgroups as different disease entities. The results have been adjusted for age, gender and origin of individuals. RESULTS: The well-known CD2, REL, TNFAIP3, IRF5, PTPRC, and CCR6 have been confirmed as RA disease associated loci together with recently discovered BACH2, RASGRP1, and IKZF3 loci, taking all RA patients as a unique phenotype. Results from both serological subgroups separately reflect the specificity of these susceptibility loci and show additional ACPA-positive specific associations for variants at IL6R, IL2RA, BLK, DDX6, IL6, and TLE3 genes. CONCLUSIONS: The results from GAPAID project are consistent with previously established RA disease associations for CD2, PTPRC, REL, CCR6, TNFAIP3, IRF5, BLK, IL2RA, and DDX6 loci. In addition, IL6R, BACH2, RASGRP1, TLE3, and IKZF3 are replicated for the first time in an independent European population and IL6 appears to be a suggestive new RA associated locus. The stratified analysis based on ACPA status provides further support for distinct genetic aetiologies of RA subsets, which might have therapeutic implications.

Influence of MTHFR C677T polymorphism on methotrexate monotherapy discontinuation in rheumatoid arthritis patients: results from the GAPAID European project.

OBJECTIVES: Methotrexate (MTX) is the most widely prescribed drug for rheumatoid arthritis (RA) patients, but 45% of them discontinue therapy within two years, either due to inefficacy or toxicity. Several authors have reported contradictory results related to C677T polymorphism in the MTHFR gene and response to MTX in RA. The purpose of this study was to further explore this genotype-response association in a European RA population. METHODS: This retrospective longitudinal study included a total of 269 RA patients from Italy and Hungary, of whom 73.2% had available data on MTX treatment (197 patients). C677T polymorphism (rs1801133) was genotyped by quantitative PCR using TaqMan assays. Genotype association analysis and Kaplan-Meier method were used for statistical comparisons between patients continuing and patients who abandoned MTX treatment. RESULTS: A total of 85 out of the 197 RA patients (43%) abandoned MTX treatment by the time of analysis. No significant genotype-MTX discontinuation association was found for the overall population, either at the end of the study (p=0.375), or during the follow-up (p=0.324). When the analysis was restricted to the 68 patients on MTX monotherapy, a borderline association (OR 3.15, 95% CI 0.93-10.67, p=0.057) was noted with the recessive genetic model. In agreement with that, a Kaplan-Meier analysis showed a significantly shorter time-to-discontinuation of MTX monotherapy for homozygous carriers of the T-allele (p=0.042). CONCLUSIONS: These results demonstrate that the C677T polymorphism in the MTHFR gene is involved in MTX monotherapy discontinuation in a multicentre European patient cohort, confirming previous results.

Stress Response to Climate Change and Postharvest Handling in Two Differently Pigmented Lettuce Genotypes: Impact on Alternaria alternata Invasion and Mycotoxin Production.

Many species of Alternaria are important pathogens that cause plant diseases and postharvest rots. They lead to significant economic losses in agriculture and affect human and animal health due to their capacity to produce mycotoxins. Therefore, it is necessary to study the factors that can result in an increase in A. alternata. In this study, we discuss the mechanism by which phenol content protects from A. alternata, since the red oak leaf cultivar (containing higher phenols) showed lower invasion than the green one, Batavia, and no mycotoxin production. A climate change scenario enhanced fungal growth in the most susceptible cultivar, green lettuce, likely because elevated temperature and CO2 levels decrease plant N content, modifying the C/N ratio. Finally, while the abundance of the fungi was maintained at similar levels after keeping the lettuces for four days at 4 °C, this postharvest handling triggered TeA and TEN mycotoxin synthesis, but only in the green cultivar. Therefore, the results demonstrated that invasion and mycotoxin production are cultivar- and temperature-dependent. Further research should be directed to search for resistant cultivars and effective postharvest strategies to reduce the toxicological risk and economic losses related to this fungus, which are expected to increase in a climate change scenario.

Recovering High-Quality Host Genomes from Gut Metagenomic Data through Genotype Imputation.

Metagenomic datasets of host-associated microbial communities often contain host DNA that is usually discarded because the amount of data is too low for accurate host genetic analyses. However, genotype imputation can be employed to reconstruct host genotypes if a reference panel is available. Here, the performance of a two-step strategy is tested to impute genotypes from four types of reference panels built using different strategies to low-depth host genome data (≈2× coverage) recovered from intestinal samples of two chicken genetic lines. First, imputation accuracy is evaluated in 12 samples for which both low- and high-depth sequencing data are available, obtaining high imputation accuracies for all tested panels (>0.90). Second, the impact of reference panel choice in population genetics statistics on 100 chickens is assessed, all four panels yielding comparable results. In light of the observations, the feasibility and application of the applied imputation strategy are discussed for different species with regard to the host DNA proportion, genomic diversity, and availability of a reference panel. This method enables leveraging insofar discarded host DNA to get insights into the genetic structure of host populations, and in doing so, facilitates the implementation of hologenomic approaches that jointly analyze host and microbial genomic data.

A short exposure to a semi-natural habitat alleviates the honey bee hive microbial imbalance caused by agricultural stress.

Honeybee health and the species' gut microbiota are interconnected. Also noteworthy are the multiple niches present within hives, each with distinct microbiotas and all coexisting, which we termed "apibiome". External stressors (e.g. anthropization) can compromise microbial balance and bee resilience. We hypothesised that (1) the bacterial communities of hives located in areas with different degrees of anthropization differ in composition, and (2) due to interactions between the multiple microbiomes within the apibiome, changes in the community of a niche would impact the bacteria present in other hive sections. We characterised the bacterial consortia of different niches (bee gut, bee bread, hive entrance and internal hive air) of 43 hives from 3 different environments (agricultural, semi-natural and natural) through 16S rRNA amplicon sequencing. Agricultural samples presented lower community evenness, depletion of beneficial bacteria, and increased recruitment of stress related pathways (predicted via PICRUSt2). The taxonomic and functional composition of gut and hive entrance followed an environmental gradient. Arsenophonus emerged as a possible indicator of anthropization, gradually decreasing in abundance from agriculture to the natural environment in multiple niches. Importantly, after 16 days of exposure to a semi-natural landscape hives showed intermediate profiles, suggesting alleviation of microbial dysbiosis through reduction of anthropization.

Population Structure and Diversity in European Honey Bees (Apismellifera L.)-An Empirical Comparison of Pool and Individual Whole-Genome Sequencing.

BACKGROUND: Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera). METHODS: We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches. RESULTS: Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers. CONCLUSIONS: If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.

Novel strategies to improve chicken performance and welfare by unveiling host-microbiota interactions through hologenomics.

Fast optimisation of farming practices is essential to meet environmental sustainability challenges. Hologenomics, the joint study of the genomic features of animals and the microbial communities associated with them, opens new avenues to obtain in-depth knowledge on how host-microbiota interactions affect animal performance and welfare, and in doing so, improve the quality and sustainability of animal production. Here, we introduce the animal trials conducted with broiler chickens in the H2020 project HoloFood, and our strategy to implement hologenomic analyses in light of the initial results, which despite yielding negligible effects of tested feed additives, provide relevant information to understand how host genomic features, microbiota development dynamics and host-microbiota interactions shape animal welfare and performance. We report the most relevant results, propose hypotheses to explain the observed patterns, and outline how these questions will be addressed through the generation and analysis of animal-microbiota multi-omic data during the HoloFood project.

Transcriptomic dataset for Sardina pilchardus: Assembly, annotation, and expression of nine tissues.

European sardine or pilchard is a planktonic small pelagic fish present from the North Sea in Europe to the coast of Senegal in the North of Africa, and across the Mediterranean sea to the Black Sea. Ecologically, sardines are an intermediary link in the trophic network, preying on plankton and being predated by larger fishes, marine mammals, and seabirds. This species is of great nutritional and economic value as a cheap but rich source of protein and fat. It is either consumed directly by humans or fed as fishmeal for aquaculture and farm animals. Despite its importance in the food basket, little is known about the molecular mechanisms involved in protein and lipid synthesis in this species. We collected nine tissues of Sardina pilchardus and reconstructed the transcriptome. In all, 198,597 transcripts were obtained, from which 68,031 are protein-coding. Quality assessment of the transcriptome was performed by back-mapping reads to the transcriptome and by searching for Single Copy Orthologs. Additionally, Gene Ontology and KEGG annotations were retrieved for most of the protein-coding genes. Finally, each library was quantified in terms of Transcripts per Million to disclose their expression patterns.

Impact of dimethylpyrazole-based nitrification inhibitors on soil-borne bacteria.

Nitrogen (N) input from fertilizers modifies the properties of agricultural soils as well as bacterial community diversity, composition and relationships. This can lead to negative impacts such as the deterioration of system multifunctionality, whose maintenance is critical to normal nutrient cycling. Synthetic nitrification inhibitors (NIs) can be combined with fertilizers to improve the efficiency of N use by reducing N losses. However, analysis of their effects on non-target bacteria are scarce. This study aimed to analyze the effect of applying the NIs DMPP and DMPSA on the whole bacterial community. Through 16S rRNA amplicon sequencing we determined the differences between samples in terms of microbial diversity, composition and co-occurrence networks. The application of DMPP and DMPSA exerted little impact on the abundance of the dominant phyla. Nevertheless, several significant shifts were detected in bacterial diversity, co-occurrence networks, and the abundance of particular taxa, where soil water content played a key role. For instance, the application of NIs intensified the negative impact of N fertilization on bacterial diversity under high water-filled pore spaces (WFPS) (>64%), reducing community diversity, whereas alpha-diversity was not affected at low WFPS (<55%). Interestingly, despite NIs are known to inhibit ammonia monooxygenase (AMO) enzyme, both NIs almost exclusively inhibited Nitrosomonas genera among AMO holding nitrifiers. Thus, Nitrosomonas showed abundance reductions of up to 47% (DMPP) and 66% (DMPSA). Nonetheless, non-target bacterial abundances also shifted with NI application. Notably, DMPSA application partially alleviated the negative effect of fertilization on soil multifunctionality. A remarkable increase in populations related to system multifunctionality, such as Armatimonadetes (up to +21%), Cyanobacteria (up to +30%) and Fibrobacteres (up to +25%) was observed when DMPSA was applied. NI application substantially influenced microbial associations by decreasing the complexity of co-occurrence networks, decreasing the total edges and node connectivity, and increasing path distances.

High-density SNP genotyping detects homogeneity of Spanish and French Basques, and confirms their genomic distinctiveness from other European populations.

A recent study reported that Basques do not constitute a genetically distinct population, and that Basques from Spanish and French provinces do not show significant genetic similarity. These conclusions disagree with numerous previous studies, and are not consistent with the historical and linguistic evidence that supports the distinctiveness of Basques. In order to further investigate this controversy, we have genotyped 83 Spanish Basque individuals and used these data to infer population structure based on more than 60,000 single nucleotide polymorphisms of several European populations. Here, we present the first high-throughput analysis including Basques from Spanish and French provinces, and show that all Basques constitute a homogeneous group that can be clearly differentiated from other European populations.

EXFI: Exon and splice graph prediction without a reference genome.

For population genetic studies in nonmodel organisms, it is important to use every single source of genomic information. This paper presents EXFI, a Python pipeline that predicts the splice graph and exon sequences using an assembled transcriptome and raw whole-genome sequencing reads. The main algorithm uses Bloom filters to remove reads that are not part of the transcriptome, to predict the intron-exon boundaries, to then proceed to call exons from the assembly, and to generate the underlying splice graph. The results are returned in GFA1 format, which encodes both the predicted exon sequences and how they are connected to form transcripts. EXFI is written in Python, tested on Linux platforms, and the source code is available under the MIT License at https://github.com/jlanga/exfi.

Digging into the Genomic Past of Swiss Honey Bees by Whole-Genome Sequencing Museum Specimens.

Historical specimens in museum collections provide opportunities to gain insights into the genomic past. For the Western honey bee, Apis mellifera L., this is particularly important because its populations are currently under threat worldwide and have experienced many changes in management and environment over the last century. Using Swiss Apis mellifera mellifera as a case study, our research provides important insights into the genetic diversity of native honey bees prior to the industrial-scale introductions and trade of non-native stocks during the 20th century-the onset of intensive commercial breeding and the decline of wild honey bees following the arrival of Varroa destructor. We sequenced whole-genomes of 22 honey bees from the Natural History Museum in Bern collected in Switzerland, including the oldest A. mellifera sample ever sequenced. We identify both, a historic and a recent migrant, natural or human-mediated, which corroborates with the population history of honey bees in Switzerland. Contrary to what we expected, we find no evidence for a significant genetic bottleneck in Swiss honey bees, and find that genetic diversity is not only maintained, but even slightly increased, most probably due to modern apicultural practices. Finally, we identify signals of selection between historic and modern honey bee populations associated with genes enriched in functions linked to xenobiotics, suggesting a possible selective pressure from the increasing use and diversity of chemicals used in agriculture and apiculture over the last century.

Response of Horticultural Soil Microbiota to Different Fertilization Practices.

Environmentally friendly agricultural production necessitates manipulation of microbe-plant interactions, requiring a better understanding of how farming practices influence soil microbiota. We studied the effect of conventional and organic treatment on soil bacterial richness, composition, and predicted functional potential. 16S rRNA sequencing was applied to soils from adjacent plots receiving either a synthetic or organic fertilizer, where two crops were grown within treatment, homogenizing for differences in soil properties, crop, and climate. Conventional fertilizer was associated with a decrease in soil pH, an accumulation of Ag, Mn, As, Fe, Co, Cd, and Ni; and an enrichment of ammonia oxidizers and xenobiotic compound degraders (e.g., Candidatus Nitrososphaera, Nitrospira, Bacillus, Pseudomonas). Soils receiving organic fertilization were enriched in Ti (crop biostimulant), N, and C cycling bacteria (denitrifiers, e.g., Azoarcus, Anaerolinea; methylotrophs, e.g., Methylocaldum, Methanosarcina), and disease-suppression (e.g., Myxococcales). Some predicted functions, such as glutathione metabolism, were slightly, but significantly enriched after a one-time manure application, suggesting the enhancement of sulfur regulation, nitrogen-fixing, and defense of environmental stressors. The study highlights that even a single application of organic fertilization is enough to originate a rapid shift in soil prokaryotes, responding to the differential substrate availability by promoting soil health, similar to recurrent applications.

A novel transcriptome-derived SNPs array for tench (Tinca tinca L.).

Tench (Tinca tinca L.) has great economic potential due to its high rate of fecundity and long-life span. Population genetic studies based on allozymes, microsatellites, PCR-RFLP and sequence analysis of genes and DNA fragments have revealed the presence of Eastern and Western phylogroups. However, the lack of genomic resources for this species has complicated the development of genetic markers. In this study, the tench transcriptome and genome were sequenced by high-throughput sequencing. A total of 60,414 putative SNPs were identified in the tench transcriptome using a computational pipeline. A set of 96 SNPs was selected for validation and a total of 92 SNPs was validated, resulting in the highest conversion and validation rate for a non-model species obtained to date (95.83%). The validated SNPs were used to genotype 140 individuals belonging to two tench breeds (Tabor and Hungarian), showing low (FST = 0.0450) but significant (<0.0001) genetic differentiation between the two tench breeds. This implies that set of validated SNPs array can be used to distinguish the tench breeds and that it might be useful for studying a range of associations between DNA sequence and traits of importance. These genomic resources created for the tench will provide insight into population genetics, conservation fish stock management, and aquaculture.

Hygroregulation, a key ability for eusocial insects: Native Western European honeybees as a case study.

Sociality has brought many advantages to various hymenoptera species, including their ability of regulating physical factors in their nest (e.g., temperature). Although less studied, humidity is known to be important for egg, larval and pupal development, and also for nectar concentration. Two subspecies of Apis mellifera of the M evolutionary lineage were used as models to test the ability of a superorganism (i.e. honeybee colony) to regulate the humidity in its nest (i.e. "hygroregulation hypothesis") in four conservation centers: two in France (A. m. mellifera) and two in Portugal (A. m. iberiensis). We investigated the ability of both subspecies to regulate the humidity in hives daily, but also during the seasons for one complete year. Our data and statistical analysis demonstrated the capacity of the bees to regulate humidity in their hive, regardless of the day, season or subspecies. Furthermore, the study showed that humidity in beehives is stable even during winter, when brood is absent, and when temperature is known to be less stable in the beehives. These results suggest that humidity is important for honeybees at every life stage, maybe because of the 'imprint' of the evolutionary history of this hymenopteran lineage.

Association between combinations of genetic polymorphisms and epidemiopathogenic forms of bovine paratuberculosis.

Control of major mycobacterial diseases affecting livestock is a challenging issue that requires different approaches. The use of genetic markers for improving resistance to Mycobacterium avium subsp. paratuberculosis infection in cattle has been explored as a promising population strategy We performed paratuberculosis epidemiopathogenic phenotypic and genotypic characterization involving 24 SNPs in six candidate genes (NOD2, CD209, SLC11A1, SP110, TLR2 and TLR4) on 502 slaughtered Friesian cows. In the current study, we investigate whether recently proposed paratuberculosis (PTB) epidemiopathogenic (EP) forms (apparently free-AF, latent-LAT and patent-PAT) could be associated with some combination of these 24 SNPs. Best EP form grouping was obtained using a combination of 5 SNPs in four genes (CD209: rs210748127; SLC11A1: rs110090506; SP110: rs136859213 and rs110480812; and TLR2: rs41830058). These groups were defined according to the level of infection progression risk to patent epidemiopathogenic forms and showed the following distributions: LOWIN (low) with 39 (8%) cases (94.9% AF/5.1% LAT/0% PAT); LATIN (low) with 17 (3%) cases (5.9% AF/94.1% LAT/0% PAT); AVERIN (average) with 413 (82%) cases (52.1% AF/38.5% LAT/9.4% PAT) and PATIN (patent) with 33 (7%) cases (36.4% AF/24.2% LAT/39.4% PAT). Age of slaughter was significantly higher for LATIN (88.3 months) compared to AVERIN (65.3 months; p = 0.0007) and PATIN (59.1 months; p = 0.0004), and for LOWIN (73.9 months) compared to PATIN (p = 0.0233), and nearly significant compared to AVERIN (p = 0.0572) These results suggest that some selected genetic polymorphisms have a potential use as markers of PTB EP forms and thus add a new tool for the control of this widespread infection.