Exploring the miRNAs Profile in Dark-Cutting Beef.

In an animal production system, different stressors may cause the depletion of muscle glycogen stores, resulting in an elevated pH at 24 h post mortem (pH24), which leads to cell metabolism alterations that affect the conversion of muscle into meat, causing meat quality defects, such as dark-cutting beef, also known as dark, firm, and dry (DFD) beef. This process may involve the alteration of small non-coding RNAs (miRNAs), which play critical regulatory roles in cellular processes. Here, we determined whether differential miRNA expression in the Longissimus thoracis et lumborum muscle from the Asturiana de los Valles breed at 24 h post mortem could serve as an early indicator of beef quality defects. Following total RNA extraction, complete miRNAome sequencing revealed 12 miRNAs that were significantly upregulated (p < 0.001) in DFD beef compared to the levels in CONTROL beef. These miRNAs are mainly involved in the cellular responses to redox imbalances and apoptosis. Among these, four miRNAs known to be related to oxidative stress (bta-miR-1246, bta-miR-2332, bta-miR-23b-5p, and bta-miR-2411-3p) were validated via quantitative RT-PCR. Some of their target proteins were also analyzed using Western blotting. High 70 kDa heat shock protein and low Caspase-9 expressions (p < 0.01) were found in DFD beef, suggesting the downregulation of apoptosis. These results suggest the importance of miRNAs in regulating stress in muscle cells during early post mortem, as differences in the abundance of some of these miRNAs are still observed at 24 h post mortem. These changes lead to an inadequate conversion of muscle into meat, resulting in meats with quality defects.

Industrial impact on sustainable dairy farms: Essential elements, hazardous metals and polycyclic aromatic hydrocarbons in forage and cow's milk.

Sustainable dairy farms are characterised by the self-production of forage for animal feed. These farms are sometimes located near industrial areas, entailing a risk of food chain contamination with hazardous metals and polycyclic aromatic hydrocarbons (PAHs). Accordingly, evaluating the impact of pollution on forage and milk is of great interest. In this study, the effects of industrial factors on sustainable forage from 43 dairy farms and possible correlations between inorganic elements and PAHs were studied. Spearman's correlation and principal component analysis (PCA) were performed for the forage and milk. Most of the inorganic elements in the forage were below the maximum residual limits for cadmium (Cd) and lead (Pb), established in EU 2013/1275 and EU 2019/1869, respectively. However, arsenic (As) and mercury (Hg) levels were above their respective limits in the forage (EU 2019/1869). No milk samples exceeded the maximum residual limits for Pb (EU 488/2014) or Cd (EU 1881/2006) in dairy products. Heavy-weight PAHs (HW-PAHs, four or more aromatic rings) were detected in forage but not in milk. In the forage samples, HW-PAHs were positively correlated with Zn and Cd. In addition, some hazardous metals (chromium (Cr), iron (Fe), As, Hg, and Pb) also were positively correlated with Zn and Cd. Interestingly, no correlations were found between forage pollutants and milk, suggesting that these pollutants have a low transfer rate to milk. The PCA results highlighted the predominant contribution of PAHs to the global variance in forage samples collected at different distances from industrial areas. In milk, the contributions of hazardous metals and PAHs were more balanced than in forages. Finally, when distances to potential pollution sources were included in the PCA of forage samples, a negative correlation was observed between the former and the concentrations of HW-PAHs, Cd, and Zn, suggesting that thermal power plants and steel factory emissions were the main sources of polluting forage in this area.

Effects of Cow's Milk Processing on MicroRNA Levels.

MicroRNAs (miRNAs) regulate gene expression and might resist adverse physicochemical conditions, which makes them potential biomarkers. They are being investigated as biomarkers of dairy production systems, based on the variations in their levels in raw milk depending on animal diet and management. Whether miRNA levels can serve as biomarkers for dairy products remains unclear, since technological or culinary treatments, such as fermentation, may alter their levels. Here, 10 cow dairy farms were sampled in Asturias (north-west Spain) and milk samples were subjected to microwave heating or used to produce yogurt or cheese. Total RNA was isolated from raw milk and three derived products, and levels of seven miRNAs, selected based on previous studies as possible milk production system biomarkers, were assessed by RT-qPCR. The treatments decreased levels of all miRNAs to some extent. These results also imply that cheesemaking increases the concentration of miRNAs in this product; raw milk and cheese supposedly may provide similar concentrations of miRNAs, higher than those of yogurt and microwaved milk. They also indicate that the content of certain miRNAs in raw milk cannot necessarily be extrapolated to other dairy products.

Hypertrophic osteopathy in a common dolphin (Delphinus delphis) with concurrent pulmonary Halocercus delphini infestation.

Dolphins are marine mammals that often live in coastal habitats. Common causes of severe skeletal disorders among wild dolphins are congenital vertebral anomalities, collisions with sea vessels, trauma, hunting-related injury, infectious diseases, environmental pollution, and tumors. A free-ranging male, 3-year-old common dolphin (Delphinus delphis) was found dead in the coast of Asturias in northern Spain. Postmortem examination revealed lordosis in the caudal vertebral column, while X-ray imaging and computer tomography showed well-organized palisade-like periosteal proliferation, appearing as florid-like accretions, along the spinous apophysis of 26 lumbar-caudal vertebrae. The transverse apophysis was affected on only a few caudal vertebrae. The cortical layer remained intact. Histology of vertebra tissue showed periosteal proliferation of cancellous bone. The animal was diagnosed with hypertrophic osteopathy. The lungs showed diffuse parasitic granulomatous bronchointerstitial pneumonia caused by Halocercus delphini, consolidation of the pulmonary tissue, congestion, and alveolar edema. The animal was also afflicted by parasitic granulomatous gastritis caused by Anisakis simplex sensu lato and tattoo skin disease. The dolphin suffered from hypertrophic osteopathy associated with pulmonary Halocercus delphini infestation. This syndrome, known as hypertrophic pulmonary osteopathy, has been described in diverse terrestrial mammals, including domestic animals, wildlife and humans, but not in dolphins. This case reports the first description of hypertrophic osteopathy associated to a pulmonary disorder in dolphin, and it provides insights into factors that can induce column malformation in dolphins, suggesting the importance of taking thoracic lesions into account during differential diagnosis.

Impact of Potentially Toxic Compounds in Cow Milk: How Industrial Activities Affect Animal Primary Productions.

Potentially toxic elements (PTEs) and polycyclic aromatic hydrocarbons (PAHs) frequently coexist in soils near industrial areas and sometimes in environmental compartments directly linked to feed (forage) and food (milk) production. However, the distribution of these pollutants along the dairy farm production chain is unclear. Here, we analyzed soil, forage, and milk samples from 16 livestock farms in Spain: several PTEs and PAHs were quantified. Farms were compared in terms of whether they were close to (<5 km) or far away from (>5 km) industrial areas. The results showed that PTEs and PAHs were enriched in the soils and forages from farms close to industrial areas, but not in the milk. In the soil, the maximum concentrations of PTEs reached 141, 46.1, 3.67, 6.11, and 138 mg kg-1 for chromium, arsenic, cadmium, mercury, and lead, respectively, while fluoranthene (172.8 µg kg-1) and benzo(b)fluoranthene (177.4 µg kg-1) were the most abundant PAHs. Principal component analysis of the soil PTEs suggested common pollution sources for iron, arsenic, and lead. In the forage, the maximum contents of chromium, arsenic, cadmium, mercury, and lead were 32.8, 7.87, 1.31, 0.47, and 7.85 mg kg-1, respectively. The PAH found in the highest concentration in the feed forage was pyrene (120 µg kg-1). In the milk, the maximum PTE levels were much lower than in the soil or the feed forages: 74.1, 16.1, 0.12, 0.28, and 2.7 µg kg-1 for chromium, arsenic, cadmium, mercury, and lead, respectively. Neither of the two milk samples exceeded the 20 µg kg-1 limit for lead set in EU 1881/2006. Pyrene was the most abundant PAH found in the milk (39.4 µg kg-1), while high molecular weight PAHs were not detected. For PTEs, the results showed that soil-forage transfer factors were higher than forage-milk ratios. Our results suggest that soils and forages around farms near industries, as well as the milk produced from those farms, have generally low levels of PTE and PAH contaminants.

Differences in the microRNAs Levels of Raw Milk from Dairy Cattle Raised under Extensive or Intensive Production Systems.

Studying microRNA (miRNAs) in certain agri-food products is attractive because (1) they have potential as biomarkers that may allow traceability and authentication of such products; and (2) they may reveal insights into the products' functional potential. The present study evaluated differences in miRNAs levels in fat and cellular fractions of tank milk collected from commercial farms which employ extensive or intensive dairy production systems. We first sequenced miRNAs in three milk samples from each production system, and then validated miRNAs whose levels in the cellular and fat fraction differed significantly between the two production systems. To accomplish this, we used quantitative PCR with both fractions of tank milk samples from another 20 commercial farms. Differences in miRNAs were identified in fat fractions: overall levels of miRNAs, and, specifically, the levels of bta-mir-215, were higher in intensive systems than in extensive systems. Bovine mRNA targets for bta-miR-215 and their pathway analysis were performed. While the causes of these miRNAs differences remain to be elucidated, our results suggest that the type of production system could affect miRNAs levels and potential functionality of agri-food products of animal origin.

Variation of miRNA Content in Cow Raw Milk Depending on the Dairy Production System.

Pasture-based milk presents several advantages over milk from intensive industrial farming in terms of human health, the environment, animal welfare, and social aspects. This highlights the need for reliable methods to differentiate milk according to its origin on the market. Here, we explored whether miRNA profiles could serve as a marker of milk production systems. We compared levels of previously described miRNAs in milk from four production systems (altogether 112 milk samples): grazing, zero grazing, grass silage or corn silage. Total RNA was extracted from the fat phase, and miRNAs levels were quantified by real-time quantitative PCR. The levels of the miRNAs bta-miR-155 and bta-miR-103 were higher in the grazing system than in corn silage farms. The levels of bta-miR-532, bta-miR-103 and bta-miR-7863 showed differences between different farm managements. The miRNAs bta-miR-155 and bta-miR-103 were predicted to participate in common functions related to fat metabolism and fatty acid elongation. All four differentially expressed miRNAs were predicted to participate in transport, cell differentiation, and metabolism. These results suggest that the dairy production system influences the levels of some miRNAs in milk fat, and that bta-miR-155 and bta-miR-103 may be potential biomarkers to identify milk from pasture-managed systems.

Wolf (Canis lupus) as canine adenovirus type 1 (CAdV-1) sentinel for the endangered cantabrian brown bear (Ursus arctos arctos).

Canine adenovirus type 1 (CAdV-1) causes infectious canine hepatitis (ICH) and has recently been described as a cause of death among endangered populations of European brown bear (Ursus arctos arctos) in the Cantabrian mountain range in Asturias, Spain. Sympatric wild and domestic carnivores can act as reservoirs of the virus and likely spread it into the environment and subsequently transmit it to brown bears. The present work investigates the prevalence and geo-temporal distribution of CAdV-1 among free-ranging wolves (Canis lupus) in Asturias from 2009 to 2018, during which three fatal cases of ICH were reported among brown bears in the region. A total of 149 wolves were analysed in this study, of which 21 (14%) were found to have CAdV-1 DNA based on real-time polymerase chain reaction (RT-PCR) of spleen samples. Prevalence of the virus was similar between males and females. All but one of the 20 CAdV-1-positive animals of estimable age were younger than 2 years, and only one of the 46 adult animals (>2 years) tested positive. Prevalence was highest in the western area of Asturias and during 2010 and 2011. Our results confirm that CAdV-1 is circulating in Asturian free-ranging wolves, supporting their possible role as virus reservoirs and sentinels in the region of this emerging disease in brown bears.

Goats Naturally Infected with the Spanish Goat Encephalitis Virus (SGEV): Pathological Features and An Outbreak.

In autumn 2011, a disease outbreak caused by Spanish goat encephalitis virus (SGEV) was reported in a herd of goats from Asturias (north-western Spain), expanding the known geographic distribution of tick-borne encephalitis in Europe. The virus was classified as a new subtype (subspecies) within the Louping-ill virus species of the mammalian tick-borne flavivirus group. The aims of the present study were to describe the pathology in goats naturally infected with SGEV, as well as discuss the pathogenesis of the disease in that outbreak. A total of 22/85 (25.88%) goats (20 adults and 2 kids) died between October 2011 and June 2012, showing neurological clinical signs. Over three years, the mortality rate in the herd reached 100%. Neuropathological lesions caused by SGEV were severe and widespread throughout the central nervous system but were more severe and numerous in the proximal cervical spinal cord, medulla oblongata, pons and cerebellar cortex. They consisted of neuron necrosis, neuronophagia, mononuclear inflammatory cell perivascular cuffs (lymphocytes, plasma cells and macrophages) and gliosis. The distribution of viral antigens was restricted to the cytoplasm of neurons in several brain areas but not associated with inflammatory foci nor inflammatory cells. SGEV should be considered a significant pathogen of goats that results in severe neurological clinical disease and high mortality.

Mortality Causes in Free-Ranging Eurasian Brown Bears (Ursus arctos arctos) in Spain 1998-2018.

This work summarizes the mortality cases of twenty-five free-ranging Eurasian wild brown bears (Ursus arctos arctos) from the Cantabrian mountain range submitted for necropsy in Asturias and Castilla y León (northwestern Spain) from 1998 to 2018. Mortality cases were classified both caused by (i) "non-human intervention" or "human intervention" causes and based on (ii) "non-infectious" or "infectious" etiology. In four cases (16%) it was not possible to determine the cause of death due to the inadequate preservation of collected specimens or insufficient tissue availability. Based on "non-human intervention" or "human intervention" causes, fourteen of the 21 (66.7%) brown bears died as a consequence of "non-human intervention" due to traumatic lesions (fights, unknown traumas or infanticide), infectious canine hepatitis, neoplasia or mushroom poisoning. In contrast, seven (33.3%) brown bears died by "human intervention" due to illegal hunting (shooting or snare), handling (during transit in an attempt to reintroduce a bear back into the wild) or strychnine poisoning. Based on "non-infectious" or "infectious" etiology, twelve of the 21 (57.1%) brown bears died due to "non-infectious" causes, namely traumatic lesions such as shooting, snare, fighting or infanticide, handling, strychnine poisoning, mushroom poisoning or neoplasia. The remaining nine (42.9%) animals died due to "infectious" diseases which included gangrenous myositis, infectious canine hepatitis or septicemia. In six of those cases traumatic lesions caused by non-human or human activities were complicated with bacterial infection (clostridiosis and septicemia) which finally caused the death of those animals. Additionally, exertional myopathy was observed in the handled animal and in one bear found in a snare. In a free-ranging population of Eurasian brown bear from the Cantabrian mountain range, main causes of death are attributed to non-human related traumatic lesions and infectious diseases (primary developed such as infectious canine hepatitis or secondary developed such as clostridiosis or septicemia) which is in contrast to previously reported data for other bear populations. These data are valuable and may help in the conservation and management of this recovering population.

The Mode of Grass Supply to Dairy Cows Impacts on Fatty Acid and Antioxidant Profile of Milk.

The optimization of milk production includes a rational use of forages, respect for the environment and offers the best quality to consumers. Milk production based on grass and forages produces healthier milk and it is widely spread throughout the Atlantic arc to maximize milk yield per hectare. However, the mode of offering the grass can have a major influence on milk composition. The aim of this study was to evaluate the effect of grass supply mode (grazing, zero-grazing or ensiling) on dairy cows' performance, with particular reference to fatty acids and fat-soluble antioxidants concentration. A three by three Latin square experiment was performed with 18 dairy cows. Experimental treatments consisted of exclusive feeding with grass silage and zero-grazing, both offered ad libitum indoors, or grazing for 24 h. The results showed that grazing cows had a higher dry matter intake and greater milk yield than cows feeding on grass silage and zero-grazing, as well as higher concentrations of protein, lactose, nonfat-solids and urea in milk than housed cows. Milk fat from grazing cows had a higher proportion of unsaturated fatty acids than from cows feeding on grass silage and zero-grazing, with significant differences in the proportion of vaccenic and rumenic acids. The 18:1 trans-11 to 18:1 trans-10 ratio is proposed as biomarker to identify the milk produced from the management system of grazing cattle. Milk from grazing cows had a greater proportion of lutein than cows eating grass silage, with the zero-grazing system having intermediate values. In conclusion, the mode of grass supply affects fatty acid and antioxidant profiles of milk.

Canine adenovirus type 1 (CAdV-1) in free-ranging European brown bear (Ursus arctos arctos): A threat for Cantabrian population?

Canine adenovirus type 1 (CAdV-1) is responsible for infectious canine hepatitis. The disease has been described in captive American black bear (Ursus americanus) and European brown bear (Ursus arctos arctos), with just one recently reported case in a cub of a free-ranging brown bear (Ursus arctos horribilis) from Alaska. The aim of this work is to summarize findings related to presence and associated mortality of CAdV-1 in 21 free-ranging Cantabrian brown bears (Ursus arctos arctos) submitted to necropsy in Asturias and Castilla y León (northwestern Spain) from 1998 to 2018. On the basis of the anatomopathological findings and laboratory results three free-ranging brown bears died due to infectious canine hepatitis, which is to our knowledge the first description of death due to this disease in free-ranging bears in Europe. Gross lesions consisted of petechial haemorrhages and congestion in different internal organs, haemorrhagic fluid in internal cavities, friable and yellowish liver and thickening of gall bladder. Microscopic lesions were observed mainly in liver, kidney and brain and consisted of multifocal necrosis of cells with presence of basophilic intranuclear inclusions. Immunohistochemical (IHC) and real-time polymerase chain reaction (qPCR) techniques were used to assess the presence of CAdV-1 in paraffin-embedded liver samples. Viral antigens were detected by IHC labelling within hepatocytes and Küppfer cells in the three animals. The presence of viral DNA was confirmed by qPCR in one of them. In order to evaluate the circulation of CAdV-1 in brown bears, a retrospective study was performed using both IHC and qPCR techniques in 11 and 12 additional brown bears, respectively. An extra brown bear was found positive by IHC. This study shows that CAdV-1 surveillance of brown bears and sympatric carnivores should be considered as major concern for the monitoring the population evolution throughout time in this endangered species.

A one-step TaqMan real-time qRT-PCR assay for the specific detection and quantitation of the Spanish goat encephalitis virus (SGEV).

Louping ill-like virus (LI) has been recently detected in two different locations in the north of Spain and separated by only around 400 km. Using molecular approaches, the viruses causing both outbreaks have been shown to be different to LI virus, but also different to each other. They have been called SSEV (Spanish sheep encephalitis virus) and SGEV (Spanish goat encephalitis virus) taking into account the species from which they were isolated. The aim of this paper was to design a quantitative TaqMan real-time RT-PCR protocol, for the specific diagnostic and quantitation of SGEV. Linearity, efficiency and dynamic range as well as reproducibility and specificity of the method has been tested and established. The method has proved to be valid for the specific detection and viral load quantitation of SGEV genome in virus isolates and tissue samples from infected animals. This assay will be a useful analytical tool in early diagnosis and epidemiological surveys.

Spatial relationships between Eurasian badgers (Meles meles) and cattle infected with Mycobacterium bovis in Northern Spain.

Recent studies suggest that badgers may be a potential reservoir of Mycobacterium bovis infection for cattle in Northern Spain. The objective of this study was to investigate potential epidemiological links between cattle and badgers. Culture and molecular typing data were available for cattle culled during the national tuberculosis (TB) eradication campaigns between 2008 and 2012, as well as from 171 necropsied badgers and 60 live animals trapped and examined over the same time period. Mycobacterium tuberculosis complex strains were isolated from pooled tissues of 14 (8.2%) necropsied badgers, of which 11 were identified as M. bovis: six different spoligotypes of M. bovis were subsequently identified. In two geographical locations where these isolates were shared between cattle and badgers, infected cattle herds and badgers lived in close contact. Although it remains unclear if badgers are a maintenance or spill-over host of M. bovis in this setting, it would appear prudent to have precautionary measures in place to reduce contact between cattle and badgers.

Authentication of male beef by multiplex fast real-time PCR.

Meat is a significant source of high-quality nutrients, which are very important in the diet. Among meat products, one of the most prized is bovine meat, of which male beef has been designated to be of a higher quality. However, because of its similarity with female beef, deliberate or unintentional substitutions can occur. To avoid this, methodology based on the fast real-time polymerase chain reaction has been developed to authenticate the species and gender origin of beef. This technique consists of two polymerase chain reactions: one bovine-specific reaction and another Y-chromosome-specific multiplex reaction. This methodology has been validated for all kinds of beef products, including those subjected to intensive processing treatments, and it has subsequently been applied to 10 commercial samples labelled as ox to determine whether they are properly labelled. This assay has been shown to have high specificity, sensitivity and rapidity, with the potential to be a powerful tool in enforcing food labelling regulations.

Louping ill in goats, Spain, 2011.

Although louping ill affects mainly sheep, a 2011 outbreak in northern Spain occurred among goats. Histopathologic lesions and molecular genetics identified a new strain of louping ill virus, 94% identical to the strain from Britain. Surveillance is needed to minimize risk to domestic and wildlife species and humans.

Prion protein gene polymorphism in four West African sheep populations.

A total of 162 individuals, belonging to three Burkinabé and one Niger sheep populations, were analysed for prion protein (PrP) gene polymorphism at codons 136, 154 and 171. The ARQ allele was the most frequent in both the Burkinabé (86.7%) and the Niger (67.5%) sheep populations. The highly sensitive allele VRQ was not found in the sampled individuals. The highly resistant ARR allele was in very low frequency in the Burkina-Sahel (4.4%) and Mossi (3.2%) populations and was not present in the Djallonké and Touareg populations. Only 4 out of 15 possible PrP genotypes were identified in the sampled individuals. No favourable ARR/ARR genotypes were found in either of the breeds. Sequencing a subgroup of the samples allowed the identification of other five polymorphisms on the PrP gene sequence at codons 116, 138, 151, 237 and 240. The very low frequency of the ARR allele in the West African sheep should dissuade the implementation of a preventive selection programme aimed to increase resistance to scrapie, to avoid an extreme erosion of the genetic stock.

Zebu cattle are an exclusive legacy of the South Asia neolithic.

Animal domestication was a major step forward in human prehistory, contributing to the emergence of more complex societies. At the time of the Neolithic transition, zebu cattle (Bos indicus) were probably the most abundant and important domestic livestock species in Southern Asia. Although archaeological evidence points toward the domestication of zebu cattle within the Indian subcontinent, the exact geographic origins and phylogenetic history of zebu cattle remains uncertain. Here, we report evidence from 844 zebu mitochondrial DNA (mtDNA) sequences surveyed from 19 Asiatic countries comprising 8 regional groups, which identify 2 distinct mitochondrial haplogroups, termed I1 and I2. The marked increase in nucleotide diversity (P < 0.001) for both the I1 and I2 haplogroups within the northern part of the Indian subcontinent is consistent with an origin for all domestic zebu in this area. For haplogroup I1, genetic diversity was highest within the Indus Valley among the three hypothesized domestication centers (Indus Valley, Ganges, and South India). These data support the Indus Valley as the most likely center of origin for the I1 haplogroup and a primary center of zebu domestication. However, for the I2 haplogroup, a complex pattern of diversity is detected, preventing the unambiguous pinpointing of the exact place of origin for this zebu maternal lineage. Our findings are discussed with respect to the archaeological record for zebu domestication within the Indian subcontinent.

Identifying the most suitable endogenous control for determining gene expression in hearts from organ donors.

BACKGROUND: Quantitative real-time reverse transcription PCR (qRT-PCR) is a useful tool for assessing gene expression in different tissues, but the choice of adequate controls is critical to normalise the results, thereby avoiding differences and maximizing sensitivity and accuracy. So far, many genes have been used as a single reference gene, without having previously verified their value as controls. This practice can lead to incorrect conclusions and recent evidence indicates a need to use the geometric mean of data from several control genes. Here, we identified an appropriate set of genes to be used as an endogenous reference for quantifying gene expression in human heart tissue. RESULTS: Our findings indicate that out of ten commonly used reference genes (GADPH, PPIA, ACTB, YWHAZ, RRN18S, B2M, UBC, TBP, RPLP and HPRT), PPIA, RPLP and GADPH show the most stable gene transcription levels in left ventricle specimens obtained from organ donors, as assessed using geNorm and Normfinder software. The expression of TBP was found to be highly regulated. CONCLUSION: We propose the use of PPIA, RPLP and GADPH as reference genes for the accurate normalisation of qRT-PCR performed on heart tissue. TBP should not be used as a control in this type of tissue.

Genetic variability and differentiation in Spanish roe deer (Capreolus capreolus): A phylogeographic reassessment within the European framework.

Here we study 109 Iberian roe deer individuals corresponding to 9 Spanish populations. Individuals were sampled at locations that are expected to have acted as refugia for the species during the 20th century. Samples were analyzed for a 436 bp mtDNA fragment from the HVI region of mitochondrial DNA and 10 microsatellites. The 109 sequences gave 31 different haplotypes that enabled identification of a new haplogroup (mainly present in Northwestern Iberian populations and representing roughly a third of our samples) that is not present in other European roe deer populations. Using microsatellites, correspondence analysis and molecular coancestry information revealed high molecular differentiation among Northwestern and Central-Southern Spain roe deer populations. Both sequence and microsatellite analysis reveal that the Spanish roe deer populations are genetically heterogeneous and have high genetic structure clearly separating the Pyrenean-reintroduced populations and two main areas for the species in Spain (Northwestern and Central-Southern) coinciding with the two main areas acting as refugia for the majority of the mammal species during glaciations. The implications of the obtained information as regards the phylogeography of the species are discussed together with suggestions as to appropriate strategies for the conservation and management of populations.