Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Heat stress effects on milk yield traits and metabolites and mitigation strategies for dairy cattle breeds reared in tropical and sub-tropical countries

Heat stress is an important problem for dairy industry in many parts of the world owing to its adverse effects on productivity and profitability. Heat stress in dairy cattle is caused by an increase in core body temperature, which affects the fat production in the mammary gland. It reduces milk yield, dry matter intake, and alters the milk composition, such as fat, protein, lactose, and solids-not-fats percentages among others. Understanding the biological mechanisms of climatic adaptation, identifying and exploring signatures of selection, genomic diversity and identification of candidate genes for heat tolerance within indicine and taurine dairy breeds is an important progression toward breeding better dairy cattle adapted to changing climatic conditions of the tropics. Identifying breeds that are heat tolerant and their use in genetic improvement programs is crucial for improving dairy cattle productivity and profitability in the tropics. Genetic improvement for heat tolerance requires availability of genetic parameters, but these genetic parameters are currently missing in many tropical countries. In this article, we reviewed the HS effects on dairy cattle with regard to (1) physiological parameters; (2) milk yield and composition traits; and (3) milk and blood metabolites for dairy cattle reared in tropical countries. In addition, mitigation strategies such as physical modification of environment, nutritional, and genetic development of heat tolerant dairy cattle to prevent the adverse effects of HS on dairy cattle are discussed. In tropical climates, a more and cost-effective strategy to overcome HS effects is to genetically select more adaptable and heat tolerant breeds, use of crossbred animals for milk production, i.e., crosses between indicine breeds such as Gir, white fulani, N'Dama, Sahiwal or Boran to taurine breeds such as Holstein-Friesian, Jersey or Brown Swiss. The results of this review will contribute to policy formulations with regard to strategies for mitigating the effects of HS on dairy cattle in tropical countries.

Genomic signatures for drylands adaptation at gene-rich regions in African zebu cattle

Background

Indigenous Sudanese cattle are mainly indicine/zebu (humped) type. They thrive in the harshest dryland environments characterised by high temperatures, long seasonal dry periods, nutritional shortages, and vector disease challenges. Here, we sequenced 60 indigenous Sudanese cattle from six indigenous breeds and analysed the data using three genomic scan approaches to unravel cattle adaptation to the African dryland region.

Results

We identified a set of gene-rich selective sweep regions, detected mostly on chromosomes 5, 7 and 19, shared across African and Gir zebu. These include genes involved in immune response, body size and conformation, and heat stress response. We also identified selective sweep regions unique to Sudanese zebu. Of these, a 250 kb selective sweep on chromosome 16 spans seven genes, including PLCH2, PEX10, PRKCZ, and SKI, which are involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism.

Conclusions

Our results suggest that environmental adaptation may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control, in zebu cattle.

•

Sudanese cattle thrive in the harshest environments of the African drylands.

•

Bos indicus shared selected genes are involved in immune response, conformation, and heat stress response.

•

Sudanese zebu-specific sweep includes genes involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism.

•

Environmental adaptation in zebu cattle may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control.

A 16S Next Generation Sequencing Based Molecular and Bioinformatics Pipeline to Identify Processed Meat Products Contamination and Mislabelling

Simple Summary

Meat adulteration and fraud encompasses the deliberate fraudulent addition or substitution of proteins of animal or plant origin in edible products primarily for economic gain. The mitochondrial 16S ribosomal (rRNA) gene was used to identify species that are present in pure and processed meat samples. The meat samples were sequenced using an Illumina sequencing platform, and bioinformatics analysis was carried out for species identification. The results indicated that pork was the major contaminant in most of the meat samples. The bioinformatics pipeline demonstrated its specificity through identification of species specific and quantification of the contamination levels across all samples. Food business operators and regulatory sectors can validate this method for food fraud checks and manage any form of mislabeling in the animal or plant protein food ecosystem.

Abstract

Processed meat is a target in meat adulteration for economic gain. This study demonstrates a molecular and bioinformatics diagnostic pipeline, utilizing the mitochondrial 16S ribosomal RNA (rRNA) gene, to determine processed meat product mislabeling through Next-Generation Sequencing. Nine pure meat samples were collected and artificially mixed at different ratios to verify the specificity and sensitivity of the pipeline. Processed meat products (n = 155), namely, minced meat, biltong, burger patties, and sausages, were collected across South Africa. Sequencing was performed using the Illumina MiSeq sequencing platform. Each sample had paired-end reads with a length of ±300 bp. Quality control and filtering was performed using BBDuk (version 37.90a). Each sample had an average of 134,000 reads aligned to the mitochondrial genomes using BBMap v37.90. All species in the artificial DNA mixtures were detected. Processed meat samples had reads that mapped to the Bos (90% and above) genus, with traces of reads mapping to Sus and Ovis (2–5%) genus. Sausage samples showed the highest level of contamination with 46% of the samples having mixtures of beef, pork, or mutton in one sample. This method can be used to authenticate meat products, investigate, and manage any form of mislabeling.

Identification and in silico characterization of structural and functional impacts of genetic variants in milk protein genes in the Zebu breeds Guzerat and Gyr

Whole genome sequencing of bovine breeds has allowed identification of genetic variants in milk protein genes. However, functional repercussion of such variants at a molecular level has seldom been investigated. Here, the results of a multistep Bioinformatic analysis for functional characterization of recently identified genetic variants in Brazilian Gyr and Guzerat breeds is described, including predicted effects on the following: (i) evolutionary conserved nucleotide positions/regions; (ii) protein function, stability, and interactions; (iii) splicing, branching, and miRNA binding sites; (iv) promoters and transcription factor binding sites; and (v) collocation with QTL. Seventy-one genetic variants were identified in the caseins (CSN1S1, CSN2, CSN1S2, and CSN3), LALBA, LGB, and LTF genes. Eleven potentially regulatory variants and two missense mutations were identified. LALBA Ile60Val was predicted to affect protein stability and flexibility, by reducing the number the disulfide bonds established. LTF Thr546Asn is predicted to generate steric clashes, which could mildly affect iron coordination. In addition, LALBA Ile60Val and LTF Thr546Asn affect exonic splicing enhancers and silencers. Consequently, both mutations have the potential of affecting immune response at individual level, not only in the mammary gland. Although laborious, this multistep procedure for classifying variants allowed the identification of potentially functional variants for milk protein genes.

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.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data

BACKGROUND

The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics.

RESULTS

A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds.

CONCLUSION

The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

Reduced representation approach for identification of genome-wide SNPs and their annotation for economically important traits in Indian Tharparkar cattle

The present study was carried out in Tharparkar cattle for identification of genome-wide SNPs and microsatellites, and then annotate the identified high-quality SNPs to milk production, fertility, carcass, adaptability and immune response of economically important traits. A total of 146,011 SNPs were identified with respect to Bos taurus reference genome which are indicus specific, out of which 10,519 SNPs were found to be novel. Similarly, a total of 87,047 SNPs were identified with respect to Bos indicus reference genome. After final annotation of SNPs identified with respect to Bos indicus reference genome, 2871 SNPs were found to be associated in 383 candidate genes having to do with milk production, fertility, carcass, immune response and adaptability traits. Following that, 2571 microsatellites were identified. The information mined from the data might be of importance for the future breed improvement programs, conservation efforts and for enhancing the SNPs density of the existing bovine SNP chips.

Genome-wide association study of milk and reproductive traits in dual-purpose Xinjiang Brown cattle

Background

Dual-purpose cattle are more adaptive to environmental challenges than single-purpose dairy or beef cattle. Balance among milk, reproductive, and mastitis resistance traits in breeding programs is therefore more critical for dual-purpose cattle to increase net income and maintain well-being. With dual-purpose Xinjiang Brown cattle adapted to the Xinjiang Region in northwestern China, we conducted genome-wide association studies (GWAS) to dissect the genetic architecture related to milk, reproductive, and mastitis resistance traits. Phenotypic data were collected for 2410 individuals measured during 1995–2017. By adding another 445 ancestors, a total of 2855 related individuals were used to derive estimated breeding values for all individuals, including the 2410 individuals with phenotypes. Among phenotyped individuals, we genotyped 403 cows with the Illumina 150 K Bovine BeadChip.

Results

GWAS were conducted with the FarmCPU (Fixed and random model circulating probability unification) method. We identified 12 markers significantly associated with six of the 10 traits under the threshold of 5% after a Bonferroni multiple test correction. Seven of these SNPs were in QTL regions previously identified to be associated with related traits. One identified SNP, BovineHD1600006691, was significantly associated with both age at first service and age at first calving. This SNP directly overlapped a QTL previously reported to be associated with calving ease. Within 160 Kb upstream and downstream of each significant SNP identified, we speculated candidate genes based on functionality. Four of the SNPs were located within four candidate genes, including CDH2, which is linked to milk fat percentage, and GABRG2, which is associated with milk protein yield.

Conclusions

These findings are beneficial not only for breeding through marker-assisted selection, but also for genome editing underlying the related traits to enhance the overall performance of dual-purpose cattle.

Genome-Wide SNPs and InDels Characteristics of Three Chinese Cattle Breeds

Simple Summary

Whole-genome resequencing is an important tool to reveal the in-depth genomic characteristics of a genome. Adaptability traits are key to the survival of the south Chinese zebu cattle. However, the potential genetic information behind these remarkable traits still remains uncertain and needs to be addressed. In the current study, we utilized a total of 15 local south Chinese cattle samples (Leiqiong (LQ), Wannan (WN), Wenshan (WS)) from one of our previous studies mapped to the old reference genome (Btau_5.0.1) and remapped them to the latest reference genome (ARS-UCD1.2) to explore potential single nucleotide polymorphisms (SNPs) and insertions-deletions (InDels) responsible for some important immune related traits. The present study emphasizes and illustrates the genetic diversity, extending our previous study. The InDel annotation show that WS cattle had more enriched genes associated with immune functions than the other two breeds. Our findings provide valuable resources for further investigation of the functions of SNP- and InDels-related genes and help to determine the molecular basis of adaptive mutations in Chinese zebu cattle.

Abstract

We report genome characterization of three native Chinese cattle breeds discovering ~34.3 M SNPs and ~3.8 M InDels using whole genome resequencing. On average, 10.4 M SNPs were shared amongst the three cattle breeds, whereas, 3.0 M, 4.9 M and 5.8 M were specific to LQ, WN and WS breeds, respectively. Gene ontology (GO)analysis revealed four immune response-related GO terms were over represented in all samples, while two immune signaling pathways were significantly over-represented in WS cattle. Altogether, we found immune related genes (PGLYRP2, ROMO1, FYB2, CD46, TSC1) in the three cattle breeds. Our study provides insights into the genetic basis of Chinese indicine adaptation to the tropic and subtropical environment, and provides a valuable resource for further investigations of genetic characteristics of the three breeds.

Genomic variants identified from whole-genome resequencing of indicine cattle breeds from Pakistan

The primary goal of cattle genomics is the identification of genome-wide polymorphism associated with economically important traits. The bovine genome sequencing project was completed in 2009. Since then, using massively parallel sequencing technologies, a large number of Bos taurus cattle breeds have been resequenced and scanned for genome-wide polymorphisms. As a result, a substantial number of single nucleotide polymorphisms (SNPs) have been discovered across European Bos taurus genomes, whereas extremely less number of SNPs are cataloged for Bos indicus breeds. In this study, we performed whole-genome resequencing, reference-based mapping, functional annotation and gene enrichment analysis of 20 sires representing eleven important Bos indicus (indicine) breeds of Pakistan. The breeds sequenced here include: Sahiwal, Red Sindhi, Tharparkar and Cholistani (tropically adapted dairy and dual purpose breeds), Achai, Bhagnari, Dajal and Lohani (high altitude adapted dual and drought purpose breeds); Dhanni, Hisar Haryana and Gabrali (dairy and light drought purpose breeds). A total of 17.4 billion QC passed reads were produced using BGISEQ-500 next generation sequencing platform to generate 9 to 27-fold genome coverage (average ~16×) for each of the 20 sequenced sires. A total of 67,303,469 SNPs were identified, of which 3,850,365 were found novel and 1,083,842 insertions-deletions (InDels) were detected across the whole sequenced genomes (491,247 novel). Comparative analysis using coding region SNPs revealed a close relationship between the best milking indicine breeds; Red Sindhi and Sahiwal. On the other hand, Bhagnari and Tharparkar being popular for their adaptation to dry and extremely hot climates were found to share the highest number of SNPs. Functional annotation identified a total of 3,194 high-impact (disruptive) SNPs and 745 disruptive InDels (in 275 genes) that may possibly affect economically important dairy and beef traits. Functional enrichment analysis was performed and revealed that high or moderate impact variants in wingless-related integration site (Wnt) and vascular smooth muscle contraction (VSMC) signaling pathways were significantly over-represented in tropically adapted heat tolerant Pakistani-indicine breeds. On the other hand, vascular endothelial growth factor (VEGF) and hypoxia-inducible factor 1 (HIF-1) signaling pathways were found over-represented in highland adapted Pakistani-indicine breeds. Similarly, the ECM-receptor interaction and Jak-STAT signaling pathway were significantly enriched in dairy and beef purpose Pakistani-indicine cattle breeds. The Toll-like receptor signaling pathway was significantly enriched in most of the Pakistani-indicine cattle. Therefore, this study provides baseline data for further research to investigate the molecular mechanisms of major traits and to develop potential genomic markers associated with economically important breeding traits, particularly in indicine cattle.

Whole-Genome Sequencing of Three Native Cattle Breeds Originating From the Northernmost Cattle Farming Regions

Northern Fennoscandia and the Sakha Republic in the Russian Federation represent the northernmost regions on Earth where cattle farming has been traditionally practiced. In this study, we performed whole-genome sequencing to genetically characterize three rare native breeds Eastern Finncattle, Western Finncattle and Yakutian cattle adapted to these northern Eurasian regions. We examined the demographic history, genetic diversity and unfolded loci under natural or artificial selection. On average, we achieved 13.01-fold genome coverage after mapping the sequencing reads on the bovine reference genome (UMD 3.1) and detected a total of 17.45 million single nucleotide polymorphisms (SNPs) and 1.95 million insertions-deletions (indels). We observed that the ancestral species (Bos primigenius) of Eurasian taurine cattle experienced two notable prehistorical declines in effective population size associated with dramatic climate changes. The modern Yakutian cattle exhibited a higher level of within-population variation in terms of number of SNPs and nucleotide diversity than the contemporary European taurine breeds. This result is in contrast to the results of marker-based cattle breed diversity studies, indicating assortment bias in previous analyses. Our results suggest that the effective population size of the ancestral Asiatic taurine cattle may have been higher than that of the European cattle. Alternatively, our findings could indicate the hybrid origins of the Yakutian cattle ancestries and possibly the lack of intensive artificial selection. We identified a number of genomic regions under selection that may have contributed to the adaptation to the northern and subarctic environments, including genes involved in disease resistance, sensory perception, cold adaptation and growth. By characterizing the native breeds, we were able to obtain new information on cattle genomes and on the value of the adapted breeds for the conservation of cattle genetic resources.

High Levels of Variation Within Gene Sequences of Olea europaea L

Gene sequence variation in cultivated olive (Olea europaea L. subsp. europaea var. europaea), the most important oil tree crop of the Mediterranean basin, has been poorly evaluated up to now. A deep sequence analysis of fragments of four genes, OeACP1, OeACP2, OeLUS and OeSUT1, in 90 cultivars, revealed a wide range of polymorphisms along all recognized allele forms and unexpected allele frequencies and genotype combinations. High linkage values among most polymorphisms were recorded within each gene fragment. The great sequence variability corresponded to a low number of alleles and, surprisingly, to a small fraction of genotype combinations. The distribution, frequency, and combination of the different alleles at each locus is possibly due to natural and human pressures, such as selection, ancestrality, or fitness. Phylogenetic analyses of allele sequences showed distant and complex patterns of relationships among cultivated olives, intermixed with other related forms, highlighting an evolutionary connection between olive cultivars and the O. europaea subspecies cuspidata and cerasiformis. This study demonstrates how a detailed and complete sequence analysis of a few gene portions and a thorough genotyping on a representative set of cultivars can clarify important issues related to sequence polymorphisms, reconstructing the phylogeny of alleles, as well as the genotype combinations. The identification of regions representing blocks of recombination could reveal polymorphisms that represent putatively functional markers. Indeed, specific mutations found on the analyzed OeACP1 and OeACP2 fragments seem to be correlated to the fruit weight.

Sequencing the mosaic genome of Brahman cattle identifies historic and recent introgression including polled

Brahman cattle have a Bos indicus and Bos taurus mosaic genome, as a result of the process used to create the breed (repeat backcrossing of Bos taurus females to Bos indicus bulls). With the aim of identifying Bos taurus segments in the Brahman genome at sequence level resolution, we sequenced the genomes of 46 influential Brahman bulls. Using 36 million variants identified in the sequences, we searched for regions close to fixation for Bos indicus or Bos taurus segments that were longer than expected by chance (from simulation of the breed formation history of Brahman cattle). Regions close to fixation for Bos indicus content were enriched for protein synthesis genes, while regions of higher Bos taurus content included genes of the G-protein coupled receptor family (including genes implicated in puberty, such as THRS). The region with the most extreme Bos taurus enrichment was on chromosome 14 surrounding PLAG1. The introgressed Bos taurus allele at PLAG1 increases stature and the high frequency of the allele likely reflects strong selection for the trait. Finally, we provide evidence that the polled mutation in Brahmans, a desirable trait under very strong recent selection, is of Celtic origin and is introgressed from Bos taurus.

Haplotype Loci Under Selection in Canadian Durum Wheat Germplasm Over 60 Years of Breeding: Association With Grain Yield, Quality Traits, Protein Loss, and Plant Height

Durum wheat was introduced in the southern prairies of western Canada in the late nineteenth century. Breeding efforts have mainly focused on improving quality traits to meet the pasta industry demands. For this study, 192 durum wheat lines were genotyped using the Illumina 90K Infinium iSelect assay, and resulted in a total of 14,324 polymorphic SNPs. Genetic diversity changed over time, declining during the first 20 years of breeding in Canada, then increased in the late 1980s and early 1990s. We scanned the genome for signatures of selection, using the total variance Fst-based outlier detection method (Lositan), the hierarchical island model (Arlequin) and the Bayesian genome scan method (BayeScan). A total of 407 outliers were identified and clustered into 84 LD-based haplotype loci, spanning all 14 chromosomes of the durum wheat genome. The association analysis detected 54 haplotype loci, of which 39% contained markers with a complete reversal of allelic state. This tendency to fixation of favorable alleles corroborates the success of the Canadian durum wheat breeding programs over time. Twenty-one haplotype loci were associated with multiple traits. In particular, hap_4B_1 explained 20.6, 17.9 and 16.6% of the phenotypic variance of pigment loss, pasta b∗ and dough extensibility, respectively. The locus hap_2B_9 explained 15.9 and 17.8% of the variation of protein content and protein loss, respectively. All these pleiotropic haplotype loci offer breeders the unique opportunity for further improving multiple traits, facilitating marker-assisted selection in durum wheat, and could help in identifying genes as functional annotations of the wheat genome become available.

Genetic variants analysis of three dromedary camels using whole genome sequencing data

Whole genome wide identification and annotation of genetic variations in camels is in its first steps. The aim of this study was the identification of genome wide variants, functional annotations of them and enrichment analysis of affected genes using whole genome sequencing data of three dromedary camels. The genomes of two Iranian female dromedary camels that mostly used to produce meat and milk were sequenced to 41.9-fold and 38.6-fold coverage. A total of 4,727,238 single-nucleotide polymorphisms (SNPs) and 692,908 indels (insertions and deletions) were found by mapping raw reads to the dromedary reference assembly (GenBank Accession: GCA_000767585.1). In-silico functional annotation of the discovered variants in under study samples revealed that most SNPs (2,305,738; 48.78%) and indels (339,756; 49.03%) were located in intergenic regions. A comparison of the identified SNPs with those of the African camel (BioProject Accession: PRJNA269274) indicated that they had 993,474 SNPs in common. We found 15,168 non-synonymous SNPs in the shared variants of the three camels that could affect gene function and protein structure. Obtained results revealed that there were 7085, 6271 and 4688 non-synonymous SNPs among the 3436, 3058 and 2882 genes in the specific gene sets of Yazd dromedary, Trod dromedary and African dromedary, respectively. The list of genes predicted to be affected by non-synonymous variants in different individuals was subjected to gene ontology (GO) enrichment analysis.

Assessing genetic architecture and signatures of selection of dual purpose Gir cattle populations using genomic information

Gir is one of the main cattle breeds raised in tropical South American countries. Strong artificial selection through its domestication resulted in increased genetic differentiation among the countries in recent years. Over the years, genomic studies in Gir have become more common. However, studies of population structure and signatures of selection in divergent Gir populations are scarce and need more attention to better understand genetic differentiation, gene flow, and genetic distance. Genotypes of 173 animals selected for growth traits and 273 animals selected for milk production were used in this study. Clear genetic differentiation between beef and dairy populations was observed. Different criteria led to genetic divergence and genetic differences in allele frequencies between the two populations. Gene segregation in each population was forced by artificial selection, promoting isolation, and increasing genetic variation between them. Results showed evidence of selective forces in different regions of the genome. A total of 282 genes were detected under selection in the test population based on the fixation index (Fst), integrated haplotype score (iHS), and cross-population extend haplotype homozygosity (XP-EHH) approaches. The QTL mapping identified 35 genes associated with reproduction, milk composition, growth, meat and carcass, health, or body conformation traits. The investigation of genes and pathways showed that quantitative traits associated to fertility, milk production, beef quality, and growth were involved in the process of differentiation of these populations. These results would support further investigations of population structure and differentiation in the Gir breed.

Signatures of positive selection in African Butana and Kenana dairy zebu cattle

Butana and Kenana are two types of zebu cattle found in Sudan. They are unique amongst African indigenous zebu cattle because of their high milk production. Aiming to understand their genome structure, we genotyped 25 individuals from each breed using the Illumina BovineHD Genotyping BeadChip. Genetic structure analysis shows that both breeds have an admixed genome composed of an even proportion of indicine (0.75 ± 0.03 in Butana, 0.76 ± 0.006 in Kenana) and taurine (0.23 ± 0.009 in Butana, 0.24 ± 0.006 in Kenana) ancestries. We also observe a proportion of 0.02 to 0.12 of European taurine ancestry in ten individuals of Butana that were sampled from cattle herds in Tamboul area suggesting local crossbreeding with exotic breeds. Signatures of selection analyses (iHS and Rsb) reveal 87 and 61 candidate positive selection regions in Butana and Kenana, respectively. These regions span genes and quantitative trait loci (QTL) associated with biological pathways that are important for adaptation to marginal environments (e.g., immunity, reproduction and heat tolerance). Trypanotolerance QTL are intersecting candidate regions in Kenana cattle indicating selection pressure acting on them, which might be associated with an unexplored level of trypanotolerance in this cattle breed. Several dairy traits QTL are overlapping the identified candidate regions in these two zebu cattle breeds. Our findings underline the potential to improve dairy production in the semi-arid pastoral areas of Africa through breeding improvement strategy of indigenous local breeds.

Climate Adaptation of Tropical Cattle

There is sustained growth in the number of tropical cattle, which represent more than half of all cattle worldwide. By and large, most research in tropical areas is still focused on breeds of cattle, their particular advantages or disadvantages in tropical areas, and the tropical forages or feeds that could be usefully fed to them. A consistent issue for adaptation to climate is the heat of tropical environments. Changing the external characteristics of the animal, such as color and coat characteristics, is one way to adapt, and there are several major genes for these traits. However, further improvement in heat tolerance and other adaptation traits will need to use the entire genome and all physical and physiological systems. Apart from the response to heat, climate forcing through methane emission identifies dry season weight loss as an important if somewhat neglected trait in climate adaptation of cattle. The use of genome-estimated breeding values in tropical areas is in its infancy and will be difficult to implement, but will be essential for rapid, coordinated genetic improvement. The difficulty of implementation cannot be exaggerated and may require major improvements in methodology.

Signatures of Selection for Environmental Adaptation and Zebu × Taurine Hybrid Fitness in East African Shorthorn Zebu

The East African Shorthorn Zebu (EASZ) cattle are ancient hybrid between Asian zebu × African taurine cattle preferred by local farmers due to their adaptability to the African environment. The genetic controls of these adaptabilities are not clearly understood yet. Here, we genotyped 92 EASZ samples from Kenya (KEASZ) with more than 770,000 SNPs and sequenced the genome of a pool of 10 KEASZ. We observe an even admixed autosomal zebu × taurine genomic structure in the population. A total of 101 and 165 candidate regions of positive selection, based on genome-wide SNP analyses (meta-SS, Rsb, iHS, and ΔAF) and pooled heterozygosity (Hp) full genome sequence analysis, are identified, in which 35 regions are shared between them. A total of 142 functional variants, one novel, have been detected within these regions, in which 30 and 26 were classified as of zebu and African taurine origins, respectively. High density genome-wide SNP analysis of zebu × taurine admixed cattle populations from Uganda and Nigeria show that 25 of these regions are shared between KEASZ and Uganda cattle, and seven regions are shared across the KEASZ, Uganda, and Nigeria cattle. The identification of common candidate regions allows us to fine map 18 regions. These regions intersect with genes and QTL associated with reproduction and environmental stress (e.g., immunity and heat stress) suggesting that the genome of the zebu × taurine admixed cattle has been uniquely selected to maximize hybrid fitness both in terms of reproduction and survivability.

Single nucleotide variants and InDels identified from whole-genome re-sequencing of Guzerat, Gyr, Girolando and Holstein cattle breeds

Whole-genome re-sequencing, alignment and annotation analyses were undertaken for 12 sires representing four important cattle breeds in Brazil: Guzerat (multi-purpose), Gyr, Girolando and Holstein (dairy production). A total of approximately 4.3 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.7 to 16.4-fold genome coverage. A total of 27,441,279 single nucleotide variations (SNVs) and 3,828,041 insertions/deletions (InDels) were detected in the samples, of which 2,557,670 SNVs and 883,219 InDels were novel. The submission of these genetic variants to the dbSNP database significantly increased the number of known variants, particularly for the indicine genome. The concordance rate between genotypes obtained using the Bovine HD BeadChip array and the same variants identified by sequencing was about 99.05%. The annotation of variants identified numerous non-synonymous SNVs and frameshift InDels which could affect phenotypic variation. Functional enrichment analysis was performed and revealed that variants in the olfactory transduction pathway was over represented in all four cattle breeds, while the ECM-receptor interaction pathway was over represented in Girolando and Guzerat breeds, the ABC transporters pathway was over represented only in Holstein breed, and the metabolic pathways was over represented only in Gyr breed. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Gyr, Girolando, Guzerat and Holstein breeding programs.

Whole genome sequencing of Guzerá cattle reveals genetic variants in candidate genes for production, disease resistance, and heat tolerance

In bovines, artificial selection has produced a large number of breeds which differ in production, environmental adaptation, and health characteristics. To investigate the genetic basis of these phenotypical differences, several bovine breeds have been sequenced. Millions of new SNVs were described at every new breed sequenced, suggesting that every breed should be sequenced. Guzerat or Guzerá is an indicine breed resistant to drought and parasites that has been the base for some important breeds such as Brahman. Here, we describe the sequence of the Guzerá genome and the in silico functional analyses of intragenic breed-specific variations. Mate-paired libraries were generated using the ABI SOLiD system. Sequences were mapped to the Bos taurus reference genome (UMD 3.1) and 87% of the reference genome was covered at a 26X. Among the variants identified, 2,676,067 SNVs and 463,158 INDELs were homozygous, not found in any database searched, and may represent true differences between Guzerá and B. taurus. Functional analyses investigated with the NGS-SNP package focused on 1069 new, non-synonymous SNVs, splice-site variants (including acceptor and donor sites, and the conserved regions at both intron borders, referred to here as splice regions) and coding INDELs (NS/SS/I). These NS/SS/I map to 935 genes belonging to cell communication, environmental adaptation, signal transduction, sensory, and immune systems pathways. These pathways have been involved in phenotypes related to health, adaptation to the environment and behavior, and particularly, disease resistance and heat tolerance. Indeed, 105 of these genes are known QTLs for milk, meat and carcass, production, reproduction, and health traits. Therefore, in addition to describing new genetic variants, our approach provided groundwork for unraveling key candidate genes and mutations.

Genome-wide analysis reveals signatures of selection for important traits in domestic sheep from different ecoregions

Background

Throughout a long period of adaptation and selection, sheep have thrived in a diverse range of ecological environments. Mongolian sheep is the common ancestor of the Chinese short fat-tailed sheep. Migration to different ecoregions leads to changes in selection pressures and results in microevolution. Mongolian sheep and its subspecies differ in a number of important traits, especially reproductive traits. Genome-wide intraspecific variation is required to dissect the genetic basis of these traits.

Results

This research resequenced 3 short fat-tailed sheep breeds with a 43.2-fold coverage of the sheep genome. We report more than 17 million single nucleotide polymorphisms and 2.9 million indels and identify 143 genomic regions with reduced pooled heterozygosity or increased genetic distance to each other breed that represent likely targets for selection during the migration. These regions harbor genes related to developmental processes, cellular processes, multicellular organismal processes, biological regulation, metabolic processes, reproduction, localization, growth and various components of the stress responses. Furthermore, we examined the haplotype diversity of 3 genomic regions involved in reproduction and found significant differences in TSHR and PRL gene regions among 8 sheep breeds.

Conclusions

Our results provide useful genomic information for identifying genes or causal mutations associated with important economic traits in sheep and for understanding the genetic basis of adaptation to different ecological environments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3212-2) contains supplementary material, which is available to authorized users.

Evolutionary Analyses of Hanwoo (Korean Cattle)-Specific Single-Nucleotide Polymorphisms and Genes Using Whole-Genome Resequencing Data of a Hanwoo Population

Advances in next generation sequencing (NGS) technologies have enabled population-level studies for many animals to unravel the relationships between genotypic differences and traits of specific populations. The objective of this study was to perform evolutionary analysis of single nucleotide polymorphisms (SNP) in genes of Korean native cattle Hanwoo in comparison to SNP data from four other cattle breeds (Jersey, Simmental, Angus, and Holstein) and four related species (pig, horse, human, and mouse) obtained from public databases through NGS-based resequencing. We analyzed population structures and differentiation levels for the five cattle breeds and estimated species-specific SNPs with their origins and phylogenetic relationships among species. In addition, we identified Hanwoo-specific genes and proteins, and determined distinct changes in protein-protein interactions among five species (cattle, pig, horse, human, mouse) in the STRING network database by additionally considering indirect protein interactions. We found that the Hanwoo population was clearly different from the other four cattle populations. There were Hanwoo-specific genes related to its meat trait. Protein interaction rewiring analysis also confirmed that there were Hanwoo-specific protein-protein interactions that might have contributed to its unique meat quality.

Identification of novel SNPs in differentially expressed genes and its association with horn cancer of Bos indicus bullocks by next-generation sequencing

The use of polymorphic markers like SNPs promises to provide comprehensive tool for analysing genome and identifying genomic regions that contribute to cancer phenotype. Horn cancer is the most common cancer among Bos indicus animals. Increased expression of some genes due to polymorphisms increases risk of HC incidence. We successfully amplified 91 SNPs located in 69 genes in 52 samples, each of HC and HN. Equimolar concentration of amplicons from 69 PCR products of each sample was pooled and subjected to sequencing using Ion Torrent PGM. Data obtained were analysed using DNASTAR software package and case control analysis using SAS software. We found SNP present in BPIFA1 gene of B. indicus shows association with event of HC which reflects its potential to be a genetic marker. Bioinformatic analysis to detect structural and functional impact nsSNP of BPIFA1 added another layer of confirmation to our result. We successfully identified SNP associated with HC as well as demonstrated efficient approach for limited number of SNP discovery and validation in targeted genomics regions in large number of samples combining PCR amplification and Ion Torrent PGM sequencing which suits small-scale laboratories with limited budget.

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.

Use of Genomic Tools to Improve Cattle Health in the Context of Infectious Diseases

Although infectious diseases impose a heavy economic burden on the cattle industry, the etiology of many disorders that affect livestock is not fully elucidated, and effective countermeasures are often lacking. The main tools available until now have been vaccines, antibiotics and antiparasitic drugs. Although these have been very successful in some cases, the appearance of parasite and microbial resistance to these treatments is a cause of concern. Next-generation sequencing provides important opportunities to tackle problems associated with pathogenic illnesses. This review describes the rapid gains achieved to track disease progression, identify the pathogens involved, and map pathogen interactions with the host. Use of novel genomic tools subsequently aids in treatment development, as well as successful creation of breeding programs aimed toward less susceptible livestock. These may be important tools for mitigating the long term effects of combating infection and helping reduce the reliance on antibiotic treatment.

Whole-genome sequencing of the endangered bovine species Gayal (Bos frontalis) provides new insights into its genetic features

Gayal (Bos frontalis) is a semi-wild and endangered bovine species that differs from domestic cattle (Bos taurus and Bos indicus), and its genetic background remains unclear. Here, we performed whole-genome sequencing of one Gayal for the first time, with one Red Angus cattle and one Japanese Black cattle as controls. In total, 97.8 Gb of sequencing reads were generated with an average 11.78-fold depth and >98.44% coverage of the reference sequence (UMD3.1). Numerous different variations were identified, 62.24% of the total single nucleotide polymorphisms (SNPs) detected in Gayal were novel, and 16,901 breed-specific nonsynonymous SNPs (BS-nsSNPs) that might be associated with traits of interest in Gayal were further investigated. Moreover, the demographic history of bovine species was first analyzed, and two population expansions and two population bottlenecks were identified. The obvious differences among their population sizes supported that Gayal was not B. taurus. The phylogenic analysis suggested that Gayal was a hybrid descendant from crossing of male wild gaur and female domestic cattle. These discoveries will provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds and may assist relevant departments with future conservation and utilization of Gayal.

Genome-Wide Study of Structural Variants in Bovine Holstein, Montbéliarde and Normande Dairy Breeds

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.

Sequence analysis of single-copy genes in two wild olive subspecies: nucleotide diversity and potential use for testing admixture

The wild olive distribution extends from the Mediterranean region to south Asia and Austral Africa. The species is also invasive, particularly in Australia. Here, we investigated the sequence variation at five nuclear single-copy genes in 41 native and invasive accessions of the Mediterranean and African olive subspecies. The nucleotide diversity was assessed and the phylogenetic relationships between alleles were depicted with haplotype networks. A Bayesian clustering method (STRUCTURE) was applied to identify the main gene pools. We found an average of 18.4 alleles per locus. Native Mediterranean and African olives only share one allele, which testifies for ancient admixture on the Red Sea hills. The presence of divergent alleles in the Mediterranean olive, as well as the identification of two main genetic clusters, suggests a complex origin with two highly differentiated gene pools from the eastern and western Mediterranean that recently admixed. In the invasive range, relatively high nucleotide diversity is observed as a consequence of the introduction of alleles from two subspecies. Our data confirm that four invasive individuals are early-generation hybrids. Finally, the utility of single-copy gene sequences in olive population genomic and phylogenetic studies is briefly discussed.

Livestock and the promise of genomics

The emergence of the middle class in countries such as Brazil, Russia, India, and China is resulting in increasing global demand for animal-based food products. This increase represents a unique opportunity for Canadian livestock producers to export their products to new markets and expand Canada's reputation as a global provider of safe and highest quality food items. This article has two major themes. First, current Canadian contributions to livestock genomics in the cattle and swine industries are outlined. Second, important future opportunities are discussed, including the high throughput collection of phenotypic data, development of environmentally friendly livestock, emergence of decision support software, and the use of Web 2.0. Through the use of genomic technologies, livestock producers can not only ensure that the nutritional demands of Canada are secured, but also play a pivotal role in ensuring the rest of the world is fed as well. Furthermore, investment through initiatives led by Genome Canada has ensured that Canada is favorably positioned to contribute cutting-edge solutions to meet this global challenge. Ultimately, genomic-based innovations will enable producers to increase efficiency, lower production costs, decrease the use of prophylactics, and limit the expenditure of resources.