Genomic analysis of inbreeding level, kinship and breed relationships in Creole cattle from South America.

The conservation of animal genetic resources refers to measures taken to prevent the loss of genetic diversity in livestock populations, including the protection of breeds from extinction. Creole cattle populations have suffered a drastic reduction in recent decades owing to absorbent crosses or replacement with commercial breeds of European or Indian origin. Genetic characterization can serve as a source of information for conservation strategies to maintain genetic variation. The objective of this work was to evaluate the levels of inbreeding and kinship through the use of genomic information. A total of 903 DNAs from 13 cattle populations from Argentina, Bolivia and Uruguay were genotyped using an SNP panel of 48 K. Also, a dataset of 76 K SNPs from Peruvian Creole was included. Two inbreeding indices (FROH and Fhat2) and kinship relationships were calculated. In addition, effective population size (Ne), linkage disequilibrium, population composition and phylogenetic relationships were estimated. In Creole cattle, FROH ranged from 0.14 to 0.03, and Fhat2 was close to zero. The inferred Ne trends exhibited a decline toward the present for all populations, whereas Creole cattle presented a lower magnitude of Ne than foreign breeds. Cluster analysis clearly differentiated the taurine and Zebu components (K2) and showed that Bolivian Creole cattle presented Zebu gene introgression. Despite the population reduction, Creole populations did not present extreme values of consanguinity and kinship and maintain high levels of genetic diversity. The information obtained in this work may be useful for planning conservation programmes for these valuable local animal genetic resources.

Polledness in Argentinean Creole cattle, five centuries surviving.

Polledness has been shown to have autosomal Mendelian inheritance, with the polled locus being dominant to the horned locus. This trait was mapped to the BTA1 centromeric end in several breeds. One of the distinctive attributes of Creole cattle, such as the Argentinean Creole, is the presence of long, lyre-shaped horns. However, polled native animals were reported before the introduction of modern selected European breeds. Here, we studied the origin of the polled mutation, either independent or introgressed, in a Creole line from the Creole cattle founder group at the IIACS-INTA Leales Experimental Station (northwest Argentina). The study sample (65 animals: 26 horned and 39 polled) was genotyped using high-density SNP microarrays and three previously reported genetic markers (P202 ID , P80kb ID and PG ). A genome-wide association study, selection signatures, linkage disequilibrium analysis and copy number variations were used to detect the responsible region and the segregating haplotypes/alleles. The interval mapped in the Leales herd (1.23-2.13 Mb) overlapped with the region previously reported in several European cattle breeds, suggesting that the same locus could be segregating in this population. The previously reported variants PF and PG were not detected, thus dismissing the Holstein-Friesian and Nellore origins of the polled phenotype in this native breed. Conversely, the presence of the Celtic variant PC suggests an almost complete co-segregation. The cluster analysis rejected the hypothesis of recent introgression, which is compatible with the historical record of polled Creole cattle in northwest Argentina.

Evidence of positive selection towards Zebuine haplotypes in the BoLA region of Brangus cattle.

The Brangus breed was developed to combine the superior characteristics of both of its founder breeds, Angus and Brahman. It combines the high adaptability to tropical and subtropical environments, disease resistance, and overall hardiness of Zebu cattle with the reproductive potential and carcass quality of Angus. It is known that the major histocompatibility complex (MHC, also known as bovine leucocyte antigen: BoLA), located on chromosome 23, encodes several genes involved in the adaptive immune response and may be responsible for adaptation to harsh environments. The objective of this work was to evaluate whether the local breed ancestry percentages in the BoLA locus of a Brangus population diverged from the estimated genome-wide proportions and to identify signatures of positive selection in this genomic region. For this, 167 animals (100 Brangus, 45 Angus and 22 Brahman) were genotyped using a high-density single nucleotide polymorphism array. The local ancestry analysis showed that more than half of the haplotypes (55.0%) shared a Brahman origin. This value was significantly different from the global genome-wide proportion estimated by cluster analysis (34.7% Brahman), and the proportion expected by pedigree (37.5% Brahman). The analysis of selection signatures by genetic differentiation (F st ) and extended haplotype homozygosity-based methods (iHS and Rsb) revealed 10 and seven candidate regions, respectively. The analysis of the genes located within these candidate regions showed mainly genes involved in immune response-related pathway, while other genes and pathways were also observed (cell surface signalling pathways, membrane proteins and ion-binding proteins). Our results suggest that the BoLA region of Brangus cattle may have been enriched with Brahman haplotypes as a consequence of selection processes to promote adaptation to subtropical environments.

Nonsynonymous changes of equine lentivirus receptor-1 (ELR1) gene in amino acids involved in the interaction with equine infectious anemia virus (EIAV).

Equine lentivirus receptor-1 (ELR1) has been characterized as the specific functional receptor that mediates equine infectious anemia virus (EIAV) entrance to horse macrophages. This receptor is tumor necrosis factor receptor superfamily member 14 (TNFRSF14). The aim of this study was to investigate the occurrence of allelic variants in the coding sequence of equine TNFRSF14 gene by screening for single-nucleotide polymorphisms (SNPs) in different equine populations. Forty seven horse samples were randomly selected from a reservoir of EIAV-seropositive and seronegative samples collected from different outbreaks and regions of Argentina. DNA samples were scanned via PCR and direct sequencing of exon 3 and exon 5 of TNFRSF14 gene. A total of 21 SNPs were identified, of which 11 were located in coding sequences. Within exon 5, four SNPs caused nonsynonymous substitutions, while two other SNPs caused synonymous substitutions in crucial residues (Ser112 and Thr114) implicated in the interaction with EIAV. Despite some of exon 5 variants occurred exclusively in EIAV-positive or EIAV-negative horses, critical residues for the function of the mature protein were conserved, accounting for selective pressures in favor of preserving the specific function of TNFRSF members and the host immune response. To our knowledge, this is the first report of the existence of allelic variations involving some crucial amino acid residues in horse ELR1. Further, it could be an initial step to test the possible functional relevance and relationship of these variants with EIAV infection and disease progression as well as to develop preventive strategies.

DNA profile of dog feces as evidence to solve a homicide.

Dog fecal samples were collected at the crime scene and from the shoes of the suspect to see whether they could be linked. DNA was genotyped using a 145bp fragment containing a 60bp hotspot region of the mitochondrial DNA (mtDNA) control region. Once the species origin was identified, sequences were aligned with the 23 canine haplotypes defined, showing that evidence and reference had 100% identity with haplotype 5. The frequency of haplotype 5 and the exclusion power of the reference population were 0.056 and 0.89, respectively. The forensic index showed that it was 20 times more likely that the evidence belonged to the reference dog than to some other unknown animal. The results support that the mtDNA hypervariable region 1 (HV1) is a good alternative for typing in trace or degraded casework samples when the STR panel fails, and demonstrate the utility of domestic animal samples to give additional information to solve human legal cases.

The importance of having zinc during in vitro maturation of cattle cumulus-oocyte complex: role of cumulus cells.

The aim of this study was to investigate the influence of zinc (Zn) on the health of cumulus-oocyte complex (COC) during in vitro maturation (IVM). Experiments were designed to evaluate the effect of Zn added to IVM medium on: DNA integrity, apoptosis, cumulus expansion and superoxide dismutase (SOD) activity of cumulus cells (CC). Also, role of CC on Zn transport during IVM was evaluated on oocyte developmental capacity. DNA damage and early apoptosis were higher in CC matured with 0 µg/ml Zn compared with 0.7, 1.1 and 1.5 µg/ml Zn (p < 0.05). Cumulus expansion did not show differences in COC matured with or without Zn supplementation (p > 0.05). Superoxide dismutase activity was higher in COC matured with 1.5 µg/ml Zn than with 0 µg/ml Zn (p < 0.05). Cleavage and blastocyst rates were recorded after IVM in three maturation systems: intact COCs, denuded oocytes with cumulus cells monolayer (DO + CC) and denuded oocytes (DO). Cleavage rates were similar when COC, DO + CC or DO were matured with 1.5 µg/ml Zn compared with control group (p > 0.05). Blastocyst rates were significantly higher in COC than in DO + CC and DO with the addition of 1.5 µg/ml Zn during IVM (p < 0.01). Blastocyst quality was enhanced in COC and DO + CC compared with DO when Zn was added to IVM medium (p < 0.001). The results of this study indicate that Zn supplementation to IVM medium (i) decreased DNA damage and apoptosis in CC; (ii) increased SOD activity in CC; (iii) did not modify cumulus expansion and cleavage rates after in vitro fertilization; (iv) improved subsequent embryo development up to blastocyst stage; and (v) enhanced blastocyst quality when CC were present either in intact COC or in coculture during IVM.

Foreign meat identification by DNA breed assignment for the Chinese market.

Methods for individual identification are usually employed for traceability, whereas breed identification is useful to detect commercial frauds. In this study, Chinese Yellow Cattle (CYC) samples plus data from six Bos taurus breeds, two Bos indicus breeds, and one composite breed were used to develop an allocation test based on 22 microsatellites. The test allowed discriminating all foreign breeds from the CYC, although some CYC individuals were wrongly allocated as Limousin or Holstein, probably due to the recent introduction of these breeds into China. In addition, CYC evidenced a previously reported Zebu cline (south-north) and a possible structure within the B. taurus component that should be confirmed. An independent test performed with meat samples of unknown breed origin from Argentina allocated 92% of them to either Angus, Hereford, or their crossbreed, but none was identified as CYC. We conclude that the test is a suitable tool to certify meat of foreign breed origin and to detect adulterations of CYC beef labeled as imported meat.

ELA-DRA polymorphisms are not associated with Equine Arteritis Virus infection in horses from Argentina.

Polymorphisms at Major Histocompatibility Complex (MHC) genes have been associated with resistance/susceptibility to infectious diseases in domestic animals. The aim of this investigation was to evaluate whether polymorphisms of the DRA gene the Equine Lymphocyte Antigen is associated with susceptibility to Equine Arteritis Virus (EAV) infection in horses in Argentina. The equine DRA gene was screened for polymorphisms using Pyrosequencing® Technology which allowed the detection of three ELA-DRA exon 2 alleles. Neither allele frequencies nor genotypic differentiation exhibited any statistically significant (P-values=0.788 and 0.745) differences between the EAV-infected and no-infected horses. Fisher's exact test and OR calculations did not show any significant association. As a consequence, no association could be established between the serological condition and ELA-DRA.

Eca20 microsatellite polymorphisms in equine viral arteritis-infected horses from Argentina.

We investigated the association of equine arteritis virus (EAV) infection and three short tandem repeat (STR) polymorphisms located within or in close proximity to equine lymphocyte antigen (ELA) region. We used a case-control design as a first approach before proceeding to select candidate genes. One hundred and sixty-five Silla Argentino horses were taken in 2002 from positive serological detections of EAV in Argentina, to determine whether STR genotypes were correlated to genetic susceptibility to EVA. Allele frequency distribution did not show significant differences between both groups (P = 0.0781). However, in particular alleles, Fisher exact test and odds ratio calculations showed significant values >1 for TKY08 and LEX52, and <1 for UM011, TKY08, LEX52 and VHL20. Interestingly, TKY08 STR is located in ELA class I region.

History and selection imprinting on genetic relationships among bovine breeds analyzed trough five genes related with marbling.

Many candidate genes have been suggested as responsible for marbling in beef cattle, for instance diacylglycerol O-acyltransferase 1, thyroglobulin, growth hormone, leptin and stearoyl CoA desaturase. The objective of the present work was to evaluate the polymorphisms of five SNPs of these candidate genes in 389 animals of 18 Bos Taurus and Bos indicus breeds. The obtained results were compared with ones previously obtained with STRs and loci related to milk production in these populations. Moreover we analyzed whether the phylogenies reconstructed using SNPs associated with marbling resulted in the known tree topology. The tree constructed with UPGMA, using genetic distance D(A), exhibit a topology partially consistent with the historical origin of breeds. The result observed in the Correspondence Analysis coincided with the topology of the UPGMA tree. This work allowed us to evaluate the five SNPs genetic diversity and to demonstrate that the grouping of the breeds may be the result of its history, selection process, or both at once.

Characterization and validation of bovine gonadotripin releasing hormone receptor (GNRHR) polymorphisms.

Gonadotropin releasing hormone and its receptor (GNRHR) play a critical role in sexual differentiation and reproduction. Available evidence shows a strong genetic component in the timing of puberty. In bovines, there are significant differences within and among beef breeds in the time when bulls reach puberty. Despite its economic importance, there are not many SNPs or genetic markers associated with this characteristic. The aims of the study were to identify DNA polymorphism in the bovine GNRHR by re-sequencing analysis, determine haplotype phases, and perform a population study in a selected tag SNP in six breeds. Eight SNPs were detected, including: one in the Upstream Regulatory Region (URR), five in the coding regions, and two in non-coding regions. This polymorphism level corresponds to one variant every 249.4bp and a global nucleotide diversity of 0.385. Two haplogroups comprising nine haplotypes and two linkage blocks were detected. Despite 5 tag SNPs were required to capture all variability, just one SNP allowed to define both haplogroups, and only two SNPs were needed to differentiate the most common haplotypes. An additional taq SNP was necessary to identify both URR variants. Allele-frequency analysis of a selected taq SNP among breeds showed a geographical cline. European Bos taurus breeds had lower frequencies of the C allele than B. indicus type cattle, while Creole cattle and Wagyu breeds had intermediate frequency. There was a significant correlation between frequency profile and timing of puberty among the studied breeds, which seems to suggest that genetic variation within bovine GNRHR gene could explain at least part of the reported variability.

Effect of consanguinity on Argentinean Angus beef DNA traceability.

Since the 1990s several authors have envisaged the use of DNA to certify meat origin. Two major parameters must be assessed before a DNA based traceability protocol can be implemented in the food chain: (i) the information content of a DNA marker set in a specific livestock breed or group of breeds; (ii) the minimum number of DNA markers needed to obtain a statistically acceptable match probability. The objective of the present work was to establish the effect of different levels of inbreeding in the matching efficiency, and the minimum number of microsatellite markers needed, in a DNA based meat traceability program, starting from an 11-microsatellite marker panel. Samples were obtained from beef production farms in South America, where animals are typically bred under pasture-based extensive conditions. Three groups of animals with different consanguinity rates were sampled. Exclusion power (Q) was higher than 0.999998 and match probability lower than 3.01E-08, for the whole set of markers within each group. Both values were affected by consanguinity. To reach a two mismatch criteria exclusion power (Q(2)) of 99.99, six markers were needed in unrelated animals whereas seven markers were needed in related animals. To reach Q(2)=99.9999, 8 and 10 microsatellite markers, respectively, were needed. In general, one or two more microsatellite markers were needed to identify consanguineous animals. This study proved the DNA marker set used to be suitable for the identification of the meat from all slaughtered animals in Argentina, per week, month, and year.

Development of an ELA-DRA gene typing method based on pyrosequencing technology.

The polymorphism of equine lymphocyte antigen (ELA) class II DRA gene had been detected by polymerase chain reaction-single-strand conformational polymorphism (PCR-SSCP) and reference strand-mediated conformation analysis. These methodologies allowed to identify 11 ELA-DRA exon 2 sequences, three of which are widely distributed among domestic horse breeds. Herein, we describe the development of a pyrosequencing-based method applicable to ELA-DRA typing, by screening samples from eight different horse breeds previously typed by PCR-SSCP. This sequence-based method would be useful in high-throughput genotyping of major histocompatibility complex genes in horses and other animal species, making this system interesting as a rapid screening method for animal genotyping of immune-related genes.

African matrilineages in American Creole cattle: evidence of two independent continental sources.

In order to clarify the historical origin and phylogeographic affinities of Creole cattle matrilineages throughout the American continent, we analysed published D-loop mtDNA sequences (n = 454) from Creole, Iberian and African cattle breeds. The Western European T3 haplogroup was the most common in American Creole cattle (63.6%), followed by the African T1 (32.4%) and the Near Eastern T2 haplogroups (4%). None of the sequences were found in Bos indicus types. Within the African T1 haplogroup there were two subclades, T1a and T1*, whose geographic distribution in America was clearly disjointed. T1a is a highly divergent clade originally reported for Creole cattle from Brazil and the Lesser Antilles, but whose geographic distribution in Africa remains unknown. In contrast, lineages attributable to T1* are restricted in America to the region colonized by the Spaniards. We propose a new hypothesis for the origins of Creole cattle that summarizes all previously published historical and genetic data. While the African T1* fraction in Creole cattle may have arrived in America through the Iberian breeds, the divergent T1a lineages may have been introduced by Portuguese and other European crowns from some unknown, not-yet-sampled African location. Additional molecular studies will be required for pinpointing the specific African regional source.

New polymorphisms for the BoLA-DRB3 upstream regulatory region.

Two new alleles, named BoLA-DRB3-P*06 and BoLA-DRB3-P*07, have been identified for the upstream regulatory region of the BoLA-DRB3 gene. The 228-bp nucleotide sequences of the promoter comprising the W, X, Y, CAAT and TATA regulatory boxes were analysed. The BoLA-DRB3-P*06 exhibits one insertion between the W and X boxes, and one transition between the X and Y boxes. On the other hand, the BoLA-DRB3-P*07 showed one insertion in the X box.

Polymorphisms of the upstream regulatory region of the major histocompatibility complex DRB genes in domestic horses.

Sequence information was obtained on the variation of the ELA-DRB upstream regulatory region (URR) after polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) cloning and sequencing of approximately 220 bp upstream of the first exon of horse DRB genes. The sequence of the proximal URR of equine DRB is composed of highly conserved sequence motifs, showing the presence of the W, X, Y, CAAT and TATA conserved boxes of major histocompatibility complex (MHC) class II promoters. Five different polymorphic horse DRB promoter sequences were detected in five horse breeds. The results demonstrate the existence of polymorphism in the nucleotide sequences of the ELA-DRB URR, located in the functionally important conserved consensus sequences, the X2 box, the Y box and the TATA box, while conservation were observed in X1 and CAAT boxes. The nucleotide diversity among horse URRs was intermediate between that seen within human and mouse DRB promoters, suggesting the existence of another important source of variability in ELA-DRB genes. In addition, phylogenetic comparisons, identity analysis and sequence organization suggested that the reported sequences would correspond to an expressed ELA-DRB locus. However, further information about the functional significance of these promoter polymorphisms will probably be acquired through expression studies on the different sequences.

Polymorphism in the bovine BOLA-DRB3 upstream regulatory regions detected through PCR-SSCP and DNA sequencing.

In the present work, we describe through polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) and DNA sequencing the polymorphism within the URR-BoLA-DRB3 in 15 cattle breeds. In total, seven PCR-SSCP defined alleles were detected. The alignment of studied sequences showed six polymorphic sites (four transitions, one transversion and one deletion) in the interconsensus regions of the BoLA-DRB3 upstream regulatory region (URR), while the consensus boxes were invariant. Five out of six detected polymorphic sites were of one nucleotide substitution in the interconsensus regions. It is expected that these mutations do not affect significantly the level of expression. In contrast, the deletion observed in the sequence between CCAAT and TATA boxes could have some effect on affinity interactions between the promoter region and the transcription factors. The URR-BoLA-DRB3 DNA analyzed sequences showed moderate level of nucleotide diversity, high level of identity among them and were grouped in the same clade in the phylogenetic tree. In addition, the phylogenetic tree, the similarity analysis and the sequence structure confirmed that the fragment analyzed in this study corresponds to the URR-BoLA-DRB3. The functional role of the observed polymorphic sites among the regulatory motifs in bovine needs to be analyzed and confirmed by means of gene expression assays.

Analysis of ELA-DQB exon 2 polymorphism in Argentine Creole horses by PCR-RFLP and PCR-SSCP.

The second exon of equine leucocyte antigen (ELA)-DQB genes was amplified from genomic DNA of 32 Argentine Creole horses by PCR. Amplified DNA was analysed by PCR-restriction fragment length polymorphism (RFLP) and PCR-single-strand conformation polymorphism (SSCP). The PCR-RFLP analysis revealed two HaeIII patterns, four RsaI patterns, five MspI patterns and two HinfI patterns. EcoRI showed no variation in the analysed sample. Additional patterns that did not account for known exon 2 DNA sequences were observed, suggesting the existence of novel ELA-DQB alleles. PCR-SSCP analysis exhibited seven different band patterns, and the number of bands per animal ranged from four to nine. Both methods indicated that at least two DQB genes are present. The presence of more than two alleles in each animal showed that the primers employed in this work are not specific for a unique DQB locus. The improvement of this PCR-RFLP method should provide a simple and rapid technique for an accurate definition of ELA-DQB typing in horses.

Analysis of genetic diversity and population structure in Argentine and Bolivian Creole cattle using five loci related to milk production

Data from five protein-coding loci related to dairy production were used to study the genetic diversity and population structure of Argentine and Bolivian Creole cattle breeds. Genomic DNA was extracted from blood samples of six Creole cattle breeds: Argentine (n = 230), Patagonian (n = 25); "Saavedreño" (n = 140), "Chaqueño Boliviano" (n = 30), "Yacumeño" (n = 27), and "Chusco" (n = 11). kappa-casein, beta-lactoglobulin, growth hormone and prolactin were measured by PCR-RFLP, while alphaS1-casein was typed by PCR-ASO. The results are discussed, focusing on: historical origin, recent differentiation and selection events, Zebu gene introgression, and population structure. This work shows that: (i) For the studied genes, the observed gene frequency profiles of Argentine and Bolivian Creole cattle breeds were close to the data reported for Iberian breeds and for other South-American Creole cattle breeds which are historically related; (ii) although Zebu gene introgression has been reported at the studied loci, these breeds seem to be far from the Zebu gene frequency profiles; and (iii) the Argentine and Bolivian Creole cattle showed significant levels of subdivision, but each population has maintained its degree of genetic variability