Africanized honey bee colonies in Costa Rica: first evidence of its management, brood nest structure and factors associated with varroa mite infestation.

Management, brood nest structure and factors associated with varroa mite infestation were studied in 60 apiaries of Africanized honey bees in the northwest region of the Central Valley of Costa Rica. Apiaries were monitored two times. The first monitoring was taken forward during the rainy season between May and November 2019. The second monitoring during the dry season between February and March 2020. Information about the beekeepers, apiaries and management was collected through a survey. Amount of open and capped brood, honey and pollen were measured in the field. The infestation rate of varroa (IRV) was quantified using standard laboratory methods. A determination of multi-residue pesticides in bee bread was made through GC-MS/MS and LC-MS/MS techniques. According to the results, most of the beekeepers produce honey (96.7%), participate in training activities (82.2%), and change the bee queens annually (70%). The first monitoring was characterized by a lower amount of capped brood and honey reserves compared to the second one. IRV was significantly higher in the first monitoring (6.0 ± 0.4) in comparison with the second one (3.0 ± 0.3) (U Mann-Whitney p < 0.001). The maximum value for the first monitoring exceeds 40%, while this value was close to 25% in the second monitoring. Mite infestation exposed significant differences in relation to the variables associated to the beekeeper's management, i.e., change of bee queen (p = 0.002) or when beekeepers monitor varroa mites (p = 0.004). Additionally, the IRV had inverse correlations (p < 0.01) with the number of comb sides with capped brood (Spearman's rho coefficient = - 0.190), and honey reserves (Spearman's rho coefficient = - 0.168). Furthermore, 23 of 60 bee bread samples presented one to five pesticide residues, being the most frequent antifungal agrochemicals.

The demographic history and adaptation of Canarian goat breeds to environmental conditions through the use of genome-wide SNP data.

BACKGROUND: The presence of goats in the Canary Islands dates back to the late 1st millennium BC, which coincides with the colonization by the Amazigh settlers. However, the exact geographic origin of Canarian goats is uncertain since the Amazigh peoples were distributed over a wide spatial range. Nowadays, three Canarian breeds (Palmera, Majorera and Tinerfeña) are officially recognized, along with two distinct South and North Tinerfeña ecotypes, with the South Tinerfeña and Majorera goats thriving in arid and dry semi-desertic environments and the Palmera and North Tinerfeña goats are adapted to humid and temperate areas that are influenced by trade winds. Genotypes for 224 Canarian goats were generated using the Illumina Goat single nucleotide polymorphism (SNP)50 BeadChip. By merging these data with the genotypes from 1007 individuals of African and Southern European ancestry, our aim was to ascertain the geographic origin of the Canarian goats and identify genes associated with adaptation to diverse environmental conditions. RESULTS: The diversity indices of the Canarian breeds align with most of those of the analyzed local breeds from Africa and Europe, except for the Palmera goats that showed lower levels of genetic variation. The Canarian breeds demonstrate a significant genetic differentiation compared to other populations, which indicates a history of prolonged geographic isolation. Moreover, the phylogenetic reconstruction indicated that the ancestry of the Canarian goats is fundamentally North African rather than West African. The ADMIXTURE and the TreeMix analyses showed no evidence of gene flow between Canarian goats and other continental breeds. The analysis of runs of homozygosity (ROH) identified 13 ROH islands while the window-based FST method detected 25 genomic regions under selection. Major signals of selection were found on Capra hircus (CHI) chromosomes 6, 7, and 10 using various comparisons and methods. CONCLUSIONS: This genome-wide analysis sheds new light on the evolutionary history of the four breeds that inhabit the Canary Islands. Our findings suggest a North African origin of the Canarian goats. In addition, within the genomic regions highlighted by the ROH and FST approaches, several genes related to body size and heat tolerance were identified.

Single and longitudinal genome-wide association studies for dairy traits available in goats with three recorded lactations.

Milk yield and composition phenotypes are systematically recorded across several lactations in goats, but the majority of genome-wide association studies (GWAS) performed so far have rather ignored the longitudinal nature of such data. Here, we have used two different GWAS approaches to analyse data from three lactations recorded in Murciano-Granadina goats. In Analysis 1, independent GWAS have been carried out for each trait and lactation, while a single longitudinal GWAS, jointly considering all data, has been performed in Analysis 2. In both analyses, genome-wide significant QTL for lactose percentage on chromosome 2 (129.77-131.01 Mb) and for milk protein percentage on the chromosome 6 (74.8-94.6 Mb) casein gene cluster region were detected. In Analysis 1, several QTL were not replicated in all three lactations, possibly due to the existence of lactation-specific genetic determinants. In Analysis 2, we identified several genome-wide significant QTL related to milk yield and protein content that were not uncovered in Analysis 1. The increased number of QTL identified in Analysis 2 suggests that the longitudinal GWAS is particularly well suited for the genetic analysis of dairy traits. Moreover, our data confirm that variability within or close to the casein complex is the main genetic determinant of milk protein percentage in Murciano-Granadina goats.

Univariate and multivariate genome-wide association studies for hematological traits in Murciano-Granadina goats.

Hematological traits are important indicators of health status, and they are frequently used as criteria for clinical diagnosis. In humans, the genomic architecture of blood traits has been investigated in depth and thousands of associations with genetic variants have been found. In contrast, the association between marker genotypes and the variation of hematological traits has not been investigated in goats yet. Herewith, we have recorded 12 hematological parameters in 882 Murciano-Granadina goats that were also genotyped with the Goat SNP50 BeadChip (Illumina). Performance of a univariate genome-wide association study (GWAS) made it possible to detect one genomic region on goat chromosome (CHI) 21 (19.2-19.5 Mb) associated, at the genome-wide level of significance, with 4 red blood cell traits. The three markers displaying the highest significances were rs268272996 (CHI21: 19225290 bp), rs268273004 (CHI21: 19565629 bp) and rs268239059 (CHI13: 9615190 bp). Consistently, a multivariate GWAS indicated that the rs268273004 marker on chromosome 21 is associated with seven blood cell traits. Interestingly, this marker maps close to the FA Complementation Group I (FANCI) gene (CHI21: 20021947-20077025 bp), which is functionally related to Fanconi anemia, a syndrome characterized by bone marrow failure, aplastic anemia, and congenital disorders. We have also uncovered additional chromosome-wide significant associations between genetic markers and erythrocyte and leukocyte traits in the univariate GWAS. These findings evidence that the phenotypic variation of hematological traits in goats is regulated, at least to some extent, by polygenic determinants distributed in multiple chromosomes.

Genomic patterns of homozygosity and inbreeding depression in Murciano-Granadina goats.

BACKGROUND: Inbreeding depression can adversely affect traits related to fitness, reproduction and productive performance. Although current research suggests that inbreeding levels are generally low in most goat breeds, the impact of inbreeding depression on phenotypes of economic interest has only been investigated in a few studies based on genealogical data. RESULTS: We genotyped 1040 goats with the Goat SNP50 BeadChip. This information was used to estimate different molecular inbreeding coefficients and characterise runs of homozygosity and homozygosity patterns. We detected 38 genomic regions with increased homozygosity as well as 8 ROH hotspots mapping to chromosomes 1, 2, 4, 6, 14, 16 and 17. Eight hundred seventeen goats with available records for dairy traits were analysed to evaluate the potential consequences of inbreeding depression on milk phenotypes. Four regions on chromosomes 8 and 25 were significantly associated with inbreeding depression for the natural logarithm of the somatic cell count. Notably, these regions contain several genes related with immunity, such as SYK, IL27, CCL19 and CCL21. Moreover, one region on chromosome 2 was significantly associated with inbreeding depression for milk yield. CONCLUSIONS: Although genomic inbreeding levels are low in Murciano-Granadina goats, significant evidence of inbreeding depression for the logarithm of the somatic cell count, a phenotype closely associated with udder health and milk yield, have been detected in this population. Minimising inbreeding would be expected to augment economic gain by increasing milk yield and reducing the incidence of mastitis, which is one of the main causes of dairy goat culling.

On the origins of American Criollo pigs: A common genetic background with a lasting Iberian signature.

American Criollo pigs are thought to descend mainly from those imported from the Iberian Peninsula starting in the late 15th century. Criollo pigs subsequently expanded throughout the Americas, adapting to very diverse environments, and possibly receiving influences from other origins. With the intensification of agriculture in the mid-20th century, cosmopolitan breeds largely replaced Criollo pigs, and the few remaining are mostly maintained by rural communities in marginal areas where they still play an important socio-economic and cultural role. In this study, we used 24 microsatellite markers in samples from 1715 pigs representing 46 breeds with worldwide distribution, including 17 American Criollo breeds, with the major focus of investigating their genetic diversity, structure and breed relationships. We also included representatives of the Iberian, Local British, Hungarian, Chinese and Commercial breeds, as well as Wild Boar, in order to investigate their possible influence in the genetic composition of Criollos. Our results show that, when compared with the other breeds, Criollo pigs present higher levels of genetic diversity, both in terms of allelic diversity and expected heterozygosity. The various analyses indicate that breed differentiation overall explains nearly 21% of the total genetic diversity. Criollo breeds showed their own identity and shared a common genetic background, tending to cluster together in various analyses, even though they differ from each other. A close relationship of Criollos with Iberian breeds was revealed by all the different analyses, and the contribution of Iberian breeds, particularly of the Celtic breeds, is still present in various Criollo breeds. No influence of Chinese breeds was detected on Criollos, but a few were influenced by Commercial breeds or by wild pigs. Our results confirm the uniqueness of American Criollo pigs and the role that Iberian breeds have played in their development.

Author Correction: The genetic ancestry of American Creole cattle inferred from uniparental and autosomal genetic markers.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A genome-wide association analysis for body, udder, and leg conformation traits recorded in Murciano-Granadina goats.

Morphological traits are of great importance to dairy goat production given their effect on phenotypes of economic interest. However, their underlying genomic architecture has not yet been extensively characterized. Herein, we aimed to identify genomic regions associated with body, udder, and leg conformation traits recorded in 825 Murciano-Granadina goats. We genotyped this resource population using the GoatSNP50 BeadChip (Illumina Inc., San Diego, CA) and performed genome-wide association analyses using the GEMMA software. We found 2 genome-wide significant associations between markers rs268273468 [Capra hircus (CHI) 16:69617700] and rs268249346 (CHI 28:18321523) and medial suspensory ligament. In contrast, we did not detect any genome-wide significant associations for body and leg traits. Moreover, we found 12, 19, and 7 chromosome-wide significant associations for udder, body, and leg traits, respectively. Comparison of our data with previous studies revealed a low level of positional concordance between regions associated with morphological traits. In addition to technical factors, this lack of concordance could be due to a substantial level of genetic heterogeneity among breeds or to the strong polygenic background of morphological traits, which makes it difficult to detect genetic factors that have small phenotypic effects.

Diversity Analysis and Genetic Relationships among Local Brazilian Goat Breeds Using SSR Markers.

The genetic diversity of six Brazilian native goats was reported using molecular markers. Hair samples of 332 animals were collected from different goat breeds (Moxotó, Canindé, Serrana Azul, Marota, Repartida, and Graúna) from five states of Northeast Brazil (Paraíba, Pernambuco, Rio Grande do Norte, Bahia, and Piauí). A panel of 27 microsatellites or single sequence repeats (SSRs) were selected and amplified using a polymerase chain reaction (PCR) technique. All populations showed an average allele number of over six. The mean observed heterozygosity for Brazilian breeds was superior to 0.50. These results demonstrated the high genetic diversity in the studied populations with values ranging from 0.53 (Serrana Azul) to 0.62 (Repartida). The expected average heterozygosity followed the same trend ranging from 0.58 (Serrana Azul) to 0.65 (Repartida), and the values obtained are very similar for all six breeds. The fixation index (Fis) had values under 10% except for the Moxotó breed (13%). The mean expected heterozygosity of all Brazilian populations was over 0.50. Results indicated a within-breed genetic variability in the Brazilian breeds based on the average number of alleles and the average observed heterozygosity. The interbreed genetic diversity values showed proper genetic differentiation among local Brazilian goat breeds.

A genome-wide analysis of copy number variation in Murciano-Granadina goats.

BACKGROUND: In this work, our aim was to generate a map of the copy number variations (CNV) segregating in a population of Murciano-Granadina goats, the most important dairy breed in Spain, and to ascertain the main biological functions of the genes that map to copy number variable regions. RESULTS: Using a dataset that comprised 1036 Murciano-Granadina goats genotyped with the Goat SNP50 BeadChip, we were able to detect 4617 and 7750 autosomal CNV with the PennCNV and QuantiSNP software, respectively. By applying the EnsembleCNV algorithm, these CNV were assembled into 1461 CNV regions (CNVR), of which 486 (33.3% of the total CNVR count) were consistently called by PennCNV and QuantiSNP and used in subsequent analyses. In this set of 486 CNVR, we identified 78 gain, 353 loss and 55 gain/loss events. The total length of all the CNVR (95.69 Mb) represented 3.9% of the goat autosomal genome (2466.19 Mb), whereas their size ranged from 2.0 kb to 11.1 Mb, with an average size of 196.89 kb. Functional annotation of the genes that overlapped with the CNVR revealed an enrichment of pathways related with olfactory transduction (fold-enrichment = 2.33, q-value = 1.61 × 10-10), ABC transporters (fold-enrichment = 5.27, q-value = 4.27 × 10-04) and bile secretion (fold-enrichment = 3.90, q-value = 5.70 × 10-03). CONCLUSIONS: A previous study reported that the average number of CNVR per goat breed was ~ 20 (978 CNVR/50 breeds), which is much smaller than the number we found here (486 CNVR). We attribute this difference to the fact that the previous study included multiple caprine breeds that were represented by small to moderate numbers of individuals. Given the low frequencies of CNV (in our study, the average frequency of CNV is 1.44%), such a design would probably underestimate the levels of the diversity of CNV at the within-breed level. We also observed that functions related with sensory perception, metabolism and embryo development are overrepresented in the set of genes that overlapped with CNV, and that these loci often belong to large multigene families with tens, hundreds or thousands of paralogous members, a feature that could favor the occurrence of duplications or deletions by non-allelic homologous recombination.

Analyzing the genomic and transcriptomic architecture of milk traits in Murciano-Granadina goats.

BACKGROUND: In this study, we aimed to investigate the molecular basis of lactation as well as to identify the genetic factors that influence milk yield and composition in goats. To achieve these two goals, we have analyzed how the mRNA profile of the mammary gland changes in seven Murciano-Granadina goats at each of three different time points, i.e. 78 d (T1, early lactation), 216 d (T2, late lactation) and 285 d (T3, dry period) after parturition. Moreover, we have performed a genome-wide association study (GWAS) for seven dairy traits recorded in the 1st lactation of 822 Murciano-Granadina goats. RESULTS: The expression profiles of the mammary gland in the early (T1) and late (T2) lactation were quite similar (42 differentially expressed genes), while strong transcriptomic differences (more than one thousand differentially expressed genes) were observed between the lactating (T1/T2) and non-lactating (T3) mammary glands. A large number of differentially expressed genes were involved in pathways related with the biosynthesis of amino acids, cholesterol, triglycerides and steroids as well as with glycerophospholipid metabolism, adipocytokine signaling, lipid binding, regulation of ion transmembrane transport, calcium ion binding, metalloendopeptidase activity and complement and coagulation cascades. With regard to the second goal of the study, the performance of the GWAS allowed us to detect 24 quantitative trait loci (QTLs), including three genome-wide significant associations: QTL1 (chromosome 2, 130.72-131.01 Mb) for lactose percentage, QTL6 (chromosome 6, 78.90-93.48 Mb) for protein percentage and QTL17 (chromosome 17, 11.20 Mb) for both protein and dry matter percentages. Interestingly, QTL6 shows positional coincidence with the casein genes, which encode 80% of milk proteins. CONCLUSIONS: The abrogation of lactation involves dramatic changes in the expression of genes participating in a broad array of physiological processes such as protein, lipid and carbohydrate metabolism, calcium homeostasis, cell death and tissue remodeling, as well as immunity. We also conclude that genetic variation at the casein genes has a major impact on the milk protein content of Murciano-Granadina goats.

Effect of olive-derived antioxidants (3,4-dihydroxyphenylethanol and 3,4 dihydroxyphenylglycol) on sperm motility and fertility in liquid ram sperm stored at 15°C or 5°C.

The aim of the present study was to evaluate the effect of the addition of two olive oil-derived antioxidants, hydroxytyrosol (3,4-dihydroxyphenylethanol, HT) and 3,4-dihydroxyphenylglycol (DHPG), on ovine semen during liquid storage at 5°C and 15°C. Semen was collected, pooled, diluted and then divided into aliquots supplemented with different concentrations (5 µg/ml, 10 µg/ml, 50 µg/ml and 100 µg/ml) of HT, DHPG and a mixture (MIX) of both antioxidants. Sperm motility characteristics were assessed in the different samples at 0, 6, 24, 48, 72 and 96 hr after cooling, and a fertility trial was also conducted. The results showed that the antioxidant addition did not significantly improve total and progressive motility in ovine cooled sperm maintained at 15° or 5°C. However, in samples stored at 5°C, LIN (48, 72, 96 hr), STR (0 hr) and WOB (0, 48, 72, 96 hr) values significantly decreased in comparison with control treatment when high antioxidant concentrations were added (MIX100 or HT100). When samples were maintained at 15°C, MIX50 showed significantly higher VCL values than the control treatment after 6 hr cooling, and MIX100 showed significantly lower VCL values at 96 hr after cooling. According to the artificial insemination trial, no significant differences were observed when antioxidants were added. In conclusion, the use of HT and DHPG showed small impact in sperm motility and fertility was not affected (nor detrimentally nor positively) when insemination was carried out using antioxidant-supplemented liquid sperm.

Effect of three commercial extenders on sperm motility and fertility in liquid ram semen stored at 15 °C or 5 °C.

The effect of different extenders on sperm motility and fertility was evaluated during liquid storage of ram semen at 5 °C and 15 °C. The semen was collected, pooled and diluted in three commercial extenders: Inra 96® (INRA) based on skimmed milk, Biladyl® A fraction (BIL) based on egg yolk, and Ovixcell® (OVIX) based on soybean lecithin. Then, sperm motility was evaluated at 0, 6, 24, 48, 72 and 96 h. In order to evaluate fertility, samples stored at 15 °C were used after dilution in INRA and OVIX. Results showed that progressive motility was significantly higher up to 72 h of storage in sperm samples maintained at 5 °C in comparison with 15 °C, similarly for each tested diluent. When samples were stored at 5 °C in OVIX, kinematic parameters such as velocity (except curvilinear velocity, VCL), trajectory [linearity (LIN), straightness (STR), wobble (WOB)], amplitude of lateral head displacement (ALH) and beat/cross frequency (BCF) were higher than in INRA and BIL. No significant differences in pregnancy rate were detected between INRA (62.6%) and OVIX (58.9%). In conclusion, liquid storage at 5 °C with OVIX extender is an interesting option since non-animal components are used, and this extender offers similar in vitro and in vivo efficacy as other extenders containing animal components.

Deciphering the Patterns of Genetic Admixture and Diversity in the Ecuadorian Creole Chicken.

Latin American Creole chickens are generally not characterized; this is the case in Ecuador, where the lack of scientific information is contributing to their extinction. Here, we developed a characterization of the genetic resources of Ecuadorian chickens located in three continental agroecosystems (Pacific coastal, Andean, and Amazonian). Blood samples of 234 unrelated animals were collected in six provinces across Ecuador: Bolívar, Chimborazo, Cotopaxi, Guayas, Morona Santiago, and Tungurahua, in order to perform a genetic characterization and population structure assessment using the AVIANDIV project microsatellites panel (30 loci) and D-loop sequences of mitochondrial DNA and comparing with reference data from other breeds or genetic lines. The results indicate that Ecuadorian Creole chickens are the result of the admixture of different genetic groups that occurred during the last five centuries. While the influence of South Spanish breeds is demonstrated in the colonial age, genetic relationships with other breeds (Leghorn, Spanish fighter cock) cannot be discarded. The geographical configuration of the country and extreme climate variability have influenced the genetic isolation of groups constituting a homogeneous genetic status into the whole population. This is not only a source of genetic variation, but also a critical point because genetic drift produces a loss of genetic variants.

The genetic ancestry of American Creole cattle inferred from uniparental and autosomal genetic markers.

Cattle imported from the Iberian Peninsula spread throughout America in the early years of discovery and colonization to originate Creole breeds, which adapted to a wide diversity of environments and later received influences from other origins, including zebu cattle in more recent years. We analyzed uniparental genetic markers and autosomal microsatellites in DNA samples from 114 cattle breeds distributed worldwide, including 40 Creole breeds representing the whole American continent, and samples from the Iberian Peninsula, British islands, Continental Europe, Africa and American zebu. We show that Creole breeds differ considerably from each other, and most have their own identity or group with others from neighboring regions. Results with mtDNA indicate that T1c-lineages are rare in Iberia but common in Africa and are well represented in Creoles from Brazil and Colombia, lending support to a direct African influence on Creoles. This is reinforced by the sharing of a unique Y-haplotype between cattle from Mozambique and Creoles from Argentina. Autosomal microsatellites indicate that Creoles occupy an intermediate position between African and European breeds, and some Creoles show a clear Iberian signature. Our results confirm the mixed ancestry of American Creole cattle and the role that African cattle have played in their development.

A genomic map of climate adaptation in Mediterranean cattle breeds.

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short-term evolutionary response to climate pressure induced by their post-domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome-wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds.

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.

Measuring and modeling for the assessment of the genetic background behind cognitive processes in donkeys.

New productive niches can offer new commercial perspectives linked to donkeys' products and human therapeutic or leisure applications. However, no assessment for selection criteria has been carried out yet. First, we assessed the animal inherent features and environmental factors that may potentially influence several cognitive processes in donkeys. Then, we aimed at describing a practical methodology to quantify such cognitive processes, seeking their inclusion in breeding and conservation programmes, through a multifactorial linear model. Sixteen cognitive process-related traits were scored on a problem-solving test in a sample of 300 Andalusian donkeys for three consecutive years from 2013 to 2015. The linear model assessed the influence and interactions of four environmental factors, sex as an animal-inherent factor, age as a covariable, and the interactions between these factors. Analyses of variance were performed with GLM procedure of SPSS Statistics for Windows, Version 24.0 software to assess the relative importance of each factor. All traits were significantly (P<0.05) affected by all factors in the model except for sex that was not significant for some of the cognitive processes, and stimulus which was not significant (P<0.05) for all of them except for the coping style related ones. The interaction between all factors within the model was non-significant (P<0.05) for almost all cognitive processes. The development of complex multifactorial models to study cognitive processes may counteract the inherent variability in behavior genetics and the estimation and prediction of related breeding parameters, key for the implementation of successful conservation programmes in apparently functionally misplaced endangered breeds.

Detecting the existence of gene flow between Spanish and North African goats through a coalescent approach.

Human-driven migrations are one of the main processes shaping the genetic diversity and population structure of domestic species. However, their magnitude and direction have been rarely analysed in a statistical framework. We aimed to estimate the impact of migration on the population structure of Spanish and African goats. To achieve this goal, we analysed a dataset of 1,472 individuals typed with 23 microsatellites. Population structure of African and Spanish goats was moderate (mean FST = 0.07), with the exception of the Canarian and South African breeds that displayed a significant differentiation when compared to goats from North Africa and Nigeria. Measurement of gene flow with Migrate-n and IMa coalescent genealogy samplers supported the existence of a bidirectional gene flow between African and Spanish goats. Moreover, IMa estimates of the effective number of migrants were remarkably lower than those calculated with Migrate-n and classical approaches. Such discrepancies suggest that recent divergence, rather than extensive gene flow, is the main cause of the weak population structure observed in caprine breeds.

Molecular Study of the Amazonian Macabea Cattle History.

Macabea cattle are the only Bos taurus breed that have adapted to the wet tropical conditions of the Amazon. This breed has integrated into the culture of the indigenous Shuar-Asuar nations probably since its origins, being one of the few European zoogenetic resources assimilated by the deep-jungle Amazon communities. Despite its potential for local endogenous sustainable development, this breed is currently endangered. The present study used molecular genetics tools to investigate the within- and between-breeds diversity, in order to characterize the breed population, define its associations with other breeds, and infer its origin and evolution. The within-breed genetic diversity showed high values, as indicated by all genetic parameters, such as the mean number of alleles (MNA = 7.25±2.03), the observed heterozygosity (Ho = 0.72±0.02) and the expected heterozygosity (He = 0.72±0.02). The between-breeds diversity analysis, which included factorial correspondence analysis, Reynolds genetic distance, neighbor-joining analysis, and genetic structure analysis, showed that the Macabea breed belongs to the group of the American Creoles, with a Southern-Spain origin. Our outcomes demonstrated that the Macabea breed has a high level of purity and null influences of exotic cosmopolitan breeds with European or Asiatic origin. This breed is an important zoogenetic resource of Ecuador, with relevant and unique attributes; therefore, there is an urgent need to develop conservation strategies for the Macabea breed.