Genome-wide characteristics of copy number variation in Polish Holstein and Polish Red cattle using SNP genotyping assay

Genome-Wide Detection of Copy Number Variations and Their Potential Association with Carcass and Meat Quality Traits in Pingliang Red Cattle

Copy number variation (CNV) serves as a significant source of genetic diversity in mammals and exerts substantial effects on various complex traits. Pingliang red cattle, an outstanding indigenous resource in China, possess remarkable breeding value attributed to their tender meat and superior marbling quality. However, the genetic mechanisms influencing carcass and meat quality traits in Pingliang red cattle are not well understood. We generated a comprehensive genome-wide CNV map for Pingliang red cattle using the GGP Bovine 100K SNP chip. A total of 755 copy number variable regions (CNVRs) spanning 81.03 Mb were identified, accounting for approximately 3.24% of the bovine autosomal genome. Among these, we discovered 270 potentially breed-specific CNVRs in Pingliang red cattle, including 143 gains, 73 losses, and 54 mixed events. Functional annotation analysis revealed significant associations between these specific CNVRs and important traits such as carcass and meat quality, reproduction, exterior traits, growth traits, and health traits. Additionally, our network and transcriptome analysis highlighted CACNA2D1, CYLD, UBXN2B, TG, NADK, and ITGA9 as promising candidate genes associated with carcass weight and intramuscular fat deposition. The current study presents a genome-wide CNV map in Pingliang red cattle, highlighting breed-specific CNVRs, and transcriptome findings provide valuable insights into the underlying genetic characteristics of Pingliang red cattle. These results offer potential avenues for enhancing meat quality through a targeted breeding program.

Analysis of the effect of polymorphisms within the CATHL7 gene on dairy performance parameters

Introduction

Antimicrobial peptides, including cathelicidins, play a significant role in farm animals, influencing animal welfare, immunity, and thus the quality of animal products.

Material and Methods

The study used amplification-created restriction site and PCR-restriction fragment length polymorphism to analyse single nucleotide polymorphisms of the CATHL7 gene encoding the BMAP-34 protein in cattle, at positions 2,383 G > C and 2,468 G > C. The material was collected from 279 Polish Black-and-White Holstein-Friesian dairy cows.

Results

There were statistically significant differences between milk performance parameters in cows with the CATHL7/HhaI and CATHL7/HinfI genotypes. In the case of the CATHL7/HhaI polymorphism, the highest milk yield and protein and lactose content and the lowest somatic cell count in milk were observed for the CC genotype, while the fat content was the highest in milk from cows with the GG genotype. In the case of the CATHL7/HinfI polymorphism, the highest protein and lactose content in milk was observed for the CC genotype.

Conclusion

The results were statistically significant, which suggests that the search for relationships can be continued, and that the results can be used to improve selection programmes supporting dairy farming.

Genome-wide detection of copy number variations in Tharparkar cattle

Copy number variations (CNVs) are major forms of genetic variation with an increasing importance in animal genomics. This study used the Illumina BovineSNP 50 K BeadChip to detect the genome-wide CNVs in the Tharparkar cattle. With the aid of PennCNV software, we noticed a total of 447 copy number variation regions (CNVRs) across the autosomal genome, occupying nearly 2.17% of the bovine genome. The average size of detected CNVRs was found to be 122.2 kb, the smallest CNVR being 50.02 kb in size, to the largest being 1,232.87 Kb. Enrichment analyses of the genes in these CNVRs gave significant associations with molecular adaptation-related Gene Ontology (GO) terms. Most CNVR genes were significantly enriched for specific biological functions; signaling pathways, sensory responses to stimuli, and various cellular processes. In addition, QTL analysis of CNVRs described them to be linked with economically essential traits in cattle. The findings here provide crucial information for constructing a more comprehensive CNVR map for the indigenous cattle genome.

Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip

Background

Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks.

Results

A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits.

Conclusions

Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5759-1) contains supplementary material, which is available to authorized users.

Landscape of copy number variations in Bos taurus: individual - and inter-breed variability

Background

The number of studies of Copy Number Variation in cattle has increased in recent years. This has been prompted by the increased availability of data on polymorphisms and their relationship with phenotypes. In addition, livestock species are good models for some human phenotypes. In the present study, we described the landscape of CNV driven genetic variation in a large population of 146 individuals representing 13 cattle breeds, using whole genome DNA sequence.

Results

A highly significant variation among all individuals and within each breed was observed in the number of duplications (P < 10⁻¹⁵) and in the number of deletions (P < 10⁻¹⁵). We also observed significant differences between breeds for duplication (P = 0.01932) and deletion (P = 0.01006) counts. The same variation CNV length - inter-individual and inter-breed differences were significant for duplications (P < 10⁻¹⁵) and deletions (P < 10⁻¹⁵). Moreover, breed-specific variants were identified, with the largest proportion of breed-specific duplications (9.57%) found for Fleckvieh and breed-specific deletions found for Brown Swiss (5.00%). Such breed-specific CNVs were predominantly located in intragenic regions, however in Simmental, one deletion present in five individuals was found in the coding sequence of a novel gene ENSBTAG00000000688 on chromosome 18. In Brown Swiss, Norwegian Red and Simmental breed-specific deletions were located within KIT and MC1R genes, which are responsible for a coat colour. The functional annotation of coding regions underlying the breed-specific CNVs showed that in Norwegian Red, Guernsey, and Simmental significantly under- and overrepresented GO terms were related to chemical stimulus involved in sensory perception of smell and the KEGG pathways for olfactory transduction. In addition, specifically for the Norwegian Red breed, the dopaminergic synapse KEGG pathway was significantly enriched within deleted parts of the genome.

Conclusions

The CNV landscape in Bos taurus genome revealed by this study was highly complex, with inter-breed differences, but also a significant variation within breeds. The former, may explain some of the phenotypic differences among analysed breeds, and the latter contributes to within-breed variation available for selection.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4815-6) contains supplementary material, which is available to authorized users.

Genotyping-by-sequencing performance in selected livestock species

Application of next generation sequencing for large scale genotyping in livestock is limited by high costs and challenging data analysis process. However, available restriction enzyme-based enrichment techniques like e.g. genotyping-by-sequencing (GBS) are promising tools allowing reduction of financial outlies by a high sample multiplexing and narrowing down the sequenced genome areas to the randomly distributed read tags. In this study, we tested the performance of standard, PstI endonuclease-adapted GBS protocol for population genetics in cattle, horse and sheep with application of different, including low-depth sequencing setups. It was found that the detected SNPs display desirable polymorphism parameters and are evenly scattered across the whole genome including gene coding regions. It was also shown that the SNPs can be successfully applied in population genetics, revealing the genetic differentiation of the studied breeds. The GBS approach represents a cost-effective alternative to existing genotyping methods which may find adoption in various research applications.

Genomic landscape of copy number variation and copy neutral loss of heterozygosity events in equine sarcoids reveals increased instability of the sarcoid genome

Although they are the most common neoplasms in equids, sarcoids are not fully characterized at the molecular level. Therefore, the objective of this study was to characterize the landscape of structural rearrangements, such as copy number variation (CNV) and copy neutral loss of heterozygosity (cnLOH), in the genomes of sarcoid tumor cells. This information will not only broaden our understanding of the characteristics of this genome but will also improve the general knowledge of this tumor and the mechanisms involved in its generation. To this end, Equine SNP64K Illumina microarrays were applied along with bioinformatics tools dedicated for signal intensity analysis. The analysis revealed increased instability of the genome of sarcoid cells compared with unaltered skin tissue samples, which was manifested by the prevalence of CNV and cnLOH events. Many of the identified CNVs overlapped with the other research results, but the simultaneously observed variability in the number and sizes of detected aberrations indicated a need for further studies and the development of more reliable bioinformatics algorithms. The functional analysis of genes co-localized with the identified aberrations revealed that these genes are engaged in vital cellular processes. In addition, a number of these genes directly contribute to neoplastic transformation. Furthermore, large numbers of cnLOH events identified in the sarcoids suggested that they may play no less significant roles than CNVs in the carcinogenesis of this tumor. Thus, our results indicate the importance of cnLOH and CNV in equine sarcoid oncogenesis and present a direction of future research.

Copy number variations and genome-wide associations reveal putative genes and metabolic pathways involved with the feed conversion ratio in beef cattle

The use of genome-wide association results combined with other genomic approaches may uncover genes and metabolic pathways related to complex traits. In this study, the phenotypic and genotypic data of 1475 Nellore (Bos indicus) cattle and 941,033 single nucleotide polymorphisms (SNPs) were used for genome-wide association study (GWAS) and copy number variations (CNVs) analysis in order to identify candidate genes and putative pathways involved with the feed conversion ratio (FCR). The GWAS was based on the Bayes B approach analyzing genomic windows with multiple regression models to estimate the proportion of genetic variance explained by each window. The CNVs were detected with PennCNV software using the log R ratio and B allele frequency data. CNV regions (CNVRs) were identified with CNVRuler and a linear regression was used to associate CNVRs and the FCR. Functional annotation of associated genomic regions was performed with the Database for Annotation, Visualization and Integrated Discovery (DAVID) and the metabolic pathways were obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG). We showed five genomic windows distributed over chromosomes 4, 6, 7, 8, and 24 that explain 12 % of the total genetic variance for FCR, and detected 12 CNVRs (chromosomes 1, 5, 7, 10, and 12) significantly associated [false discovery rate (FDR) < 0.05] with the FCR. Significant genomic regions (GWAS and CNV) harbor candidate genes involved in pathways related to energetic, lipid, and protein metabolism. The metabolic pathways found in this study are related to processes directly connected to feed efficiency in beef cattle. It was observed that, even though different genomic regions and genes were found between the two approaches (GWAS and CNV), the metabolic processes covered were related to each other. Therefore, a combination of the approaches complement each other and lead to a better understanding of the FCR.

Next-generation sequencing as a powerful motor for advances in the biological and environmental sciences

Next-generation sequencing (NGS) provides unprecedented insight into (meta)genomes, (meta)transcriptomes (cDNA) and (meta)barcodes of individuals, populations and communities of Archaea, Bacteria and Eukarya, as well as viruses. This special issue combines reviews and original papers reporting technical and scientific advances in genomics and transcriptomics of non-model species, as well as quantification and functional analyses of biodiversity using NGS technologies of the second and third generations. In addition, certain papers also exemplify the transition from Sanger to NGS barcodes in molecular taxonomy.