Molecular tools and analytical approaches for the characterization of farm animal genetic diversity

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

Population diversity, admixture, and demographic trend of the Sumba Ongole cattle based on genomic data

OBJECTIVE

Sumba Ongole (SO) cattle are valuable breed due to their important role in the development of Indonesian cattle. Despite rapid advances in molecular technology, no genomic studies on SO cattle have been conducted to date. The aim of this study is to provide genomic profile related to the population diversity, admixture, and demographic trends of SO cattle.

METHODS

Genomic information was gathered from 79 SO cattle using the Illumina Bovine SNP50 v3 Beadchip, and for comparative purposes, additional genotypes from 209 cattle populations worldwide were included. The expected and observed heterozygosity, inbreeding coefficient, pairwise fixation indices between-population, and Nei's genetic distance were examined. Multidimensional scaling, admixture, and treemix analyses were used to investigate the population structure. Based on linkage disequilibrium and effective population size calculations, the demographic trend was observed.

RESULTS

The findings indicated that the genetic diversity of SO cattle was similar to that of other indicine breeds. SO cattle were genetically related to indicines but not to taurines or Bali cattle. The study further confirmed the close relationship between SO, Ongole, and Nellore cattle. Additionally, a small portion of the Ongole mixture were identified dominant in the SO population at the moment. The study also discovered that SO and Bali cattle (Bos javanicus) could have been ancestors in the development of Ongole Grade cattle, which corresponds to the documented history of Ongolization. Our finding indicate that SO cattle have maintained stability and possess unique traits separate from their ancestors.

CONCLUSION

In conclusion, the genetic diversity of the SO cattle has been conserved as a result of the growing significance of the present demographic trend. Consistent endeavors are necessary to uphold the fitness of the breed.

Exploring the Effects of Robertsonian Translocation 1/29 (Rob (1;29)) on Genetic Diversity in Minor Breeds of Spanish Berrenda Cattle via Genome-Wide Analysis

Most of the previous studies on the genetic variability in Spanish "Berrenda" breeds have been carried out using DNA microsatellites. The present work aimed to estimate the genetic diversity, population structure, and potential genetic differences among individuals of both Berrenda breeds and groups based on the presence of the Robertsonian chromosomal translocation, rob (1;29). A total of 373 samples from animals belonging to the two breeds, including 169 cases diagnosed as rob (1;29)-positive, were genotyped using an SNP50K chip. The genetic diversity at the breed level did not show significant differences, but it was significantly lower in those subpopulations containing the rob (1;29). Runs of homozygosity identified a region of homozygosity on chromosome 6, where the KIT (KIT proto-oncogene, receptor tyrosine kinase) gene, which determines the typical spotted coat pattern in both breeds, is located. The four subpopulations considered showed minor genetic differences. The regions of the genome that most determined the differences between the breeds were observed on chromosomes 4, 6, 18, and 22. The presence of this Robertsonian translocation did not result in sub-structuring within each of the breeds considered. To improve the reproductive performance of Berrenda breeds, it would be necessary to implement strategies considering the involvement of potential breeding stock carrying rob (1;29).

Structure and genetic variability of the Costa Rican Paso horse

The Costa Rican Paso Horse (CPC) is a breed developed in Costa Rica. The objectives were to estimate the genetic structure and evaluate the levels of genetic variability of the population. The genotypes of 14 microsatellites in 3654 records (2052 females and 1602 males) were analyzed. Expected (He) and observed (Ho) heterozygosity, polymorphic information content (CIP), fixation index (FIS), Shannon index, as well as Hardy-Weinberg disequilibrium (DHW) were evaluated. Kinship relationships (Rij) were estimated throughout the entire population. The effective population size (Ne) was calculated, alternating allele frequencies less than 0.05, 0.02 and 0.01. The Bayesian clustering study was carried out to infer how many lines are appropriate from the analysis of genotypes using multiple loci. The number of alleles per locus ranged from 7 to 17, with an average value of 9.6; nine loci presented DHW (P < 0.05); two loci presented negative FIS values, the same as Ho > He; the average of CIP, Ho and He was 0.254, 0.756 and 0.785, respectively. At the 12 loci where He > Ho, the differences ranged from 0.002 to 0.341 (0.036 on average). For Ne, the estimates were 201.9, 230.1, and 241.5. In the Rij, 54.86% of the estimates were in the interval of 0.01 to 77.7%. The number of lines that define the population corresponds to three, with an approximate composition of 33.1%, 32.4% and 34.5%, respectively. The CPC, as a subdivided population with DHW and a reduction in heterozygotes may be associated with possible Wahlund effects. Keywords: Wahlund effect, equines, genetic markers, synthetic breed, Hardy Weinberg.

The novel IMAGE001 genotyping array as a valuable alternative for genetic diversity screening in chicken: a demonstration in a local chicken breed in Belgium

Screening for genetic diversity in livestock species breeds is of utmost importance, especially for local, small populations that are at the risk of extinction. Luckily, recent developments in technology increase access to genotyping, also for numerically small breeds. One of these new technologies is the IMAGE001 single nucleotide polymorphism genotyping array that includes markers for 6 different species (cow, pig, sheep, chicken, horse and goat). For our current study, we studied the Turkey-headed Malines chicken, a local chicken breed in Belgium, for the first time. A total of 110 animals were genotyped, together with 29 samples from 4 supposedly related breeds. The genotypes were used to assess the genetic diversity of this local breed. Our analysis revealed an average inbreeding coefficient of 0.20 through runs of homozygosity analysis, and effective population size estimation based on linkage disequilibrium indicated a low genetic diversity (Ne = 34). Moreover, a principal component analysis and a genetic differentiation study (FST) were performed using these marker data to position the Turkey-headed Malines relative to the 4 other indigenous Belgian chicken breeds. Finally, we discussed the practical implications of the overlap between the IMAGE001 array and other existing chicken genotyping arrays. This study is the first use of the novel IMAGE001 array to evaluate a local chicken breed, and demonstrates it as a viable option for genomic characterization a breed. Moreover, with this research, we are able to provide a good basis for further evaluation of the Belgian chicken heritage.

Genome-wide mapping of signatures of selection using a high-density array identified candidate genes for growth traits and local adaptation in chickens

BACKGROUND

Availability of single nucleotide polymorphism (SNP) genotyping arrays and progress in statistical analyses have allowed the identification of genomic regions and genes under selection in chicken. In this study, SNP data from the 600 K Affymetrix chicken array were used to detect signatures of selection in 23 local Italian chicken populations. The populations were categorized into four groups for comparative analysis based on live weight (heavy vs light) and geographical area (Northern vs Southern Italy). Putative signatures of selection were investigated by combining three extended haplotype homozygosity (EHH) statistical approaches to quantify excess of haplotype homozygosity within (iHS) and between (Rsb and XP-EHH) groups. Presence of runs of homozygosity (ROH) islands was also analysed for each group.

RESULTS

After editing, 541 animals and 313,508 SNPs were available for statistical analyses. In total, 15 candidate genomic regions that are potentially under selection were detected among the four groups: eight within a group by iHS and seven by combining the results of Rsb and XP-EHH, which revealed divergent selection between the groups. The largest overlap between genomic regions identified to be under selection by the three approaches was on chicken chromosome 8. Twenty-one genomic regions were identified with the ROH approach but none of these overlapped with regions identified with the three EHH-derived statistics. Some of the identified regions under selection contained candidate genes with biological functions related to environmental stress, immune responses, and disease resistance, which indicate local adaptation of these chicken populations.

CONCLUSIONS

Compared to commercial lines, local populations are predominantly reared as backyard chickens, and thus, may have developed stronger resistance to environmental challenges. Our results indicate that selection can play an important role in shaping signatures of selection in local chicken populations and can be a starting point to identify gene mutations that could have a useful role with respect to climate change.

Genetic Diversity and Population Structure of Selected Ethiopian Indigenous Cattle Breeds Using Microsatellite Markers

Background

In Ethiopia, livestock contributes 45% of agricultural GDP. Despite the economic role played by the sector, there have been little efforts to genetically improve the indigenous cattle. Morphological characterization of selected Ethiopian indigenous cattle has been made for (Bonga, Jimma, and Kerayu) cattle types. But, the selected indigenous cattle were not characterized at molecular level (genetic diversity information). Hence, this work was initiated to detect and determine the genetic diversity and population structure of selected Ethiopian indigenous cattle ecotypes using microsatellite markers.

Results

Different alleles were identified (131) and thirty-three of these alleles were unique to specific ecotypes. All loci used were informative with PIC values ranging from 0.5 (TGLA126) to 0.84 (ETH10) with a mean of 0.70 per locus. The Shannon information index ranged from (I = 1.02) ILST006 to (I = 1.63) ETH10 with an average of 1.28 revealing there is genetic diversity. Moreover, analysis of molecular variance (AMOVA) revealed 84% genetic variation within a population and 13% variation among populations. The value of F-statistics (Fst) (0.129 = 13%) indicated that there was moderate genetic differentiation among ecotypes. The (UPGMA) revealed, Bonga and Jimma clustered together while Kerayu cattle were relatively distinct, Principal coordinates analysis (PCOA) and structure analysis grouped the individual into different clusters confirming the presence of ecotype admixture due to geographical origins and uncontrolled mating.

Conclusion

In general, this study has successfully characterized the genetic diversity and population structure of Bonga, Jimma, and Kerayu cattle ecotypes using high polymorphic/informative microsatellite markers. According to this study, Kerayu cattle have high AR and PA when compared to Bonga and Jimma cattle populations. So, the Kerayu population is more diverse than others and it is the hotspot for genetic diversity study. The generated information is very relevant for breeder and genetic conservation.

Genome-wide analysis for the melatonin trait associated genes and SNPs in dairy goat (Capra hircus) as the molecular breeding markers

Previous studies have reported that the endogenous melatonin level is positively associated with the quality and yield of milk of cows. In the current study, a total of 34,921 SNPs involving 1,177 genes were identified in dairy goats by using the whole genome resequencing bulked segregant analysis (BSA) analysis. These SNPs have been used to match the melatonin levels of the dairy goats. Among them, 3 SNPs has been identified to significantly correlate with melatonin levels. These 3 SNPs include CC genotype 147316, GG genotype 147379 and CC genotype 1389193 which all locate in the exon regions of ASMT and MT2 genes. Dairy goats with these SNPs have approximately 5-fold-higher melatonin levels in milk and serum than the average melatonin level detected in the current goat population. If the melatonin level impacts the milk production in goats as in cows, the results strongly suggest that these 3 SNPs can serve as the molecular markers to select the goats having the improved milk quality and yield. This is a goal of our future study.

Inbreeding and pedigree analysis of the European red dairy cattle

BACKGROUND

Red dairy cattle breeds have an important role in the European dairy sector because of their functional characteristics and good health. Extensive pedigree information is available for these breeds and provides a unique opportunity to examine their population structure, such as effective population size, depth of the pedigree, and effective number of founders and ancestors, and inbreeding levels. Animals with the highest genetic contributions were identified. Pedigree data included 9,073,403 animals that were born between 1900 and 2019 from Denmark, Finland, Germany, Latvia, Lithuania, the Netherlands, Norway, Poland, and Sweden, and covered 32 breeds. The numerically largest breeds were Red Dairy Cattle and Meuse-Rhine-Yssel.

RESULTS

The deepest average complete generation equivalent (9.39) was found for Red Dairy Cattle in 2017. Mean pedigree completeness ranged from 0.6 for Finncattle to 7.51 for Red Dairy Cattle. An effective population size of 166 animals was estimated for the total pedigree and ranged from 35 (Rotes Höhenvieh) to 226 (Red Dairy Cattle). Average generation intervals were between 5 and 7 years. The mean inbreeding coefficient for animals born between 1960 and 2018 was 1.5%, with the highest inbreeding coefficients observed for Traditional Angler (4.2%) and Rotes Höhenvieh (4.1%). The most influential animal was a Dutch Meuse-Rhine-Yssel bull born in 1960. The mean inbreeding level for animals born between 2016 and 2018 was 2% and highest for the Meuse-Rhine-Yssel (4.64%) and Rotes Hohenvieh breeds (3.80%).

CONCLUSIONS

We provide the first detailed analysis of the genetic diversity and inbreeding levels of the European red dairy cattle breeds. Rotes Höhenvieh and Traditional Angler have high inbreeding levels and are either close to or below the minimal recommended effective population size, thus it is necessary to implement tools to monitor the selection process in order to control inbreeding in these breeds. Red Dairy Cattle, Vorderwälder, Swedish Polled and Hinterwälder hold more genetic diversity. Regarding the Meuse-Rhine-Yssel breed, given its decreased population size, increased inbreeding and low effective population size, we recommend implementation of a breeding program to prevent further loss in its genetic diversity.

Genetic Architecture and Signatures of Selection in the Caqueteño Creole (Colombian Native Cattle)

Evolutionary mechanisms have shaped the genomic architecture of Colombian Creole cattle breeds. The mating and selection processes have impacted several traits, promoting differences within and between populations. Studies of population structure and selection signatures in Colombian Creole breeds are scarce, and need more attention to better understand genetic differentiation, gene flow, and genetic distance. This study aimed to analyze the population structure and identify selection imprints in the Criollo Caqueteño (CAQ) population. It used 127 CAQ animals genotyped with Chip HD 777,000 SNPs. The population structure analyses used discriminant principal component analysis (DAPC), integrated haplotype scoring (iHS), and index-fixing (Fst) methodologies to detect selection signals. We can highlight SNP regions on the genes TMPRSS15, PGAM2, and EGFR, identified by the Fst method. Additionally, the iHS regions for cluster 1 identified candidate genes on BTA 3 (CMPK1 and FOXD2), BTA 11 (RCAN1), and BTA 22 (ARPP21). In group 2, we can highlight the genes on BTA 4 (SLC13A4, BRAF), BTA 9 (ULBP), BTA 14 (CSMD3) and BTA 19 (KRTAP9-2). These candidate genes have been associated with fertility traits, precocity, growth, and environmental and disease resistance, indicating a genetic potential in CAQ animals. All this promotes a better understanding of the diversity and genetic structure in the CAQ population. Based on that, our study can significantly assist the sustainable development and conservation of the breed in the Colombian Amazon.

Genomic Population Structure of the Main Historical Genetic Lines of Spanish Merino Sheep

Simple Summary

Historical documentation shows that the Spanish Merino sheep was selected over many centuries due to the quality of wool, following which it was used to originate all other Merino breeds around the world, mainly by crossbreeding with local breeds. Today, the historical genetic lines that originated the Spanish Merino are still preserved in several closed herds in which they have been bred for nearly 200 years, maintaining their original genetic purity. Our study demonstrates, using a genomic approach, the exceptional genetic richness and variability of these lines, which are clearly differentiated from modern Merino breeds, and must therefore be protected to safeguard the large genetic pool they represent.

Abstract

According to historiographical documentation, the Romans first began to select Merino sheep in the Iberian Peninsula during the first century, with the aim of obtaining a breed appreciated for the quality of its wool. This process continued locally during the Middle Ages, when Spanish sheep were protected, and their export to foreign countries was banned. It was during the 16th century when individual Merino sheep were allowed to spread around the world to be used to improve the wool quality of local breeds. However, the wool crisis of the 1960s shifted the selection criteria of the Merino breed towards meat production at the expenses of wool. Consequently, individuals that display the genetic and phenotypic characteristics of those sheep originally bred in the kingdom of Spain in the Middle Ages are extremely difficult to find in commercial herds. In this study, we characterized the genetic basis of 403 individuals from the main historical Spanish Merino genetic lines (Granda, Hidalgo, Lopez-Montenegro, Maeso, Donoso and Egea), which were bred in isolation over the last 200 years, using a genomic approach based on genotyping data from the Axiom™ Ovine 50K SNP Genotyping Array. Our analysis included measuring population structure, genomic differentiation indexes, runs of homozygosity (ROH) patterns, and an analysis of molecular variance (AMOVA). The results showed large genetic differences between the historical lines, even though they belong to the same breed. In addition, ROH analysis showed differences due to increased inbreeding among the ancient generations compared with the modern Merino lines, confirming the breed’s ancestral and closed origin. However, our results also showed a high variability and richness within the Spanish historical Merino lines from a genetic viewpoint. This fact, together with their great ability to produce high-quality wool, suggests that ancestral Merino lines from Spain should be considered a valuable genetic population to be maintained as a resource for the improvement of wool-producing sheep breeds all around the world.

High Genetic Diversity but Absence of Population Structure in Local Chickens of Sri Lanka Inferred by Microsatellite Markers

Local chicken populations belonging to five villages in two geographically separated provinces of Sri Lanka were analyzed using 20 microsatellite markers to determine the genetic diversity of local chickens. Population genetic parameters were estimated separately for five populations based on geographic locations and for eight populations based on phenotypes, such as naked neck, long legged, crested or crown, frizzle feathered, Giriraj, commercial layer, crossbreds, and non-descript chicken. The analysis revealed that there was a high genetic diversity among local chickens with high number of unique alleles, mean number of alleles per locus (MNA), and total number of alleles per locus per population. A total of 185 microsatellite alleles were detected in 192 samples, indicating a high allelic diversity. The MNA ranged from 8.10 (non-descript village chicken) to 3.50 (Giriraj) among phenotypes and from 7.30 (Tabbowa) to 6.50 (Labunoruwa) among village populations. In phenotypic groups, positive inbreeding coefficient (F IS) values indicated the existence of population substructure with evidence of inbreeding. In commercial layers, a high expected heterozygosity He = 0.640 ± 0.042) and a negative F IS were observed. The positive F IS and high He estimates observed in village populations were due to the heterogeneity of samples, owing to free mating facilitated by communal feeding patterns. Highly admixed nature of phenotypes was explained as a result of rearing many phenotypes by households (58%) and interactions of chickens among neighboring households (53%). A weak substructure was evident due to the mating system, which disregarded the phenotypes. Based on genetic distances, crown chickens had the highest distance to other phenotypes, while the highest similarity was observed between non-descript village chickens and naked neck birds. The finding confirms the genetic wealth conserved within the populations as a result of the breeding system commonly practiced by chicken owners. Thus, the existing local chicken populations should be considered as a harbor of gene pool, which can be readily utilized in developing locally adapted and improved chicken breeds in the future.

Exploring genetic diversity and population structure of Punjab goat breeds using Illumina 50 K SNP bead chip

Pakistan has 35 goat breeds. Moreover, the province of Punjab has highest goat population constituting 37% of country's total population with seven goat breeds including Beetal, Daira Deen Panah, Nachi, Barbari, Teddi, Pahari, and Pothwari. The diversity study of breeds warrants the documentation of breeds particularly using genome wide panel of markers, i.e., SNP chip. The objective of the current study was to fill this gap of information. Therefore, in current study we collected total of 879 unrelated goat blood samples along with data on body weight measurements; genomic DNA was extracted, and genotyping was carried out using 50 K SNP bead chip. Quality control measures were performed in Plink 1.07. Genetic diversity was observed among studied populations using heterozygosity and pairwise F_ST estimates, principal component analysis, admixture analysis in Plink software with visualization in Clumpak, and constructing phylogenetic tree in Mega 7 software. Moderate to high level of heterozygosity was observed among the studied populations. Coefficient of inbreeding varied from 0.0186 ± 0.0327 in Pahari to 0.183 ± 0.0715 in Barbari. Barbari and Daira Deen Panah had quite higher level of inbreeding coefficient as compared to all other breeds with value of 0.183 ± 0.0715 and 0.1378 ± 0.0741, respectively. PCA identified three steps of subdividing the seven goat breeds at various levels of K. All the seven breeds made independent clusters at various levels of PCA. Admixture analysis revealed the distinctness of Teddi and Barbari breeds. Genetic sub-structuring was observed in the admixture patterns of Beetal breed. Moreover, high level of genetic admixture was observed in Nachi, Pahari, Pothwari, and Daira Deen Panah breeds. Admixture results were further interpreted by calculating pairwise F_ST values. Our results provided first insights about genetic diversity of Pakistani goat breeds based on genomic data. To conclude, the enriched goat breed diversity in Pakistan could provide valuable genetic reservoir for national breeding schemes.

Assessing the genetic background and genomic relatedness of red cattle populations originating from Northern Europe

Background

Local cattle breeds need special attention, as they are valuable reservoirs of genetic diversity. Appropriate breeding decisions and adequate genomic management of numerically smaller populations are required for their conservation. At this point, the analysis of dense genome-wide marker arrays provides encompassing insights into the genomic constitution of livestock populations. We have analyzed the genetic characterization of ten cattle breeds originating from Germany, The Netherlands and Denmark belonging to the group of red dairy breeds in Northern Europe. The results are intended to provide initial evidence on whether joint genomic breeding strategies of these populations will be successful.

Results

Traditional Danish Red and Groningen White-Headed were the most genetically differentiated breeds and their populations showed the highest levels of inbreeding. In contrast, close genetic relationships and shared ancestry were observed for the populations of German Red and White Dual-Purpose, Dutch Meuse-Rhine-Yssel, and Dutch Deep Red breeds, reflecting their common histories. A considerable amount of gene flow from Red Holstein to German Angler and to German Red and White Dual-Purpose was revealed, which is consistent with frequent crossbreeding to improve productivity of these local breeds. In Red Holstein, marked genomic signatures of selection were reported on chromosome 18, suggesting directed selection for important breeding goal traits. Furthermore, tests for signatures of selection between Red Holstein, Red and White Dual-Purpose, and Meuse-Rhine-Yssel uncovered signals for all investigated pairs of populations. The corresponding genomic regions, which were putatively under different selection pressures, harboured various genes which are associated with traits such as milk and beef production, mastitis and female fertility.

Conclusions

This study provides comprehensive knowledge on the genetic constitution and genomic connectedness of divergent red cattle populations in Northern Europe. The results will help to design and optimize breeding strategies. A joint genomic evaluation including some of the breeds studied here seems feasible.

Current genetic conservation of Chinese indigenous horses revealed with Y-chromosomal and mitochondrial DNA polymorphisms

To investigate the genetic diversity of Chinese indigenous horses and determine the genetic status of extant horse breeds, novel Y chromosomal microsatellite markers and known Y chromosomal SNPs and mtDNA loop sequences, were employed to study the genetic diversity levels of 13 Chinese indigenous horse populations and four introduced breeds. Sixteen Y-chromosomal microsatellite markers, including seven newly identified loci, were used in the genotyping. The results showed that 4 out of the 16 loci were highly polymorphic in Chinese indigenous horse populations, in which the polymorphisms of 3 loci, ECAYP12, ECAYP13, and ECAYCAU3, were first reported in the present study. The polymorphic Y chromosomal microsatellite markers result in 19 haplotypes in the studied horses and formed 24 paternal lines when merged with the 14 Y chromosomal SNPs reported previously. The haplotypes CHT18 and SS24 harboring AMELY gene mutation were the ancestral haplotypes, and other haplotypes were derived from them by one or more mutation steps. The horse populations in mountainous and remote areas of southwestern China have the most ancient paternal lines, which suggests that ancient paternal lines preserved in local populations attributed to less human interventions. Our results also showed that the northern local breeds had higher mtDNA diversity than the southern ones in China. The frequency of haplogroup B, F, and G of mtDNA in Chinese indigenous horses has declined in recent years, and some breeds are in endangered status mainly due to small population sizes. Urgent actions should be taken to conserve the genetic diversity of the indigenous horse populations, especially the rare paternal lines. Our findings help to elucidate the genetic diversity and evolutionary history of Chinese domestic horses, which will facilitate the conservation of the indigenous horses in the future.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Assessment of genetic diversity and conservation priorities in some Turkish indigenous Hair goat populations by microsatellite loci

Genetic diversity in livestock breeds is required for breeding studies, response to environmental changes and sustainable production. The aim of this study was to evaluate the genetic diversity in Hair goats reared in 9 districts of Antalya province and to determine the populations that have the highest contribution to the total genetic diversity. For this purpose, 180 samples from 9 districts (Korkuteli-KRK, Elmalý-ELM, Kaþ-KAS, Demre-DMR, Manavgat- MNG, Gündoðmuþ-GND, Ýbradý-IBR, Akseki-AKS and Gazipaþa-GZP) of Antalya province were genotyped by 20 microsatellite loci. The mean number of alleles per locus for each population ranged from 8.45 (GND) to 9.25 (MNG), while mean number of effective allele varied between 5.40 (GND) and 6.22 (MNG). The lowest average observed heterozygosity was in the ELM populations (0.71) while the highest Ho value detected in KAS populations (0.78). Mean expected heterozygosity values varied from 0.80 (GND) to 0.84 (DMR, MNG). Mean PIC values ranged from 0.77 (GND, AKS) to 0.80 (DMR, MNG) in populations. Inbreeding coefficients were detected between 0.05 (KAS) and 0.13 (ELM) in district populations. According to two different methods, the highest contribution to the total genetic diversity comes from KAS (-0.244) and AKS populations (0.482). In conclusion, high genetic diversity and low level of inbreeding were determined in Turkish indigenous Hair goats. Hair goats have great potential for breeding studies and for adaptation to the environmental conditions that will possibly change in the future. Especially, genetic variation in KAS and AKS populations should be conserved.

A Genetic Window on Sardinian Native Horse Breeds through Uniparental Molecular Systems

Simple Summary

The horse is a mammalian species showing a high variation among maternal lineages but a limited variability in the paternal inheritance. The female phylogenetic history is commonly investigated by analyzing the maternally transmitted mitochondrial DNA (mtDNA), while the male perspective is provided by the paternally inherited portion of the Y chromosome (NRY). Here we explored the variation of both non-recombining genetic systems in three horse breeds present in Sardinia: Giara, Sarcidano, and Sardinian Anglo-Arab. The analysis of 34 stallions revealed three differentially distributed NRY types: (i) the first and most ancestral one is typical of Sarcidano; (ii) the second is well represented in Giara and seems to derive from Neapolitan/Oriental stallions; (iii) the third confirms the Thoroughbred influence in the Sardinian Anglo-Arab breed. By extending the analysis to 178 mtDNAs, we observed a common maternal origin for Giara and Sarcidano. Contrarily, the outlier behavior of the Sardinian Anglo-Arab is due to its higher mitochondrial variability, testifying for multiple maternal lineages in its current population. Our preliminary findings highlight the importance of a parallel molecular screening of NRYs and mtDNAs to reconstruct both paternal and maternal phylogenetic histories and to fully evaluate the extent of autochthonous genetic resources in the island.

Abstract

Sardinia, an island located to the west of Italy in the Mediterranean Sea, boasts three native horse breeds: Giara, Sarcidano, and Sardinian Anglo-Arab. Here, we have investigated for the first time three loci of the non-recombining region of the Y chromosome (NRY) in 34 stallions from these breeds and performed a phylogenetic analysis of the maternal relationships among 178 previously published mitochondrial control regions. We found that the current NRY diversity of Sardinian horse breeds is linked to three haplotypes (HT), all identified within Sarcidano. Each breed showed a typical HT: HT1 (ancestral) was the most represented in Sarcidano, HT2 (Neapolitan/Oriental wave) in Giara, and HT3 (Thoroughbred wave) in Sardinian Anglo-Arab. The specificity of each haplotype suggests the influence of independent breeding strategies and the effect of genetic drift in each Sardinian population. The female counterpart, extended to 178 horses, showed a low genetic variability and a common maternal origin for Giara and Sarcidano. The higher variability of the Sardinian Anglo-Arab indicates multiple mare lineages in its current population. Further genetic analyses will be crucial to understand the paternal history of male horses, preserve the endangered mares’ and stallions’ lineages, and improve the enhancement of autochthonous genetic resources on this island.

Estimating genetic diversity and population structure of 22 chicken breeds in Asia using microsatellite markers

Objective

Estimating the genetic diversity and structures, both within and among chicken breeds, is critical for the identification and conservation of valuable genetic resources. In chickens, microsatellite (MS) marker polymorphisms have previously been widely used to evaluate these distinctions. Our objective was to analyze the genetic diversity and relationships among 22 chicken breeds in Asia based on allelic frequencies.

Methods

We used 469 genomic DNA samples from 22 chicken breeds from eight Asian countries (South Korea, KNG, KNB, KNR, KNW, KNY, KNO; Laos, LYO, LCH, LBB, LOU; Indonesia, INK, INS, ING; Vietnam, VTN, VNH; Mongolia, MGN; Kyrgyzstan, KGPS; Nepal, NPS; Sri Lanka, SBC) and three imported breeds (RIR, Rhode Island Red; WLG, White Leghorn; CON, Cornish). Their genetic diversity and phylogenetic relationships were analyzed using 20 MS markers.

Results

In total, 193 alleles were observed across all 20 MS markers, and the number of alleles ranged from 3 (MCW0103) to 20 (LEI0192) with a mean of 9.7 overall. The NPS breed had the highest expected heterozygosity (H_exp, 0.718±0.027) and polymorphism information content (PIC, 0.663±0.030). Additionally, the observed heterozygosity (H_obs) was highest in LCH (0.690±0.039), whereas WLG showed the lowest H_exp (0.372±0.055), H_obs (0.384±0.019), and PIC (0.325±0.049). Nei’s DA genetic distance was the closest between VTN and VNH (0.086), and farthest between KNG and MGN (0.503). Principal coordinate analysis showed similar results to the phylogenetic analysis, and three axes explained 56.2% of the variance (axis 1, 19.17%; 2, 18.92%; 3, 18.11%). STRUCTURE analysis revealed that the 22 chicken breeds should be divided into 20 clusters, based on the highest ΔK value (46.92).

Conclusion

This study provides a basis for future genetic variation studies and the development of conservation strategies for 22 chicken breeds in Asia.

Population genetic structure and maternal lineage of South African crossbred Nguni cattle using the cytochrome b gene in mtDNA

The Nguni cattle breed predominates South Africa and is endowed with traits favourable against environmental stressors such as heat stress and resistance to diseases. Interventions to improve production have led to the erosion of the genetic integrity of local breeds and the introduction of exotic breeds has proved ineffective as they fail to perform well due to different climatic conditions and production systems. In this study, the genetic structure and genetic lineage of Nguni crossbreds from 6 populations were assessed using the mitochondrial cytochrome b gene. Twelve polymorphic sites were detected resulting in 11 haplotypes with haplotype and nucleotide diversities of 0.550 ± 0.135 and 0.0019 ± 0.0011, respectively. Only 2 of the 6 populations displayed recent population expansion events, whereas the majority adhered to neutral evolution. The basal haplotype contained approximately 60% of the studied populations and there were four unique haplotypes that were revealed. A possible Nguni descript haplotype was uncovered, and this haplotype was found in all populations but was however devoid of individuals from around the world. The genetic structure of the populations was rather low (average pairwise F_ST = 0.066 and Slatkins F_ST = 0.094), and approximately 96% of the total genetic variation was accounted for by differences within populations. Phylogenetic analyses supported the clustering of all the samples within the Bos taurus clade and no Bos indicus haplotype was detected. Furthermore, no intermediate haplotype of taurine and indicine was detected. Overall, the maternal lineage of the crossbreds points to a taurine origin and the low genetic diversity depicts the retention of the Nguni genetic pool and possibly its superior adaptive traits.

Genetic structure and admixture in sheep from terminal breeds in the United States

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright’s F_ST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright’s F_ST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture, and based on cross‐validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.

Asian water buffalo: domestication, history and genetics

The domestic Asian water buffalo (Bubalus bubalis) is found on all five continents, with a global population of some 202 million. The livelihoods of more people depend on this species than on any other domestic animal. The two distinct types (river and swamp) descended from different wild Asian water buffalo (Bubalus arnee) populations that diverged some 900 kyr BP and then evolved in separate geographical regions. After domestication in the western region of the Indian subcontinent (ca. 6300 years BP), the river buffalo spread west as far as Egypt, the Balkans and Italy. Conversely, after domestication in the China/Indochina border region ca. 3000-7000 years BP, swamp buffaloes dispersed through south-east Asia and China as far as the Yangtze River valley. Molecular and morphological evidence indicates that swamp buffalo populations have strong geographic genetic differentiation and a lack of gene flow, but strong phenotypic uniformity. In contrast, river buffalo populations show a weaker phylogeographic structure, but higher phenotypic diversity (i.e. many breeds). The recent availability of a high-quality reference genome and of a medium-density marker panel for genotyping has triggered a number of genome-wide investigations on diversity, evolutionary history, production traits and functional elements. The growing molecular knowledge combined with breeding programmes should pave the way to improvements in production, environmental adaptation and disease resistance in water buffalo populations worldwide.

Population structure of indigenous southern African goats based on the Illumina Goat50K SNP panel

In this study, the genetic structure of indigenous Tswana and Swazi goats using the Illumina Goat50K SNP array was investigated. Two South African commercial goat breeds were included to investigate admixture with the indigenous populations in southern Africa. A total of 144 DNA samples including Boer goats (n = 24), Kalahari Red (n = 24), Swazi (n = 48), and Tswana goats (n = 48) were genotyped. Statistical analysis was performed using PLINK version 1.07. Genetic diversity, measured as expected heterozygosity, was estimated at 0.390, 0.398, 0.413, and 0.387 for Boer, Kalahari Red, Tswana, and Swazi goats, respectively. The individual inbreeding coefficient varied from 0.019 ± 0.05 to 0.011 ± 0.06 for the Tswana and Swazi goats, respectively. The Principal component analysis clustered the populations according to geographical origin and breed type. Linkage disequilibrium (LD) for shorter intervals (0-10 kb) ranged from 0.44 to 0.56 and commercial breeds had higher values. Effective population sizes decreased with generations and at the 13th generation ranged between 87 for Boer to 266 for Tswana goats. The Tswana population exhibited the highest level of genetic variation and effective population size, which holds potential for improved production in marginal regions. A national strategy is required to maintain genetic diversity in communal goat production systems through well-structured breeding and conservation programs.

Sustainable Utilization of Indigenous Goats in Southern Africa

Goats have a key role in ensuring food security and economic livelihood to smallholder farmers in rural areas. Women play a vital role in goat rearing, promoting economic autonomy within households. Indigenous goats dominate and are of high significance due to their adaptive traits that are relevant for climate change and low maintenance. However, lack of emphasis on farmer-centered technology development and proper breed characterization remains a hitch to sustainable utilization and breed development of indigenous goats. This can be over come through proper linkage between market and production, workable regional and national agricultural policies, community breeding programs, collaborative research work within the region, and consistent government support.

Runs of Homozygosity and NetView analyses provide new insight into the genome-wide diversity and admixture of three German cattle breeds

Angler (RVA) and Red-and-White dual-purpose (RDN) cattle were in the past decades crossed with influential Red Holstein (RH) sires. However, genome-wide diversity studies in these breeds are lacking. The objective of the present study was to elucidate the genome-wide diversity and population structure of the three German cattle breeds. Using 40,851 single nucleotide polymorphism markers scored in 337 individuals, runs of homozygosity (ROH) were analysed in each breed. Clustering and a high-resolution network visualisation analyses were performed on an extended dataset that included 11 additional (outgroup) breeds. Genetic diversity levels were high with observed heterozygosity above 0.35 in all three breeds. Only RVA had a recent past effective population size (Ne) estimate above 100 at 5 generations ago. ROH length distribution followed a similar pattern across breeds and the majority of ROH were found in the length class of >5 to 10 Mb. Estimates of average inbreeding calculated from ROH (FROH) were 0.021 (RVA), 0.045 (RDN) and 0.053 (RH). Moderate to high positive correlations were found between FROH and pedigree inbreeding (FPED) and between FROH and inbreeding derived from the excess of homozygosity (FHOM), while the intercept of the regression of FROH on FPED was above zero. The population structure analysis showed strong evidence of admixture between RVA and RH. Introgression of RDN with RH genes was minimally detected and for the first time, the study uncovered Norwegian Red Cattle ancestry in RVA. Highly heterogeneous genetic background was found for RVA and RH and as expected, the breeds of the extended dataset effectively differentiated mostly based on geographical origin, validating our findings. The results of this study confirm the impact of RH sires on RVA and RDN populations. Furthermore, a close monitoring is suggested to curb further reduction of Ne in the breeds.

Signatures of selection and environmental adaptation across the goat genome post-domestication

BACKGROUND

Since goat was domesticated 10,000 years ago, many factors have contributed to the differentiation of goat breeds and these are classified mainly into two types: (i) adaptation to different breeding systems and/or purposes and (ii) adaptation to different environments. As a result, approximately 600 goat breeds have developed worldwide; they differ considerably from one another in terms of phenotypic characteristics and are adapted to a wide range of climatic conditions. In this work, we analyzed the AdaptMap goat dataset, which is composed of data from more than 3000 animals collected worldwide and genotyped with the CaprineSNP50 BeadChip. These animals were partitioned into groups based on geographical area, production uses, available records on solid coat color and environmental variables including the sampling geographical coordinates, to investigate the role of natural and/or artificial selection in shaping the genome of goat breeds.

RESULTS

Several signatures of selection on different chromosomal regions were detected across the different breeds, sub-geographical clusters, phenotypic and climatic groups. These regions contain genes that are involved in important biological processes, such as milk-, meat- or fiber-related production, coat color, glucose pathway, oxidative stress response, size, and circadian clock differences. Our results confirm previous findings in other species on adaptation to extreme environments and human purposes and provide new genes that could explain some of the differences between goat breeds according to their geographical distribution and adaptation to different environments.

CONCLUSIONS

These analyses of signatures of selection provide a comprehensive first picture of the global domestication process and adaptation of goat breeds and highlight possible genes that may have contributed to the differentiation of this species worldwide.

Next-generation metrics for monitoring genetic erosion within populations of conservation concern

Genetic erosion is a major threat to biodiversity because it can reduce fitness and ultimately contribute to the extinction of populations. Here, we explore the use of quantitative metrics to detect and monitor genetic erosion. Monitoring systems should not only characterize the mechanisms and drivers of genetic erosion (inbreeding, genetic drift, demographic instability, population fragmentation, introgressive hybridization, selection) but also its consequences (inbreeding and outbreeding depression, emergence of large-effect detrimental alleles, maladaptation and loss of adaptability). Technological advances in genomics now allow the production of data the can be measured by new metrics with improved precision, increased efficiency and the potential to discriminate between neutral diversity (shaped mainly by population size and gene flow) and functional/adaptive diversity (shaped mainly by selection), allowing the assessment of management-relevant genetic markers. The requirements of such studies in terms of sample size and marker density largely depend on the kind of population monitored, the questions to be answered and the metrics employed. We discuss prospects for the integration of this new information and metrics into conservation monitoring programmes.

Whole-genome SNP analysis elucidates the genetic structure of Russian cattle and its relationship with Eurasian taurine breeds

BACKGROUND

The origin of native and locally developed Russian cattle breeds is linked to the historical, social, cultural, and climatic features of the diverse geographical regions of Russia. In the present study, we investigated the population structure of nine Russian cattle breeds and their relations to the cattle breeds from around the world to elucidate their origin. Genotyping of single nucleotide polymorphisms (SNPs) in Bestuzhev (n = 26), Russian Black-and-White (n = 21), Kalmyk (n = 14), Kholmogor (n = 25), Kostromsky (n = 20), Red Gorbatov (n = 23), Suksun (n = 20), Yakut (n = 25), and Yaroslavl cattle breeds (n = 21) was done using the Bovine SNP50 BeadChip. SNP profiles from an additional 70 breeds were included in the analysis as references.

RESULTS

The observed heterozygosity levels were quite similar in eight of the nine studied breeds (H_O = 0.337-0.363) except for Yakut (Ho = 0.279). The inbreeding coefficients F_IS ranged from -0.028 for Kalmyk to 0.036 for Russian Black-and-White and were comparable to those of the European breeds. The nine studied Russian breeds exhibited taurine ancestry along the C1 axis of the multidimensional scaling (MDS)-plot, but Yakut was clearly separated from the European taurine breeds on the C2 axis. Neighbor-Net and admixture analyses, discriminated three groups among the studied Russian breeds. Yakut and Kalmyk were assigned to a separate group because of their Turano-Mongolian origin. Russian Black-and-White, Kostromsky and Suksun showed transboundary European ancestry, which originated from the Holstein, Brown Swiss, and Danish Red breeds, respectively. The lowest level of introgression of transboundary breeds was recorded for the Kholmogor, Yaroslavl, Red Gorbatov and Bestuzhev breeds, which can be considered as an authentic genetic resource.

CONCLUSIONS

Whole-genome SNP analysis revealed that Russian native and locally developed breeds have conserved authentic genetic patterns in spite of the considerable influence of Eurasian taurine cattle. In this paper, we provide fundamental genomic information that will contribute to the development of more accurate breed conservation programs and genetic improvement strategies.

Genetic structure in Red Junglefowl (Gallus gallus) populations: Strong spatial patterns in the wild ancestors of domestic chickens in a core distribution range

Red Junglefowl (Gallus gallus) are among the few remaining ancestors of an extant domesticated livestock species, the domestic chicken, that still occur in the wild. Little is known about genetic diversity, population structure, and demography of wild Red Junglefowl in their natural habitats. Extinction threats from habitat loss or genetic alteration from domestic introgression exacerbate further the conservation status of this progenitor species. In a previous study, we reported extraordinary adaptive genetic variation in the MHC B-locus in wild Red Junglefowl and no evidence of allelic introgression between wild and domestic chickens was observed. In this study, we characterized spatial genetic variation and population structure in naturally occurring populations of Red Junglefowl in their core distribution range in South Central Vietnam. A sample of 212 Red Junglefowl was obtained from geographically and ecologically diverse habitats across an area of 250 × 350 km. We used amplified fragment-length polymorphism markers obtained from 431 loci to determine whether genetic diversity and population structure varies. We found that Red Junglefowl are widely distributed but form small and isolated populations. Strong spatial genetic patterns occur at both local and regional scales. At local scale, population stratification can be identified to approximately 5 km. At regional scale, we identified distinct populations of Red Junglefowl in the southern lowlands, northern highlands, and eastern coastal portions of the study area. Both local and long-distance genetic patterns observed in wild Red Junglefowl may reflect the species' ground-dwelling and territorial characteristics, including dispersal barriers imposed by the Annamite Mountain Range. Spatially explicit analyses with neutral genetic markers can be highly informative and here elevates the conservation profile of the wild ancestors of domesticated chickens.

Conservation status and historical relatedness of Italian cattle breeds

Background

In the last 50 years, the diversity of cattle breeds has experienced a severe contraction. However, in spite of the growing diffusion of cosmopolite specialized breeds, several local cattle breeds are still farmed in Italy. Genetic characterization of breeds represents an essential step to guide decisions in the management of farm animal genetic resources. The aim of this work was to provide a high-resolution representation of the genome-wide diversity and population structure of Italian local cattle breeds using a medium-density single nucleotide polymorphism (SNP) array.

Results

After quality control filtering, the dataset included 31,013 SNPs for 800 samples from 32 breeds. Our results on the genetic diversity of these breeds agree largely with their recorded history. We observed a low level of genetic diversity, which together with the small size of the effective populations, confirmed that several breeds are threatened with extinction. According to the analysis of runs of homozygosity, evidence of recent inbreeding was strong in some local breeds, such as Garfagnina, Mucca Pisana and Pontremolese. Patterns of genetic differentiation, shared ancestry, admixture events, and the phylogenetic tree, all suggest the presence of gene flow, in particular among breeds that originate from the same geographical area, such as the Sicilian breeds. In spite of the complex admixture events that most Italian cattle breeds have experienced, they have preserved distinctive characteristics and can be clearly discriminated, which is probably due to differences in genetic origin, environment, genetic isolation and inbreeding.

Conclusions

This study is the first exhaustive genome-wide analysis of the diversity of Italian cattle breeds. The results are of significant importance because they will help design and implement conservation strategies. Indeed, efforts to maintain genetic diversity in these breeds are needed. Improvement of systems to record and monitor inbreeding in these breeds may contribute to their in situ conservation and, in view of this, the availability of genomic data is a fundamental resource.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0406-x) contains supplementary material, which is available to authorized users.

Population structure and genetic diversity of 25 Russian sheep breeds based on whole-genome genotyping

Background

Russia has a diverse variety of native and locally developed sheep breeds with coarse, fine, and semi-fine wool, which inhabit different climate zones and landscapes that range from hot deserts to harsh northern areas. To date, no genome-wide information has been used to investigate the history and genetic characteristics of the extant local Russian sheep populations. To infer the population structure and genome-wide diversity of Russian sheep, 25 local breeds were genotyped with the OvineSNP50 BeadChip. Furthermore, to evaluate admixture contributions from foreign breeds in Russian sheep, a set of 58 worldwide breeds from publicly available genotypes was added to our data.

Results

We recorded similar observed heterozygosity (0.354–0.395) and allelic richness (1.890–1.955) levels across the analyzed breeds and they are comparable with those observed in the worldwide breeds. Recent effective population sizes estimated from linkage disequilibrium five generations ago ranged from 65 to 543. Multi-dimensional scaling, admixture, and neighbor-net analyses consistently identified a two-step subdivision of the Russian local sheep breeds. A first split clustered the Russian sheep populations according to their wool type (fine wool, semi-fine wool and coarse wool). The Dagestan Mountain and Baikal fine-fleeced breeds differ from the other Merino-derived local breeds. The semi-fine wool cluster combined a breed of Romanian origin, Tsigai, with its derivative Altai Mountain, the two Romney-introgressed breeds Kuibyshev and North Caucasian, and the Lincoln-introgressed Russian longhaired breed. The coarse-wool group comprised the Nordic short-tailed Romanov, the long-fat-tailed outlier Kuchugur and two clusters of fat-tailed sheep: the Caucasian Mountain breeds and the Buubei, Karakul, Edilbai, Kalmyk and Tuva breeds. The Russian fat-tailed breeds shared co-ancestry with sheep from China and Southwestern Asia (Iran).

Conclusions

In this study, we derived the genetic characteristics of the major Russian local sheep breeds, which are moderately diverse and have a strong population structure. Pooling our data with a worldwide genotyping set gave deeper insight into the history and origin of the Russian sheep populations.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0399-5) contains supplementary material, which is available to authorized users.

Mitochondrial DNA variants of Podolian cattle breeds testify for a dual maternal origin

BACKGROUND

Over the past 15 years, 300 out of 6000 breeds of all farm animal species identified by the Food and Agriculture Organization of the United Nations (FAO) have gone extinct. Among cattle, many Podolian breeds are seriously endangered in various European areas. Podolian cattle include a group of very ancient European breeds, phenotypically close to the aurochs ancestors (Bos primigenius). The aim of the present study was to assess the genetic diversity of Podolian breeds and to reconstruct their origin.

METHODOLOGY

The mitochondrial DNA (mtDNA) control-regions of 18 Podolian breeds have been phylogenetically assessed. Nine non-Podolian breeds have been also included for comparison.

CONCLUSION

The overall analysis clearly highlights some peculiarities in the mtDNA gene pool of some Podolian breeds. In particular, a principal component analysis point to a genetic proximity between five breeds (Chianina, Marchigiana, Maremmana, Podolica Italiana and Romagnola) reared in Central Italy and the Turkish Grey. We here propose the suggestive hypothesis of a dual ancestral contribution to the present gene pool of Podolian breeds, one deriving from Eastern European cattle; the other arising from the arrival of Middle Eastern cattle into Central Italy through a different route, perhaps by sea, ferried by Etruscan boats. The historical migration of Podolian cattle from North Eastern Europe towards Italy has not cancelled the mtDNA footprints of this previous ancient migration.

Genomic diversity and population structure of three autochthonous Greek sheep breeds assessed with genome-wide DNA arrays

In the present study, genome-wide genotyping was applied to characterize the genetic diversity and population structure of three autochthonous Greek breeds: Boutsko, Karagouniko and Chios. Dairy sheep are among the most significant livestock species in Greece numbering approximately 9 million animals which are characterized by large phenotypic variation and reared under various farming systems. A total of 96 animals were genotyped with the Illumina's OvineSNP50K microarray beadchip, to study the population structure of the breeds and develop a specialized panel of single-nucleotide polymorphisms (SNPs), which could distinguish one breed from the others. Quality control on the dataset resulted in 46,125 SNPs, which were used to evaluate the genetic structure of the breeds. Population structure was assessed through principal component analysis (PCA) and admixture analysis, whereas inbreeding was estimated based on runs of homozygosity (ROHs) coefficients, genomic relationship matrix inbreeding coefficients (F_GRM) and patterns of linkage disequilibrium (LD). Associations between SNPs and breeds were analyzed with different inheritance models, to identify SNPs that distinguish among the breeds. Results showed high levels of genetic heterogeneity in the three breeds. Genetic distances among breeds were modest, despite their different ancestries. Chios and Karagouniko breeds were more genetically related to each other compared to Boutsko. Analysis revealed 3802 candidate SNPs that can be used to identify two-breed crosses and purebred animals. The present study provides, for the first time, data on the genetic background of three Greek indigenous dairy sheep breeds as well as a specialized marker panel that can be applied for traceability purposes as well as targeted genetic improvement schemes and conservation programs.

Population structure and genomic inbreeding in nine Swiss dairy cattle populations

BACKGROUND

Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF).

RESULTS

The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26.

CONCLUSIONS

The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.

New insights into the origin and the genetic status of the Balkan donkey from Serbia

The Balkan donkey (Equus asinus L.) is commonly regarded as a large-sized, unselected, unstructured and traditionally managed donkey breed. We assessed the current genetic status of the three largest E. asinus populations in the central Balkans (Serbia) by analysing the variability of nuclear microsatellites and the mitochondrial (mtDNA) control region of 77 and 49 individuals respectively. We further analysed our mtDNA dataset along with 209 published mtDNA sequences of ancient and modern individuals from 19 European and African populations to provide new insights into the origin and the history of the Balkan donkey. Serbian donkey populations are highly genetically diverse at both the nuclear and mtDNA levels despite severe population decline. Traditional Balkan donkeys in Serbia are rather heterogeneous; we found two groups of individuals with similar phenotypic features, somewhat distinct nuclear backgrounds and different proportions of mtDNA haplotypes belonging to matrilineal Clades 1 and 2. Another group, characterized by larger body size, different coat colour, distinct nuclear gene pool and predominantly Clade 2 haplotypes, was delineated as the Banat donkey breed. The maternal landscape of the large Balkan donkey population is highly heterogeneous and more complex than previously thought. Given the two independent domestication events in donkeys, multiple waves of introductions into the Balkans from Greece are hypothesized. Clade 2 donkeys probably appeared in Greece prior to those belonging to Clade 1, whereas expansion and diversification of Clade 1 donkeys within the Balkans predated that of Clade 2 donkeys.

Evaluation of Genetic Diversity, Population Structure, and Relationship Between Legendary Vechur Cattle and Crossbred Cattle of Kerala State, India

The legendary Vechur cattle of Kerala, described as a very short breed, and the crossbred (CB) Sunandini cattle population exhibited great phenotypic variation; hence, the present study attempted to analyze the genetic diversity existing between them. A set of 14 polymorphic microsatellites were chosen from FAO-ISAG panel and amplified from genomic DNA isolated from blood samples of 30 Vechur and 64 unrelated crossbred cattle, using fluorescent labeled primers. Both populations revealed high genetic diversity as evidenced from high observed number of alleles, Polymorphic Information Content and expected heterozygosity. Observed heterozygosity was lesser (0.699) than expected (0.752) in Vechur population which was further supported by positive F_IS value of 0.1149, indicating slight level of inbreeding in Vechur population. Overall, F_ST value was 0.065, which means genetic differentiation between crossbred and Vechur population was 6.5%, indicating that the crossbred cattle must have differentiated into a definite population that is different from the indigenous Vechur cows. Structure analysis indicated that the two populations showed distinct differences, with two underlying clusters. The present study supports the separation between Taurine and Zebu cattle and throws light onto the genetic diversity and relationship between native Vechur and crossbred cattle populations in Kerala state.

Genetic diversity and patterns of population structure in Creole goats from the Americas

Biodiversity studies are more efficient when large numbers of breeds belonging to several countries are involved, as they allow for an in-depth analysis of the within- and between-breed components of genetic diversity. A set of 21 microsatellites was used to investigate the genetic composition of 24 Creole goat breeds (910 animals) from 10 countries to estimate levels of genetic variability, infer population structure and understand genetic relationships among populations across the American continent. Three commercial transboundary breeds were included in the analyses to investigate admixture with Creole goats. Overall, the genetic diversity of Creole populations (mean number of alleles = 5.82 ± 1.14, observed heterozygosity = 0.585 ± 0.074) was moderate and slightly lower than what was detected in other studies with breeds from other regions. The Bayesian clustering analysis without prior information on source populations identified 22 breed clusters. Three groups comprised more than one population, namely from Brazil (Azul and Graúna; Moxotó and Repartida) and Argentina (Long and shorthair Chilluda, Pampeana Colorada and Angora-type goat). Substructure was found in Criolla Paraguaya. When prior information on sample origin was considered, 92% of the individuals were assigned to the source population (threshold q ≥ 0.700). Creole breeds are well-differentiated entities (mean coefficient of genetic differentiation = 0.111 ± 0.048, with the exception of isolated island populations). Dilution from admixture with commercial transboundary breeds appears to be negligible. Significant levels of inbreeding were detected (inbreeding coefficient > 0 in most Creole goat populations, P < 0.05). Our results provide a broad perspective on the extant genetic diversity of Creole goats, however further studies are needed to understand whether the observed geographical patterns of population structure may reflect the mode of goat colonization in the Americas.

Association of SNP variants of MHC Class II DRB gene with thermo-physiological traits in tropical goats

Host defense in vertebrates depend on many secreted regulatory proteins such as major histocompatibility complex (MHC) class II which provide important regulatory and effector functions of T cells. Gene polymorphism in the second exon of Capra-DRB gene in three major Nigerian goat breeds [West African Dwarf (WAD), Red Sokoto (RS), and Sahel (SH)] was analyzed by restriction fragment length polymorphisms (RFLP). Four restriction enzymes, BsaHI, AluI, HaeIII, and SacII, were utilized. The association between the polymorphic sites and some heat tolerance traits were also investigated in a total of 70 WAD, 90 RS, and 50 SH goats. Fourteen different types of alleles identified in the Nigerian goats, four of which were found in the peptide coding region (A57G, Q89R, G104D, and T112I), indicate a high degree of polymorphism at the DRB locus in this species. An obvious excess (P < 0.01) of non-synonymous substitutions than synonymous (dN/dS) in this locus is a reflection of adaptive evolution and positive selection. The phylogenetic trees revealed largely species-wise clustering in DRB gene. BsaHI, AluI, HaeIII, and SacII genotype frequencies were in Hardy-Weinberg equilibrium (P > 0.05), except AluI in RS goats and HaeIII in WAD goats (P < 0.05). The expected heterozygosity (H), which is a measure of gene diversity in the goat populations, ranged from 0.16 to 0.50. Genotypes AA (BsaHI), GG, GC and CC (AluI) and GG, GA, AA (HaeIII) appeared better in terms of heat tolerance. The heat-tolerant ability of SH and RS goats to the hot and humid tropical environment of Nigeria seemed better than that of the WAD goats. Sex effect (P < 0.05) was mainly on pulse rate and heat stress index, while there were varying interaction effects on heat tolerance. Variation at the DRB locus may prove to be important in possible selection and breeding for genetic resistance to heat stress in the tropics.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

Microsatellite diversity of the Nordic type of goats in relation to breed conservation: how relevant is pure ancestry?

In the last decades, several endangered breeds of livestock species have been re-established effectively. However, the successful revival of the Dutch and Danish Landrace goats involved crossing with exotic breeds and the ancestry of the current populations is therefore not clear. We have generated genotypes for 27 FAO-recommended microsatellites of these landraces and three phenotypically similar Nordic-type landraces and compared these breeds with central European, Mediterranean and south-west Asian goats. We found decreasing levels of genetic diversity with increasing distance from the south-west Asian domestication site with a south-east-to-north-west cline that is clearly steeper than the Mediterranean east-to-west cline. In terms of genetic diversity, the Dutch Landrace comes next to the isolated Icelandic breed, which has an extremely low diversity. The Norwegian coastal goat and the Finnish and Icelandic landraces are clearly related. It appears that by a combination of mixed origin and a population bottleneck, the Dutch and Danish Land-races are separated from the other breeds. However, the current Dutch and Danish populations with the multicoloured and long-horned appearance effectively substitute for the original breed, illustrating that for conservation of cultural heritage, the phenotype of a breed is more relevant than pure ancestry and the genetic diversity of the original breed. More in general, we propose that for conservation, the retention of genetic diversity of an original breed and of the visual phenotype by which the breed is recognized and defined needs to be considered separately.

Genetic variability of two Italian indigenous chicken breeds inferred from microsatellite marker analysis

The objective of this study was to determine the genetic structure and variability of Bionda Piemontese and Bianca di Saluzzo (Piedmont, Northwest Italy) using an international set of microsatellite loci (AVIANDIV-FAO). Differences compared with commercial lines and other Italian breeds were verified to justify the implementation of conservation programmes. Flock contribution to genetic variability was assessed following the approach implemented in the MolKin software. Comparison was performed using the fixation index and the Reynolds genetic distance. The most likely number of different populations was estimated using the clustering procedure implemented in STRUCTURE. The molecular information suggests that management practices could have prevented random mating and produced inbreeding and heterogeneity across flocks. In this respect, Bionda and Bianca show substructuring and are more similar to British breeds than other continental European breeds. Bionda and Bianca fit into the European breeds provided with the highest number of alleles and expected heterozygosity. There is a clear distinction between the Piedmont breeds and the other populations. The Piedmont poultry differ from both commercial lines and other Italian breeds and retain a high level of genetic variability. As for other indigenous breeds, Bionda and Bianca could make an original contribution to the industry in the future. A collective planned approach to restoration is essential, because the flocks are managed with poor regulation. Enhancing connection between breeders with an efficient replacement interchange and mating plan is the right way of controlling inbreeding, preventing substructuring and increasing variability within the flocks.