Wild Boar (Sus scrofa)-Fascioloides magna Interaction from the Perspective of the MHC Genes.

Fascioloidosis is a parasitic disease caused by a trematode Fascioloides magna. Since major histocompatibility complex (MHC) genes play an important role in the immune response, the aim of this study was to compare the potential differences in MHC class II SLA-DRB1 exon 2 genes between wild boar populations from infected (cases) and non-infected areas (controls). During the winter of 2021, a total of 136 wild boar tissue samples were collected, 39 cases and 97 controls. DNA was extracted and sequenced using the Illumina platform. Differences in distributions of allele combinations were calculated using the Chi-Square test for homogeneity and between proportions using the large-sample test and Fisher-Irwin test. Analysis revealed 19 previously described swine leucocyte antigen (SLA) alleles. The number of polymorphic sites was 79 (29.6%), with 99 mutations in total. Nucleotide diversity π was estimated at 0.11. Proportions of the alleles SLA-DRB1*12:05 (p = 0.0008379) and SLA-DRB1*0101 (p = 0.0002825) were statistically significantly higher in controls, and proportions of the SLA-DRB1*0602 (p = 0.006059) and SLA-DRB1*0901 (p = 0.0006601) in cases. Alleles SLA-DRB1*04:09, SLA-DRB1*0501, SLA-DRB1*11:09, and SLA-DRB1*1301 were detected only in cases, while SLA-DRB1*0404, SLA-DRB1*0701, SLA-DRB1*02:10, and SLA-DRB1*04:08 were present only in controls. We did not confirm the existence of specific alleles that could be linked to F. magna infection. Detected high variability of the MHC class II SLA-DRB1 exon 2 genes indicate high resistance potential against various pathogens.

Performance Comparison of Different Approaches in Genotyping MHC-DRB: The Contrast between Single-Locus and Multi-Locus Species.

Major histocompatibility complex (MHC) genes are widely recognised as valuable markers for wildlife genetic studies given their extreme polymorphism and functional importance in fitness-related traits. Newly developed genotyping methods, which rely on the use of next-generation sequencing (NGS), are gradually replacing traditional cloning and Sanger sequencing methods in MHC genotyping studies. Allele calling in NGS methods remains challenging due to extreme polymorphism and locus multiplication in the MHC coupled with allele amplification bias and the generation of artificial sequences. In this study, we compared the performance of molecular cloning with Illumina and Ion Torrent NGS sequencing in MHC-DRB genotyping of single-locus species (roe deer) and species with multiple DRB loci (red deer) in an attempt to adopt a reliable and straightforward method that does not require complex bioinformatic analyses. Our results show that all methods work similarly well in roe deer, but we demonstrate non-consistency in results across methods in red deer. With Illumina sequencing, we detected a maximum number of alleles in 10 red deer individuals (42), while other methods were somewhat less accurate as they scored 69-81% of alleles detected with Illumina sequencing.

Genetic sub-structuring of Croatian island populations in the Southeastern European context: a meta-analysis.

AIM: To use the method of meta-analysis to assess the influence of island population isolation on the sub-structuring of the Croatian population, as well as the influence of regional population groups on the sub-structuring of the Southeastern European population with regard to basic population genetic statistical parameters calculated by using STR locus analysis. METHODS: Bio-statistical analyses were performed for 2877 unrelated participants of both sexes from Southeastern Europe. Nine autosomal STR loci (D3S1358, vWA, FGA, TH01, TPOX, CSF1PO, D5S818, D13S317, and D7S82) were analyzed by using standard F-statistics and population structure analysis (Structure software). RESULTS: Genetic differentiation of Croatian subpopulations assessed with the FST method was higher at the level of the Croatian population (0.005) than at the level of Southeastern Europe (0.002). The island of Vis showed the most pronounced separation in the Croatian population, and Albanians from Kosovo in the population of Southeast Europe, followed by Croatia, Bosnia and Herzegovina, and Hungary. CONCLUSION: The higher structure of Croatian subpopulations in relation to Southeastern Europe suggest a certain degree of genetic isolation, most likely due to the influence of endogamy within rural island populations.

Population structure and genetic diversity of non-native aoudad populations.

The aoudad (Ammotragus lervia Pallas 1777) is an ungulate species, native to the mountain ranges of North Africa. In the second half of the twentieth century, it was successfully introduced in some European countries, mainly for hunting purposes, i.e. in Croatia, the Czech Republic, Italy, and Spain. We used neutral genetic markers, the mitochondrial DNA control region sequence and microsatellite loci, to characterize and compare genetic diversity and spatial pattern of genetic structure on different timeframes among all European aoudad populations. Four distinct control region haplotypes found in European aoudad populations indicate that the aoudad has been introduced in Europe from multiple genetic sources, with the population in the Sierra Espuña as the only population in which more than one haplotype was detected. The number of detected microsatellite alleles within all populations (< 3.61) and mean proportion of shared alleles within all analysed populations (< 0.55) indicates relatively low genetic variability, as expected for new populations funded by a small number of individuals. In STRUCTURE results with K = 2-4, Croatian and Czech populations cluster in the same genetic cluster, indicating joined origin. Among three populations from Spain, Almeria population shows as genetically distinct from others in results, while other Spanish populations diverge at K = 4. Maintenance of genetic diversity should be included in the management of populations to sustain their viability, specially for small Czech population with high proportion of shared alleles (0.85) and Croatian population that had the smallest estimated effective population size (Ne = 5.4).

Combining phylogenetic and demographic inferences to assess the origin of the genetic diversity in an isolated wolf population.

The survival of isolated small populations is threatened by both demographic and genetic factors. Large carnivores declined for centuries in most of Europe due to habitat changes, overhunting of their natural prey and direct persecution. However, the current rewilding trends are driving many carnivore populations to expand again, possibly reverting the erosion of their genetic diversity. In this study we reassessed the extent and origin of the genetic variation of the Italian wolf population, which is expanding after centuries of decline and isolation. We genotyped wolves from Italy and other nine populations at four mtDNA regions (control-region, ATP6, COIII and ND4) and 39 autosomal microsatellites. Results of phylogenetic analyses and assignment procedures confirmed in the Italian wolves a second private mtDNA haplotype, which belongs to a haplogroup distributed mostly in southern Europe. Coalescent analyses showed that the unique mtDNA haplotypes in the Italian wolves likely originated during the late Pleistocene. ABC simulations concordantly showed that the extant wolf populations in Italy and in south-western Europe started to be isolated and declined right after the last glacial maximum. Thus, the standing genetic variation in the Italian wolves principally results from the historical isolation south of the Alps.

Genetic structure in village dogs reveals a Central Asian domestication origin.

Dogs were the first domesticated species, originating at least 15,000 y ago from Eurasian gray wolves. Dogs today consist primarily of two specialized groups--a diverse set of nearly 400 pure breeds and a far more populous group of free-ranging animals adapted to a human commensal lifestyle (village dogs). Village dogs are more genetically diverse and geographically widespread than purebred dogs making them vital for unraveling dog population history. Using a semicustom 185,805-marker genotyping array, we conducted a large-scale survey of autosomal, mitochondrial, and Y chromosome diversity in 4,676 purebred dogs from 161 breeds and 549 village dogs from 38 countries. Geographic structure shows both isolation and gene flow have shaped genetic diversity in village dog populations. Some populations (notably those in the Neotropics and the South Pacific) are almost completely derived from European stock, whereas others are clearly admixed between indigenous and European dogs. Importantly, many populations--including those of Vietnam, India, and Egypt-show minimal evidence of European admixture. These populations exhibit a clear gradient of short--range linkage disequilibrium consistent with a Central Asian domestication origin.

Genetic structure and admixture between Bayash Roma from northwestern Croatia and general Croatian population: evidence from Bayesian clustering analysis.

The European Roma represent a transnational mosaic of minority population groups with different migration histories and contrasting experiences in their interactions with majority populations across the European continent. Although historical genetic contributions of European lineages to the Roma pool were investigated before, the extent of contemporary genetic admixture between Bayash Roma and non-Romani majority population remains elusive. The aim of this study was to assess the genetic structure of the Bayash Roma population from northwestern Croatia and the general Croatian population and to investigate the extent of admixture between them. A set of genetic data from two original studies (100 Bayash Roma from northwestern Croatia and 195 individuals from the general Croatian population) was analyzed by Bayesian clustering implemented in STRUCTURE software. By re-analyzing published data we intended to focus for the first time on genetic differentiation and structure and in doing so we clearly pointed to the importance of considering social phenomena in understanding genetic structuring. Our results demonstrated that two population clusters best explain the genetic structure, which is consistent with social exclusion of Roma and the demographic history of Bayash Roma who have settled in NW Croatia only about 150 years ago and mostly applied rules of endogamy. The presence of admixture was revealed, while the percentage of non-Croatian individuals in general Croatian population was approximately twofold higher than the percentage of non-Romani individuals in Roma population corroborating the presence of ethnomimicry in Roma.

First evidence of hybridization between golden jackal (Canis aureus) and domestic dog (Canis familiaris) as revealed by genetic markers.

Interspecific hybridization is relatively frequent in nature and numerous cases of hybridization between wild canids and domestic dogs have been recorded. However, hybrids between golden jackals (Canis aureus) and other canids have not been described before. In this study, we combined the use of biparental (15 autosomal microsatellites and three major histocompatibility complex (MHC) loci) and uniparental (mtDNA control region and a Y-linked Zfy intron) genetic markers to assess the admixed origin of three wild-living canids showing anomalous phenotypic traits. Results indicated that these canids were hybrids between golden jackals and domestic dogs. One of them was a backcross to jackal and another one was a backcross to dog, confirming that golden jackal-domestic dog hybrids are fertile. The uniparental markers showed that the direction of hybridization, namely females of the wild species hybridizing with male domestic dogs, was common to most cases of canid hybridization. A melanistic 3bp-deletion at the K locus (ß-defensin CDB103 gene), that was absent in reference golden jackal samples, but was found in a backcross to jackal with anomalous black coat, suggested its introgression from dogs via hybridization. Moreover, we demonstrated that MHC sequences, although rarely used as markers of hybridization, can be also suitable for the identification of hybrids, as long as haplotypes are exclusive for the parental species.

Identification of Coxiella burnetii genotypes in Croatia using multi-locus VNTR analysis.

Although Q fever affects humans and animals in Croatia, we are unaware of genotyping studies of Croatian strains of the causative pathogen Coxiella burnetii, which would greatly assist monitoring and control efforts. Here 3261 human and animal samples were screened for C. burnetii DNA by conventional PCR, and 335 (10.3%) were positive. Of these positive samples, 82 were genotyped at 17 loci using the relatively new method of multi-locus variable number tandem repeat analysis (MLVA). We identified 13 C. burnetii genotypes not previously reported anywhere in the world. Two of these 13 genotypes are typical of the continental part of Croatia and share more similarity with genotypes outside Croatia than with genotypes within the country. The remaining 11 novel genotypes are typical of the coastal part of Croatia and show more similarity to one another than to genotypes outside the country. Our findings shed new light on the phylogeny of C. burnetii strains and may help establish MLVA as a standard technique for Coxiella genotyping.

Multilocus detection of wolf x dog hybridization in italy, and guidelines for marker selection.

Hybridization and introgression can impact the evolution of natural populations. Several wild canid species hybridize in nature, sometimes originating new taxa. However, hybridization with free-ranging dogs is threatening the genetic integrity of grey wolf populations (Canis lupus), or even the survival of endangered species (e.g., the Ethiopian wolf C. simensis). Efficient molecular tools to assess hybridization rates are essential in wolf conservation strategies. We evaluated the power of biparental and uniparental markers (39 autosomal and 4 Y-linked microsatellites, a melanistic deletion at the ß-defensin CBD103 gene, the hypervariable domain of the mtDNA control-region) to identify the multilocus admixture patterns in wolf x dog hybrids. We used empirical data from 2 hybrid groups with different histories: 30 presumptive natural hybrids from Italy and 73 Czechoslovakian wolfdogs of known hybrid origin, as well as simulated data. We assessed the efficiency of various marker combinations and reference samples in admixture analyses using 69 dogs of different breeds and 99 wolves from Italy, Balkans and Carpathian Mountains. Results confirmed the occurrence of hybrids in Italy, some of them showing anomalous phenotypic traits and exogenous mtDNA or Y-chromosome introgression. Hybridization was mostly attributable to village dogs and not strictly patrilineal. The melanistic ß-defensin deletion was found only in Italian dogs and in putative hybrids. The 24 most divergent microsatellites (largest wolf-dog FST values) were equally or more informative than the entire panel of 39 loci. A smaller panel of 12 microsatellites increased risks to identify false admixed individuals. The frequency of F1 and F2 was lower than backcrosses or introgressed individuals, suggesting hybridization already occurred some generations in the past, during early phases of wolf expansion from their historical core areas. Empirical and simulated data indicated the identification of the past generation backcrosses is always uncertain, and a larger number of ancestry-informative markers is needed.

Extensive polymorphism of the major histocompatibility complex DRA gene in Balkan donkeys: perspectives on selection and genealogy.

The major histocompatibility complex (MHC) contains genes important for immune response in mammals, and these genes exhibit high polymorphism and diversity. The DRA gene, a member of the MHC class II family, is highly conserved across a large number of mammalian species, but it displays exceptionally rich sequence variations in Equidae members. We analyzed allelic polymorphism of the DRA locus in 248 donkeys sampled across the Balkan Peninsula (Albania, Bulgaria, Croatia, Macedonia, Greece and Montenegro). Five known alleles and two new alleles were identified. The new allele Eqas-DRA*0601 was found to carry a synonymous mutation, and new allele Eqas-DRA*0701, a non-synonymous mutation. We further analyzed the historical selection and allele genealogy at the DRA locus in equids. Signals of positive selection obtained by various tests were ambiguous. A conservative conclusion is that DRA polymorphism occurred relatively recently and that positive selection has been acting on the DRA locus for a relatively brief period.

Repetitive sequences in Eurasian lynx (Lynx lynx L.) mitochondrial DNA control region.

Mitochondrial DNA (mtDNA) control region (CR) of numerous species is known to include up to five different repetitive sequences (RS1-RS5) that are found at various locations, involving motifs of different length and extensive length heteroplasmy. Two repetitive sequences (RS2 and RS3) on opposite sides of mtDNA central conserved region have been described in domestic cat (Felis catus) and some other felid species. However, the presence of repetitive sequence RS3 has not been detected in Eurasian lynx (Lynx lynx) yet. We analyzed mtDNA CR of 35 Eurasian lynx (L. lynx L.) samples to characterize repetitive sequences and to compare them with those found in other felid species. We confirmed the presence of 80 base pairs (bp) repetitive sequence (RS2) at the 5' end of the Eurasian lynx mtDNA CR L strand and for the first time we described RS3 repetitive sequence at its 3' end, consisting of an array of tandem repeats five to ten bp long. We found that felid species share similar RS3 repetitive pattern and fundamental repeat motif TACAC.

Distribution of microbiological indicators of fecal pollution in the riverine substrates.

Distribution of fecal microorganisms in water, periphyton, and sediment was studied along the Zrnovnica river (Croatia) over a 1.5-year period. It was found that periphyton was inhabited by the highest number of investigated bacteria, while lower numbers of them were found in sediment and the lowest in surface water of the river. The concentrations of fecal microorganisms in periphyton and partly in sediment were found to be significantly higher in the middle of the river course, near the town of Zrnovnica, while according to the analysis of surface water the highest degree of pollution was reached on its estuary. The results were explained with respect to bacterial-algal associations. Considering the fact that most of the river microorganisms are associated with periphyton and sediment particles and only a small number of them is in the free-living form, microbiological analysis of both periphyton and sediment together with water samples has been suggested when fecal pollution of a river is concerned.

Isolation of Clostridium tertium from a striped dolphin (Stenella coeruleoalba) in the Adriatic Sea.

Bacteriologic examination of an abscess found between blubber and musculature of a striped dolphin (Stenella coeruleoalba) found dead on 30 March 2002 in the Adriatic Sea, Croatia, was performed and an aerotolerant, nontoxogenic Clostridium tertium was isolated in pure culture. National Collections of Industrial Food and Marine Bacteria (NCIMB Ltd., Aberdeen, Scotland, UK) confirmed the results. Sequencing results showed it to be C. tertium with 100% similarity. The strain was named Clostridium tertium strain Zagreb, acceded to the culture collection and assigned the accession number NCIMB 13970. This is the first record of C. tertium in marine mammals.