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.

Analysis of Doberman Pinscher and Toy Poodle samples with targeted next-generation sequencing.

Next-generation sequencing (NGS) technologies have enabled the identification of many causal variants of genetic disorders, the development of parentage tests and the analysis of multiple traits in domestic animals. In this study, we evaluated the performance of a Canine Targeted Genotyping-by-Sequencing (GBS) custom panel (Thermo Fisher Scientific, Waltham, Ma, USA) in a cohort of 95 dog DNA samples, comprising 76 Doberman Pinschers and 19 Toy Poodles from Argentina. The used panel included 383 targets (228 parentage SNVs, 137 genetic disorder markers and 18 trait markers). While paternity analysis showed correct duo (97.4%; LOD > 2.98E+13) and trio (100%; LOD > 2.20E+15) parentage assignment, the panel resulted still insufficient for excluding close relatives in inbred populations. In this sense, close relatives were wrongly assigned as parents in 12.6% of duos and 0.3% of trios. We detected 17 polymorphic markers (genetic disorders, n = 4; hair type, n = 3; coat color, n = 10) and estimated their allele frequencies in the studied breeds. The accuracy of targeted GBS results were evaluated for three markers that were associated with Progressive rod-cone degeneration, von Willebrand disease type 1 and dilated cardiomyopathy by pyrosequencing and Sanger sequencing genotyping, showing 94-100% concordance among assays. The targeted GBS custom panel resulted cost-effective strategy to study the prevalence of genetic disorders and traits in a large number of samples and to analyze genetic interactions between previously reported variants. Once assays based on AgriSeq technology were standardized, their uses are a good strategy for large-scale routine genetic evaluation of animal populations.

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.

The copper transporter (SLC31A1/CTR1) is expressed in bovine spermatozoa and oocytes: Copper in IVF medium improves sperm quality.

Adequate dietary intake of copper (Cu) is required for normal reproductive performance in cattle. The objective of this study was to investigate the pregnancy rates from cattle with deficient, marginal and adequate Cu plasma concentration at the beginning of artificial insemination protocol. Moreover, we determined Cu concentrations present in bovine oviductal fluid (OF), and the effects of Cu on fertilizing ability of bovine spermatozoa. Also, the presence of Cu transporter, SLC31A1 (also known as CTR1), in spermatozoa and in vitro matured oocyte were investigated. We found no differences in pregnancy rates among animals with adequate, marginal, and deficient Cu concentrations measured in plasma at the beginning of fixed-time artificial insemination (FTAI) protocol. Copper concentrations in OF were 38.3 ± 2.17 µg/dL (mean ± SEM) regardless of cupremia levels. The addition of 40 µg/dL Cu to IVF medium enhanced total and progressive motility, sperm viability, functional sperm membrane integrity (HOST), sperm-zona binding, and pronuclear formation. On the other hand, the presence of Cu in IVF medium did not modify acrosome integrity and cleavage rates after IVF, but impaired blastocyst rates. Cu transporter SLC31A1 was detected in bovine spermatozoa in the apical segment of acrosome, and in the oocyte matured in vitro. In conclusion, the results obtained in the present study determined that cupremia levels at the beginning of FTAI protocol did not influence the pregnancy rates at 60 d after insemination. The presence of CTR1 in bovine mature oocyte and spermatozoa, as well as the beneficial effect of Cu on sperm quality would suggest an important role of this mineral during the fertilization process.

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.

Sex chromosomal abnormalities associated with equine infertility: validation of a simple molecular screening tool in the Purebred Spanish Horse.

Chromosomal abnormalities in the sex chromosome pair (ECAX and ECAY) are widely associated with reproductive problems in horses. However, a large proportion of these abnormalities remains undiagnosed due to the lack of an affordable diagnostic tool that allows for avoiding karyotyping tests. Hereby, we developed an STR (single-tandem-repeat)-based molecular method to determine the presence of the main sex chromosomal abnormalities in horses in a fast, cheap and reliable way. The frequency of five ECAX-linked (LEX026, LEX003, TKY38, TKY270 and UCDEQ502) and two ECAY-linked (EcaYH12 and SRY) markers was characterized in 261 Purebred Spanish Horses to determine the efficiency of the methodology developed to be used as a chromosomal diagnostic tool. All the microsatellites analyzed were highly polymorphic, with a sizeable number of alleles (polymorphic information content > 0.5). Based on this variability, the methodology showed 100% sensitivity and 99.82% specificity to detect the most important sex chromosomal abnormalities reported in horses (chimerism, Turner's syndrome and sex reversal syndromes). The method was also validated with 100% efficiency in 10 individuals previously diagnosed as chromosomally aberrant. This STR screening panel is an efficient and reliable molecular-cytogenetic tool for the early detection of sex chromosomal abnormalities in equines that could be included in breeding programs to save money, effort and time of veterinary practitioners and breeders.

First case of sterility associated with sex chromosomal abnormalities in a jenny.

Chromosomal abnormalities are one of the main causes of genetic infertility in horses. Currently, their detection rate is rising due to the use of new diagnostic tools employing molecular markers linked to the sex chromosome pair. Despite genetic similarities, there are no previous reports of sterility associated with chromosomal abnormalities in the domestic donkey (Equus asinus). Hereby, we determined the presence of a chromosomal mosaicism in a female donkey with reproductive problems using molecular methodologies developed for horses. A two-and-a-half-year-old jenny characterized by morphological abnormalities of the reproductive tract was cytogenetically analysed using conventional and fluorescent techniques and a group of microsatellite markers (short tandem repeat, STR). At the same time, five ultrasound measures of the reproductive tract were taken and compared with eight contemporary jennies of the same breed. After slaughter, morphological examinations showed that the case study had a blind vaginal vestibule defining an empty pouch that covered the entrance of the cervical os. Histopathological studies demonstrated that this abnormal structure was compatible with a remnant hymen. Molecular markers, STR and fluorescent in situ hybridization determinations revealed that the animal was a 62, XX/61,X mosaic and, therefore, the first case of chromosomal abnormalities in the sex pair reported in donkeys.

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.

Effectiveness of a 95 SNP panel for the screening of breed label fraud in the Chinese meat market.

Breed assignment has proved to be useful to control meat trade and protect the value of special productions. Meat-related frauds have been detected in China; therefore, 95 SNPs selected from the ISAG core panel were evaluated to develop an automated and technologically updated tool to screen breed label fraud in the Chinese meat market. A total of 271 animals from four Chinese yellow cattle (CYC) populations, six Bos taurus breeds, two Bos indicus and one composite were used. The allocation test distinguished European, Japanese and Zebu breeds, and two Chinese genetic components. It correctly allocated Japanese Black, Zebu and British breeds in 100, 90 and 89% of samples, respectively. CYC evidenced the Zebu, Holstein and Limousin introgression. The test did not detect CYC components in any of the 25 samples from Argentinean butchers. The method could be useful to certify Angus, Hereford and Japanese Black meat, but a modification in the panel would be needed to differentiate other breeds.

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.

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.

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.

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.