Genome-wide copy number variant analysis reveals variants associated with 10 diverse production traits in Holstein cattle

Uncovering structural variants in Creole cattle from Guadeloupe and their impact on environmental adaptation through whole genome sequencing

Structural variants play an important role in evolutionary processes. Besides, they constitute a large source of inter individual genetic variation that might represent a major factor in the aetiology of complex, multifactorial traits. Their importance in adaptation is becoming increasingly evident in literature. Yet, the characterization of the genomic landscape of structural variants in local breeds remains scarce to date. Herein, we investigate patterns and gene annotation of structural variants in the Creole cattle from Guadeloupe breed using whole genome sequences from 23 bulls representative of the population. In total, we detected 32821 ascertained SV defining 15258 regions, representing ~ 17% of the Creole cattle genome. Among these, 6639 regions have not been previously reported in the Database of Genomic Variants archive. Average number of structural variants detected per individual in the studied population is in the same order of magnitude of that observed in indicine populations and higher than that reported in taurine breeds. We observe an important within-individual variability where approximately half of the detected structural variants have low frequency (MAF < 0.25). Most of the detected structural variants (55%) occurred in intergenic regions. Genic structural variants overlapped with 7793 genes and the predicted effect of most of them is ranked as "modifier". Among the structural variants that were predicted to have a high functional impact on the protein, a 5.5 Kb in length, highly frequent deletion on chromosome 2, affects ALPI, a gene associated with the interaction between gut microbiota and host immune system. The 6639 newly identified structural variants regions include three deletions and three duplications shared by more than 80% of individuals that are significantly enriched for genes related to tRNA threonylcarbamoyladenosine metabolic process, important for temperature adaptation in thermophilic organisms, therefore suggesting a potential role in the thermotolerance of Creole cattle from Guadeloupe cattle to tropical climate. Overall, highly frequent structural variants that are specific to the Creole cattle population encompass olfactory receptor and immunity genes as well as genes involved in muscle tone, muscle development and contraction. Beyond mapping and characterizing structural variants in the Creole cattle from Guadeloupe breed, this study provides valuable information for a better understanding of the potential role of chromosomal rearrangements in adaptive traits in cattle.

A Pilot Detection and Associate Study of Gene Presence-Absence Variation in Holstein Cattle

Presence-absence variations (PAVs) are important structural variations, wherein a genomic segment containing one or more genes is present in some individuals but absent in others. While PAVs have been extensively studied in plants, research in cattle remains limited. This study identified PAVs in 173 Holstein bulls using whole-genome sequencing data and assessed their associations with 46 economically important traits. Out of 28,772 cattle genes (from the longest transcripts), a total of 26,979 (93.77%) core genes were identified (present in all individuals), while variable genes included 928 softcore (present in 95-99% of individuals), 494 shell (present in 5-94%), and 371 cloud genes (present in <5%). Cloud genes were enriched in functions associated with hormonal and antimicrobial activities, while shell genes were enriched in immune functions. PAV-based genome-wide association studies identified associations between gene PAVs and 16 traits including milk, fat, and protein yields, as well as traits related to health and reproduction. Associations were found on multiple chromosomes, illustrating important associations on cattle chromosomes 7 and 15, involving olfactory receptor and immune-related genes, respectively. By examining the PAVs at the population level, the results of this research provided crucial insights into the genetic structures underlying the complex traits of Holstein cattle.

Identify Candidate Genes Associated with the Weight and Egg Quality Traits in Wenshui Green Shell-Laying Chickens by the Copy Number Variation-Based Genome-Wide Association Study

Copy number variation (CNV), as an essential source of genetic variation, can have an impact on gene expression, genetic diversity, disease susceptibility, and species evolution in animals. To better understand the weight and egg quality traits of chickens, this paper aimed to detect CNVs in Wenshui green shell-laying chickens and conduct a copy number variation regions (CNVRs)-based genome-wide association study (GWAS) to identify variants and candidate genes associated with their weight and egg quality traits to support related breeding efforts. In our paper, we identified 11,035 CNVRs in Wenshui green shell-laying chickens, which collectively spanned a length of 13.1 Mb, representing approximately 1.4% of its autosomal genome. Out of these CNVRs, there were 10,446 loss types, 491 gain types, and 98 mixed types. Notably, two CNVRs showed significant correlations with egg quality, while four CNVRs exhibited significant associations with body weight. These significant CNVRs are located on chromosome 4. Further analysis identified potential candidate genes that influence weight and egg quality traits, including FAM184B, MED28, LAP3, ATOH8, ST3GAL5, LDB2, and SORCS2. In this paper, the CNV map of the Wenshui green shell-laying chicken genome was constructed for the first time through population genotyping. Additionally, CNVRs can be employed as molecular markers to genetically improve chickens' weight and egg quality traits.

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

Unveiling the Genetic Landscape of Feed Efficiency in Holstein Dairy Cows: Insights into Heritability, Genetic Markers, and Pathways via Meta-Analysis

Improving the feeding efficiency of dairy cows is a key component to improve the utilization of land resources and meet the demand for high-quality protein. Advances in genomic methods and omics techniques have made it possible to breed more efficient dairy cows through genomic selection. The aim of this review is to obtain a comprehensive understanding of the biological background of feed efficiency (FE) complex traits in purebred Holstein dairy cows including heritability estimate, and genetic markers, genes, and pathways participating in FE regulation mechanism. Through a literature search, we systematically reviewed the heritability estimation, molecular genetic markers, genes, biomarkers, and pathways of traits related to feeding efficiency in Holstein dairy cows. A meta-analysis based on a random-effects model was performed to combine reported heritability estimates of FE complex. The heritability of residual feed intake, dry matter intake, and energy balance was 0.20, 0.34, and 0.22, respectively, which proved that it was reasonable to include the related traits in the selection breeding program. For molecular genetic markers, a total of 13 single-nucleotide polymorphisms and copy number variance loci, associated genes, and functions were reported to be significant across populations. A total of 169 reported candidate genes were summarized on a large scale, using a higher threshold (adjusted P value < 0.05). Then, the subsequent pathway enrichment of these genes was performed. The important genes reported in the articles were included in a gene list and the gene list was enriched by gene ontology (GO):biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis. Three GO:BP terms and four KEGG terms were statistically significant, which mainly focused on adenosine triphosphate (ATP) synthesis, electron transport chain, and OXPHOS pathway. Among these pathways, involved genes such as ATP5MC2, NDUFA, COX7A2, UQCR, and MMP are particularly important as they were previously reported. Twenty-nine reported biological mechanisms along with involved genes were explained mainly by four biological pathways (insulin-like growth factor axis, lipid metabolism, oxidative phosphorylation pathways, tryptophan metabolism). The information from this study will be useful for future studies of genomic selection breeding and genetic structures influencing animal FE. A better understanding of the underlying biological mechanisms would be beneficial, particularly as it might address genetic antagonism.

Genome-wide association studies for economically important traits in mink using copy number variation

Copy number variations (CNVs) are structural variants consisting of duplications and deletions of DNA segments, which are known to play important roles in the genetics of complex traits in livestock species. However, CNV-based genome-wide association studies (GWAS) have remained unexplored in American mink. Therefore, the purpose of the current study was to investigate the association between CNVs and complex traits in American mink. A CNV-based GWAS was performed with the ParseCNV2 software program using deregressed estimated breeding values of 27 traits as pseudophenotypes, categorized into traits of growth and feed efficiency, reproduction, pelt quality, and Aleutian disease tests. The study identified a total of 10,137 CNVs (6968 duplications and 3169 deletions) using the Affymetrix Mink 70K single nucleotide polymorphism (SNP) array in 2986 American mink. The association analyses identified 250 CNV regions (CNVRs) associated with at least one of the studied traits. These CNVRs overlapped with a total of 320 potential candidate genes, and among them, several genes have been known to be related to the traits such as ARID1B, APPL1, TOX, and GPC5 (growth and feed efficiency traits); GRM1, RNASE10, WNT3, WNT3A, and WNT9B (reproduction traits); MYO10, and LIMS1 (pelt quality traits); and IFNGR2, APEX1, UBE3A, and STX11 (Aleutian disease tests). Overall, the results of the study provide potential candidate genes that may regulate economically important traits and therefore may be used as genetic markers in mink genomic breeding programs.

kGWASflow: a modular, flexible, and reproducible Snakemake workflow for k-mers-based GWAS

Genome-wide association studies (GWAS) have been widely used to identify genetic variation associated with complex traits. Despite its success and popularity, the traditional GWAS approach comes with a variety of limitations. For this reason, newer methods for GWAS have been developed, including the use of pan-genomes instead of a reference genome and the utilization of markers beyond single-nucleotide polymorphisms, such as structural variations and k-mers. The k-mers-based GWAS approach has especially gained attention from researchers in recent years. However, these new methodologies can be complicated and challenging to implement. Here, we present kGWASflow, a modular, user-friendly, and scalable workflow to perform GWAS using k-mers. We adopted an existing kmersGWAS method into an easier and more accessible workflow using management tools like Snakemake and Conda and eliminated the challenges caused by missing dependencies and version conflicts. kGWASflow increases the reproducibility of the kmersGWAS method by automating each step with Snakemake and using containerization tools like Docker. The workflow encompasses supplemental components such as quality control, read-trimming procedures, and generating summary statistics. kGWASflow also offers post-GWAS analysis options to identify the genomic location and context of trait-associated k-mers. kGWASflow can be applied to any organism and requires minimal programming skills. kGWASflow is freely available on GitHub (https://github.com/akcorut/kGWASflow) and Bioconda (https://anaconda.org/bioconda/kgwasflow).

Finding Predictors of Leg Defects in Pigs Using CNV-GWAS

One of the most important areas of modern genome research is the search for meaningful relationships between genetic variants and phenotypes. In the livestock field, there has been research demonstrating the influence of copy number variants (CNVs) on phenotypic variation. Despite the wide range in the number and size of detected CNVs, a significant proportion differ between breeds and their functional effects are underestimated in the pig industry. In this work, we focused on the problem of leg defects in pigs (lumps/growths in the area of the hock joint on the hind legs) and focused on searching for molecular genetic predictors associated with this trait for the selection of breeding stock. The study was conducted on Large White pigs using three CNV calling tools (PennCNV, QuantiSNP and R-GADA) and the CNVRanger association analysis tool (CNV-GWAS). As a result, the analysis identified three candidate CNVRs associated with the formation of limb defects. Subsequent functional analysis suggested that all identified CNVs may act as potential predictors of the hock joint phenotype of pigs. It should be noted that the results obtained indicate that all significant regions are localized in genes (CTH, SRSF11, MAN1A1 and LPIN1) responsible for the metabolism of amino acids, fatty acids, glycerolipids and glycerophospholipids, thereby related to the immune response, liver functions, content intramuscular fat and animal fatness. These results are consistent with previously published studies, according to which a predisposition to the formation of leg defects can be realized through genetic variants associated with the functions of the liver, kidneys and hematological characteristics.

Association between copy number variation of SERPINA3-1 gene and growth traits in Chinese cattle

Serpin peptidase inhibitor, clade A (alpha-1 antiproteinase, antitrypsin), member 3 (SERPINA3) belongs to the serine protease inhibitor family A subtype, and contains 8 genes from SERPINA3-1 to SERPINA3-8. Although the regulatory effects of these 8 genes have been revealed one by one in recent years, the related effects of SERPINA3-1 gene on cattle growth is still unclear. This study used quantitative Real time PCR (qPCR) to detect the type of copy number variation (CNV) of SERPINA3-1 gene in a total of 542 Chinese cattle, and expression of SERPINA3-1 gene in different tissues of Qinchuan cattles (adult) on mRNA level. Then association analysis was conducted between the detection results and cattle growth traits. The results showed that the Duplication type in SERPINA3-1 gene performed better on the growth traits and the CNV was significantly correlated with multiple growth traits (p < 0.05). In addition, SERPINA3-1 gene has different expression conditions in different tissues, results showed that SERPINA3-1 gene has a low expression in muscle. In conclusion, we speculate that the SERPINA3-1 gene can be used as a molecular marker and the result of this study could be a basic material for candidate functional genes for beef cattle growth and development.

Genome wide copy number variations using Porcine 60K SNP Beadchip in Landlly pigs

In the present study, Porcine 60K SNP genotype data from 69 Landlly pigs were used to explore Copy Number Variations (CNVs) across the autosomes. A total of 386 CNVs were identified using Hidden Markov Model (HMM) in PennCNV software, which were subsequently aggregated to 115 CNV regions (CNVRs). Among the total detected CNVRs, 58 gain, 49 were loss type while remaining 8 events were both gain and loss types. Identified CNVRs covered 12.5 Mb (0.55%) of Sus scrofa reference 11.1 genome. Comparison of our results with previous investigations on pigs revealed that approximately 75% CNVRs were novel, which may be due to differences in genetic background, environment and implementation of artificial selection in Landlly pigs. Functional annotation and pathway analysis showed the significant enrichment of 267 well-annotated Sus scrofa genes in CNVRs. These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity. Additionally, KIT and FUT1 were two major genes detected on CNVR in our population. This investigation provided a comprehensive overview of CNV distribution in the Indian porcine genome for the first time, which may be useful for further investigating the association of important quantitative traits in Landlly pigs.Highlights115 CNVRs were identified in 69 Landlly pig population.Approximately 75% detected CNVRs were novel for Landlly population.Significant enrichment of 267 well-annotated Sus scrofa genes observed in these CNVRs.These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity.Comprehensive CNV map in the Indian porcine genome developed for the first time.

Heritability and variance component estimation for feed and water intake behaviors of feedlot cattle

Feed and water intake are two important aspects of cattle production that greatly impact the profitability, efficiency, and sustainability of producers. Feed and, to a lesser degree, water intake have been studied previously; however, there is little research on their associated animal behaviors and there is a lack of standardized phenotypes for these behaviors. Feed and water intakes obtained with an Insentec system (Hokofarm Group, The Netherlands) from 830 crossbred steers were used to compute five intake behaviors for both feed and water: daily sessions (DS), intake rate (IR), session size (SS), time per session (TS), and session interval (SI). Variance components and heritabilities were estimated for each trait. Heritabilities for feed intake behaviors were 0.50 ± 0.12, 0.63 ± 0.12, 0.40 ± 0.13, 0.35 ± 0.12, and 0.60 ± 0.12 for DS, IR, SS, TS, and SI, respectively. Heritabilities for water intake behaviors were 0.56 ± 0.11, 0.88 ± 0.07, 0.70 ± 0.11, 0.54 ± 0.12, and 0.80 ± 0.10 for NS, IR, SS, TS, and SI, respectively. Daily dry matter intake (DDMI) and daily water intake (DWI) had heritabilities of 0.57 ± 0.11 and 0.44 ± 0.11. Phenotypic correlations varied between pairs of traits (-0.83 to 0.82). Genetic correlations between DDMI and feed intake behaviors were moderate to high, while genetic correlations between DWI and water intake behaviors were low to moderate. Several significant single nucleotide polymorphisms (SNP) were identified for the feed and water intake behaviors. Genes and previously reported quantitative trait loci near significant SNPs were evaluated. The results indicated that feed and water intake behaviors are influenced by genetic factors and are heritable, providing one additional route to evaluate or manipulate feed and water intake.

Genome-wide elucidation of CNV regions and their association with production and reproduction traits in composite Vrindavani cattle

The present study was aimed to analyze the genome-wide copy number variations (CNVs) in Vrindavani composite cattle and concatenate them into CNV regions (CNVRs), and finally test the association of CNVRs with different production and reproduction traits. Genotypic data, generated on BovineSNP50 Beadchip (v3) array for 96 Vrindavani animals, was used to elucidate the CNVs at the genome level. Intensity data covering over 53,218 SNP genotypes on bovine genome was used. Algorithm based on Hidden Markov Model was employed in PennCNV program to detect, normalize and filter CNVs across the genome. 252 putative CNVs, detected via PennCNV program, in different individuals were concatenated into 71 CNV regions (CNVRs) using CNVRuler program. Association of CNVRs with important (re)production traits in Vrindavani animals was assessed using linear regression. Five CNVRs were found to be significantly associated with ten important (re)production traits. The genes harbored in these regions provided useful insights into the association of CNVRs with genes and ultimately the variation at phenotype level. Important genes that overlapped with CNVRs included WASHC4, HS6ST3, MBNL2, TOLLIP, PIDD1 and TSPAN4. Furthermore, the CNVRs were found to overlap with important QTLs available in AnimalQTL database which affect milk yield and composition along with reproduction and immune function traits. The copy number states of three enes were validated using digital droplet PCR technique. The results from the present study significantly enhance the understanding about CNVs in Vrindavani cattle and should help establish its CNV map. The study will also enable further investigation on association of these variants with important traits of economic interest including disease incidence.

Copy Number Variation (CNV): A New Genomic Insight in Horses

Simple Summary

This study aimed to contribute to our knowledge of CNVs, a type of genomic marker in equines, by producing, for the first time, a fine-scale characterization of the CNV regions (CNVRs) in the Pura Raza Española horse breed. We found not only the existence of a unique pattern of genomic regions enriched in CNVs in the PRE in comparison with the data available from other breeds but also the incidence of CNVs across the entire genome. Since these regions could affect the structure and dose of the genes involved, we also performed a gene ontology analysis which revealed that most of the genes overlapping in CNVRs were related to the olfactory pathways and immune response.

Abstract

Copy number variations (CNVs) are a new-fangled source of genetic variation that can explain changes in the phenotypes in complex traits and diseases. In recent years, their study has increased in many livestock populations. However, the study and characterization of CNVs in equines is still very limited. Our study aimed to investigate the distribution pattern of CNVs, characterize CNV regions (CNVRs), and identify the biological pathways affected by CNVRs in the Pura Raza Española (PRE) breed. To achieve this, we analyzed high-density SNP genotyping data (670,804 markers) from a large cohort of 654 PRE horses. In total, we identified 19,902 CNV segments and 1007 CNV regions in the whole population. The length of the CNVs ranged from 1.024 kb to 4.55 Mb, while the percentage of the genome covered by CNVs was 4.4%. Interestingly, duplications were more abundant than deletions and mixed CNVRs. In addition, the distribution of CNVs across the chromosomes was not uniform, with ECA12 being the chromosome with the largest percentage of its genome covered (19.2%), while the highest numbers of CNVs were found in ECA20, ECA12, and ECA1. Our results showed that 71.4% of CNVRs contained genes involved in olfactory transduction, olfactory receptor activity, and immune response. Finally, 39.1% of the CNVs detected in our study were unique when compared with CNVRs identified in previous studies. To the best of our knowledge, this is the first attempt to reveal and characterize the CNV landscape in PRE horses, and it contributes to our knowledge of CNVs in equines, thus facilitating the understanding of genetic and phenotypic variations in the species. However, further research is still needed to confirm if the CNVs observed in the PRE are also linked to variations in the specific phenotypical differences in the breed.

Genome-wide association study of beef bull semen attributes

Background

Cattle production is dependent upon fertility because it results in producing offspring to offset production costs. A number of semen attributes are believed to affect fertility and are frequently measured as part of routine breeding soundness exams or semen collection procedures. The objective of this study was to perform a single-step genome-wide association study (ssGWAS) for beef bull semen attributes. Beef bull fertility phenotypes including volume (VOL), concentration (CONC), number of spermatozoa (NSP), initial motility (IMot), post-thaw motility (PTMot), three-hour post-thaw motility (3HRPTMot), percentage of normal spermatozoa (%NORM), primary abnormalities (PRIM), and secondary abnormalities (SEC) were obtained from two artificial insemination (AI) centers. A total of 1819 Angus bulls with 50,624 collection records were used for ssGWAS. A five-generation pedigree was obtained from the American Angus Association and consisted of 6521 sires and 17,136 dams. Genotypes on 1163 bulls were also obtained from the American Angus Association and utilized in ssGWAS.

Results

A multi-trait animal model was used for the estimation of single nucleotide polymorphism (SNP) effects. Significant SNP were those with a -log10P-value threshold greater than 4.0. Volume, CONC, NSP, IMot, PTMot, 3HRPTMot, %NORM, PRIM, and SEC have five, three, six, seven, two, six, six, and two genome-wide significant SNP, respectively.

Conclusions

Several significant SNP were determined to be near or within quantitative trait loci (QTL) associated with beef bull semen attributes. In addition, genes associated with fertility were found to contain or be near the significant SNP found in the study. The results indicate there are regions of the genome that impact fertility, proving inclusion of genomic information into genetic evaluation should be advantageous for genetic improvement of male fertility traits.

Copy number variations of the KAT6A gene are associated with body measurements of Chinese sheep breeds

Copy number variation (CNV) is one kind of genomic structure variations and presents as gains and losses of genomic fragments. More recently, we have made an atlas of CNV maps for livestock. In the future, it is a primary focus to determine the phenotypic effects of candidate CNVs. Lysine Acetyltransferase 6 A (KAT6A) is a protein coding gene and plays a critical role in many cellular processes. However, the effects of KAT6A CNVs on sheep body measurements remains unknown. In this study, we performed quantitative polymerase chain reaction (qPCR) to detect the presences and distributions of three CNV regions within KAT6A gene in 672 sheep from four Chinese breeds. Association analysis indicated that the three CNVs of KAT6A gene were significantly associated with body measurement(s) in Small-tailed Han sheep (STH) and Hu sheep (HU) (p < 0.05), while no effects on Large-tailed Han sheep (LTH) were observed (p > 0.05) were observed. Additionally, only one CNV was significantly associated with body measurement (body length) in Chaka sheep (CK) (p < 0.05). Our study provided evidence that the CNV(s) of KAT6A gene could be used as candidate marker(s) for molecular breedings of STH, HU, and CK breeds.

Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep

The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.

Abnormal cleavage is involved in the self-correction of bovine preimplantation embryos

Chromosome instability leading to aneuploidy during early cleavage is well known in humans and cattle. Partial compaction (PC), which occurs only in some blastomeres, is suggested as a self-correction mechanism through which human embryos avoid aneuploid mosaicism. Partially compacted embryos show abnormal cleavages more frequently during early development; however, the mechanism by which blastomeres are excluded has not been elucidated. Here, we confirmed PC in approximately half of the tested bovine embryos, similar to that in human embryos. DNA sequencing of single-cell and intact embryos revealed that the morulae that excluded some blastomeres had euploidy, but many of the excluded blastomeres had aneuploidy. Time-lapse imaging of zygotes without the zona pellucida revealed that the excluded blastomeres underwent reverse and direct cleavages, which are abnormal cleavages, more frequently than the blastomeres involved in compaction. These results suggest the potential role of abnormal cleavage in the self-correction mechanism during the development of mammalian preimplantation embryos.

Genome-wide association study between copy number variants and hoof health traits in Holstein dairy cattle

Genome-wide association studies based on SNP have been completed for multiple traits in dairy cattle; however, copy number variants (CNV) could add genomic information that has yet to be harnessed. The objectives of this study were to identify CNV in genotyped Holstein animals and assess their association with hoof health traits using deregressed estimated breeding values as pseudophenotypes. A total of 23,256 CNV comprising 1,645 genomic regions were identified in 5,845 animals. Fourteen genomic regions harboring structural variations, including 9 deletions and 5 duplications, were associated with at least 1 of the studied hoof health traits. This group of traits included digital dermatitis, interdigital dermatitis, heel horn erosion, sole ulcer, white line lesion, sole hemorrhage, and interdigital hyperplasia; no regions were associated with toe ulcer. Twenty candidate genes overlapped with the regions associated with these traits including SCART1, NRXN2, KIF26A, GPHN, and OR7A17. In this study, an effect on infectious hoof lesions could be attributed to the PRAME (Preferentially Expressed Antigen in Melanoma) gene. Almost all genes detected in association with noninfectious hoof lesions could be linked to known metabolic disorders. The knowledge obtained considering information of associated CNV to the traits of interest in this study could improve the accuracy of estimated breeding values. This may further increase the genetic gain for these traits in the Canadian Holstein population, thus reducing the involuntary animal losses due to lameness.

Gene Banks as Reservoirs to Detect Recent Selection: The Example of the Asturiana de los Valles Bovine Breed

Gene banks, framed within the efforts for conserving animal genetic resources to ensure the adaptability of livestock production systems to population growth, income, and climate change challenges, have emerged as invaluable resources for biodiversity and scientific research. Allele frequency trajectories over the few last generations contain rich information about the selection history of populations, which cannot be obtained from classical selection scan approaches based on present time data only. Here we apply a new statistical approach taking advantage of genomic time series and a state of the art statistic (nSL) based on present time data to disentangle both old and recent signatures of selection in the Asturiana de los Valles cattle breed. This local Spanish originally multipurpose breed native to Asturias has been selected for beef production over the last few generations. With the use of SNP chip and whole-genome sequencing (WGS) data, we detect candidate regions under selection reflecting the effort of breeders to produce economically valuable beef individuals, e.g., by improving carcass and meat traits with genes such as MSTN, FLRT2, CRABP2, ZNF215, RBPMS2, OAZ2, or ZNF609, while maintaining the ability to thrive under a semi-intensive production system, with the selection of immune (GIMAP7, GIMAP4, GIMAP8, and TICAM1) or olfactory receptor (OR2D2, OR2D3, OR10A4, and 0R6A2) genes. This kind of information will allow us to take advantage of the invaluable resources provided by gene bank collections from local less competitive breeds, enabling the livestock industry to exploit the different mechanisms fine-tuned by natural and human-driven selection on different populations to improve productivity.

Genetic mechanisms underlying feed utilization and implementation of genomic selection for improved feed efficiency in dairy cattle

The economic importance of genetically improving feed efficiency has been recognized by cattle producers worldwide. It has the potential to considerably reduce costs, minimize environmental impact, optimize land and resource use efficiency, and improve the overall cattle industry’s profitability. Feed efficiency is a genetically complex trait that can be described as units of product output (e.g., milk yield) per unit of feed input. The main objective of this review paper is to present an overview of the main genetic and physiological mechanisms underlying feed utilization in ruminants and the process towards implementation of genomic selection for feed efficiency in dairy cattle. In summary, feed efficiency can be improved via numerous metabolic pathways and biological mechanisms through genetic selection. Various studies have indicated that feed efficiency is heritable, and genomic selection can be successfully implemented in dairy cattle with a large enough training population. In this context, some organizations have worked collaboratively to do research and develop training populations for successful implementation of joint international genomic evaluations. The integration of “-omics” technologies, further investments in high-throughput phenotyping, and identification of novel indicator traits will also be paramount in maximizing the rates of genetic progress for feed efficiency in dairy cattle worldwide.

Insight into the Possible Formation Mechanism of the Intersex Phenotype of Lanzhou Fat-Tailed Sheep Using Whole-Genome Resequencing

Simple Summary

Individuals with hermaphroditism are a serious hazard to animal husbandry and production due to their abnormal fertility. The molecular mechanism of sheep intersex formation was unclear. This study was the first to locate the homologous sequence of the goat polled intersex syndrome (PIS) region in sheep and found that the intersex traits of Lanzhou fat-tailed sheep were not caused by the lack of this region. By detecting the selective sweep regions, vital genes associated with androgen biosynthesis and the follicle stimulating hormone response entry were found, including steroid 5 alpha-reductase 2 (SRD5A2), and pro-apoptotic WT1 regulator (PAWR). Additionally, the copy number variations of the four regions on chr9, chr1, chr4, and chr16 may affect the expression of the gonadal development genes, zinc finger protein, FOG family member 2 (ZFPM2), LIM homeobox 8 (LHX8), inner mitochondrial membrane peptidase subunit 2 (IMMP2L) and slit guidance ligand 3 (SLIT3), respectively, and further affect the formation of intersex traits.

Abstract

Intersex, also known as hermaphroditism, is a serious hazard to animal husbandry and production. The mechanism of ovine intersex formation is not clear. Therefore, genome-wide resequencing on the only two intersex and two normal Lanzhou fat-tailed (LFT) sheep, an excellent but endangered Chinese indigenous sheep breed, was performed. Herein, the deletion of homologous sequences of the goat polled intersex syndrome (PIS) region (8787 bp, 247747059–247755846) on chromosome 1 of the LFT sheep was not the cause of the ovine intersex trait. By detecting the selective sweep regions, we found that the genes related to androgen biosynthesis and follicle stimulating hormone response items, such as steroid 5 alpha-reductase 2 (SRD5A2), steroid 5 alpha-reductase 3 (SRD5A3), and pro-apoptotic WT1 regulator (PAWR), may be involved in the formation of intersex traits. Furthermore, the copy number variations of the four regions, chr9: 71660801–71662800, chr1: 50776001–50778000, chr4: 58119201–58121600, and chr16: 778801–780800, may affect the expression of the zinc finger protein, FOG family member 2 (ZFPM2), LIM homeobox 8 (LHX8), inner mitochondrial membrane peptidase subunit 2 (IMMP2L) and slit guidance ligand 3 (SLIT3) genes, respectively, which contribute to the appearance of intersex traits. These results may supply a theoretical basis for the timely detection and elimination of intersex individuals in sheep, which could accelerate the healthy development of animal husbandry.

High confidence copy number variants identified in Holstein dairy cattle from whole genome sequence and genotype array data

Multiple methods to detect copy number variants (CNV) relying on different types of data have been developed and CNV have been shown to have an impact on phenotypes of numerous traits of economic importance in cattle, such as reproduction and immunity. Further improvements in CNV detection are still needed in regard to the trade-off between high-true and low-false positive variant identification rates. Instead of improving single CNV detection methods, variants can be identified in silico with high confidence when multiple methods and datasets are combined. Here, CNV were identified from whole-genome sequences (WGS) and genotype array (GEN) data on 96 Holstein animals. After CNV detection, two sets of high confidence CNV regions (CNVR) were created that contained variants found in both WGS and GEN data following an animal-based (n = 52) and a population-based (n = 36) pipeline. Furthermore, the change in false positive CNV identification rates using different GEN marker densities was evaluated. The population-based approach characterized CNVR, which were more often shared among animals (average 40% more samples per CNVR) and were more often linked to putative functions (48 vs 56% of CNVR) than CNV identified with the animal-based approach. Moreover, false positive identification rates up to 22% were estimated on GEN information. Further research using larger datasets should use a population-wide approach to identify high confidence CNVR.

Copy Number Variations and Expression Levels of Guanylate-Binding Protein 6 Gene Associated with Growth Traits of Chinese Cattle

Association studies have indicated profound effects of copy number variations (CNVs) on various phenotypes in different species. In this study, we identified the CNV distributions and expression levels of guanylate-binding protein 6 (GBP6) associated with the growth traits of Chinese cattle. The results showed that the phenotypic values of body size and weight of Xianan (XN) cattle were higher than those of Nanyang (NY) cattle. The medium CNV types were mostly identified in the XN and NY breeds, but their CNV distributions were significantly different (adjusted p < 0.05). The association analysis revealed that the body weight, cannon circumference and chest circumference of XN cattle had significantly different values in different CNV types (p < 0.05), with CNV gain types (Log22-ΔΔCt > 0.5) displaying superior phenotypic values. We also found that transcription levels varied in different tissues (p < 0.001) and the CNV gain types showed the highest relative gene expression levels in the muscle tissue, consistent with the highest phenotypic values of body weight and cannon circumference among the three CNV types. Consequently, our results suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.

Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred From Genome-Wide Copy Number Variations in Chinese Indigenous Cattle

Copy number variations (CNVs) have been demonstrated as crucial substrates for evolution, adaptation and breed formation. Chinese indigenous cattle breeds exhibit a broad geographical distribution and diverse environmental adaptability. Here, we analyzed the population structure and adaptation to high altitude of Chinese indigenous cattle based on genome-wide CNVs derived from the high-density BovineHD SNP array. We successfully detected the genome-wide CNVs of 318 individuals from 24 Chinese indigenous cattle breeds and 37 yaks as outgroups. A total of 5,818 autosomal CNV regions (683 bp-4,477,860 bp in size), covering ~14.34% of the bovine genome (UMD3.1), were identified, showing abundant CNV resources. Neighbor-joining clustering, principal component analysis (PCA), and population admixture analysis based on these CNVs support that most Chinese cattle breeds are hybrids of Bos taurus taurus (hereinafter to be referred as Bos taurus) and Bos taurus indicus (Bos indicus). The distribution patterns of the CNVs could to some extent be related to the geographical backgrounds of the habitat of the breeds, and admixture among cattle breeds from different districts. We analyzed the selective signatures of CNVs positively involved in high-altitude adaptation using pairwise Fst analysis within breeds with a strong Bos taurus background (taurine-type breeds) and within Bos taurus×Bos indicus hybrids, respectively. CNV-overlapping genes with strong selection signatures (at top 0.5% of Fst value), including LETM1 (Fst = 0.490), TXNRD2 (Fst = 0.440), and STUB1 (Fst = 0.420) within taurine-type breeds, and NOXA1 (Fst = 0.233), RUVBL1 (Fst = 0.222), and SLC4A3 (Fst=0.154) within hybrids, were potentially involved in the adaptation to hypoxia. Thus, we provide a new profile of population structure from the CNV aspects of Chinese indigenous cattle and new insights into high-altitude adaptation in cattle.

Functional and population genetic features of copy number variations in two dairy cattle populations

BACKGROUND

Copy Number Variations (CNVs) are gain or loss of DNA segments that are known to play a role in shaping a wide range of phenotypes. In this study, we used two dairy cattle populations, Holstein Friesian and Jersey, to discover CNVs using the Illumina BovineHD Genotyping BeadChip aligned to the ARS-UCD1.2 assembly. The discovered CNVs were investigated for their functional impact and their population genetics features.

RESULTS

We discovered 14,272 autosomal CNVs, which were aggregated into 1755 CNV regions (CNVR) from 451 animals. These CNVRs together cover 2.8% of the bovine autosomes. The assessment of the functional impact of CNVRs showed that rare CNVRs (MAF < 0.01) are more likely to overlap with genes, than common CNVRs (MAF ≥ 0.05). The Population differentiation index (Fst) based on CNVRs revealed multiple highly diverged CNVRs between the two breeds. Some of these CNVRs overlapped with candidate genes such as MGAM and ADAMTS17 genes, which are related to starch digestion and body size, respectively. Lastly, linkage disequilibrium (LD) between CNVRs and BovineHD BeadChip SNPs was generally low, close to 0, although common deletions (MAF ≥ 0.05) showed slightly higher LD (r2 = ~ 0.1 at 10 kb distance) than the rest. Nevertheless, this LD is still lower than SNP-SNP LD (r2 = ~ 0.5 at 10 kb distance).

CONCLUSIONS

Our analyses showed that CNVRs detected using BovineHD BeadChip arrays are likely to be functional. This finding indicates that CNVs can potentially disrupt the function of genes and thus might alter phenotypes. Also, the population differentiation index revealed two candidate genes, MGAM and ADAMTS17, which hint at adaptive evolution between the two populations. Lastly, low CNVR-SNP LD implies that genetic variation from CNVs might not be fully captured in routine animal genetic evaluation, which relies solely on SNP markers.

Copy Number Variations of Four Y-Linked Genes in Swamp Buffaloes

Simple Summary

The amplification of the male-specific region of the Y chromosome was a unique phenomenon during mammalian sex chromosome evolution. The Y-linked copy number variations of many species have been confirmed. However, the Y-linked copy number variations (CNVs) in water buffalo are still unknown. In this study, we investigated the copy number variations of four Y-linked genes (SRY, UTY, DBY, and OFD1Y) in buffalo. Our results showed that UTY was a single-copy gene in buffalo, while DBY, OFD1Y, and SRY exhibited copy number variations in buffalo.

Abstract

Copy number variation (CNV), a significant source of genetic diversity in the mammalian Y chromosome, is associated with the development of many complex phenotypes, such as spermatogenesis and male fertility. The contribution of Y-linked CNVs has been studied in various species, however, water buffalo has not been explored in this area and the genetic information still remains unknown. The aim of the current study was to investigate the CNVs of four Y-linked genes, including, sex determining Region of Y-Chromosome (SRY), ubiquitously transcribed tetratricopeptide repeat gene protein on the chromosome Y (UTY), DEAD-box helicase 3 Y-linked (DDX3Y, also known as DBY), and oral-facial-digital syndrome 1 Y-linked (OFD1Y) in 254 swamp buffaloes from 15 populations distributed across China, Vietnam, and Laos using quantitative real-time PCR (qPCR). Our results revealed the prevalence of a single-copy UTY gene in buffaloes. The DBY and OFD1Y represented CNVs among and within different buffalo breeds. The SRY showed CNVs only in Vietnamese and Laotian buffaloes. In conclusion, this study indicated that DBY, OFD1Y, and SRY showed CNVs, while the UTY was a single-copy gene in swamp buffaloes.

Genome-wide association study and genomic evaluation of feed efficiency traits in Japanese Black cattle using single-step genomic best linear unbiased prediction method

The objectives of this study were to better understand the genetic architecture and the possibility of genomic evaluation for feed efficiency traits by (i) performing genome-wide association studies (GWAS), and (ii) assessing the accuracy of genomic evaluation for feed efficiency traits, using single-step genomic best linear unbiased prediction (ssGBLUP)-based methods. The analyses were performed in residual feed intake (RFI), residual body weight gain (RG), and residual intake and body weight gain (RIG) during three different fattening periods. The phenotypes from 4,578 Japanese Black steers, which were progenies of 362 progeny-tested bulls and the genotypes from the bulls were used in this study. The results of GWAS showed that a total of 16, 8, and 12 gene ontology terms were related to RFI, RG, and RIG, respectively, and the candidate genes identified in RFI and RG were involved in olfactory transduction and the phosphatidylinositol signaling system, respectively. The realized reliabilities of genomic estimated breeding values were low to moderate in the feed efficiency traits. In conclusion, ssGBLUP-based method can lead to understand some biological functions related to feed efficiency traits, even with small population with genotypes, however, an alternative strategy will be needed to enhance the reliability of genomic evaluation.

FABP4 gene has a very large effect on feed efficiency in lactating Israeli Holstein cows

Improving feed efficiency (FE) is a major goal for the livestock industry. Previously, we have identified 48 SNP markers distributed over 32 genes significantly associated with residual feed intake (RFI) in Israeli Holstein male calves, the most significant of which are located in the bovine FABP4 gene. In the present study, we tested associations of eight of the FABP4 markers with RFI and feed conversion ratio (FCR), along with milk composition and feeding behavioral traits, in 114 lactating Israeli Holstein cows. Large allele effects were found, along with large contributions of FABP4 markers to the phenotypic variation [mean contribution of all significant markers (P < 0.05), 15.4 and 12.0% for RFI and FCR, respectively] and genotypic variation [means of all significant markers (P < 0.05), 75.7 and 32.4% in RFI and FCR, respectively]. However, the association of all significant FABP4 markers with FE and milk content traits was found in opposite directions, such that improved FE was accompanied by decreased milk content. Hence, before inclusion in breeding programs, the gain in FE must be economically balanced with the loss in milk contents. On the other hand, these findings imply that in any current improvement program concentrated on milk traits alone, without taking into account the effect on FE, the progress in milk composition is probably accompanied by deterioration of FE. These results, if confirmed in other populations and breeds, set FABP4 as a prime candidate in any marker-assisted selection program targeting FE as a whole and RFI in particular.

DNA Copy Number Variations as Markers of Mutagenic Impact

DNA copy number variation (CNV) occurs due to deletion or duplication of DNA segments resulting in a different number of copies of a specific DNA-stretch on homologous chromosomes. Implications of CNVs in evolution and development of different diseases have been demonstrated although contribution of environmental factors, such as mutagens, in the origin of CNVs, is poorly understood. In this review, we summarize current knowledge about mutagen-induced CNVs in human, animal and plant cells. Differences in CNV frequencies induced by radiation and chemical mutagens, distribution of CNVs in the genome, as well as adaptive effects in plants, are discussed. Currently available information concerning impact of mutagens in induction of CNVs in germ cells is presented. Moreover, the potential of CNVs as a new endpoint in mutagenicity test-systems is discussed.

SNV discovery and functional candidate gene identification for milk composition based on whole genome resequencing of Holstein bulls with extremely high and low breeding values

We have sequenced the whole genomes of eight proven Holstein bulls from the four half-sib or full-sib families with extremely high and low estimated breeding values (EBV) for milk protein percentage (PP) and fat percentage (FP) using Illumina re-sequencing technology. Consequently, 2.3 billion raw reads were obtained with an average effective depth of 8.1×. After single nucleotide variant (SNV) calling, total 10,961,243 SNVs were identified, and 57,451 of them showed opposite fixed sites between the bulls with high and low EBVs within each family (called as common differential SNVs). Next, we annotated the common differential SNVs based on the bovine reference genome, and observed that 45,188 SNVs (78.70%) were located in the intergenic region of genes and merely 11,871 SNVs (20.67%) located within the protein-coding genes. Of them, 13,099 common differential SNVs that were within or close to protein-coding genes with less than 5 kb were chosen for identification of candidate genes for milk compositions in dairy cattle. By integrated analysis of the 2,657 genes with the GO terms and pathways related to protein and fat metabolism, and the known quantitative trait loci (QTLs) for milk protein and fat traits, we identified 17 promising candidate genes: ALG14, ATP2C1, PLD1, C3H1orf85, SNX7, MTHFD2L, CDKN2D, COL5A3, FDX1L, PIN1, FIG4, EXOC7, LASP1, PGS1, SAO, GPLD1 and MGEA5. Our findings provided an important foundation for further study and a prompt for molecular breeding of dairy cattle.

Array CGH-based detection of CNV regions and their potential association with reproduction and other economic traits in Holsteins

BACKGROUND

Copy number variations (CNVs) are structural variants consisting of large-scale insertions and deletions of genomic fragments. Exploring CNVs and estimating their effects on phenotypes are useful for genome selection but remain challenging in the livestock.

RESULTS

We identified 1043 CNV regions (CNVRs) from array comparative genomic hybridization (CGH) data of 47 Holstein bulls. Using a probe-based CNV association approach, we detected 87 CNVRs significantly (Bonferroni-corrected P value < 0.05) associated with at least one out of 41 complex traits. Within them, 39 CNVRs were simultaneously associated with at least 2 complex traits. Notably, 24 CNVRs were markedly related to daughter pregnancy rate (DPR). For example, CNVR661 containing CYP4A11 and CNVR213 containing CTR9, respectively, were associated with DPR and other traits related to reproduction, production, and body conformation. CNVR758 was also significantly related to DPR, with a nearby gene CAPZA3, encoding one of F-actin-capping proteins which play a role in determining sperm architecture and male fertility. We corroborated these CNVRs by examining their overlapped quantitative trait loci and comparing with previously published CNV results.

CONCLUSION

To our knowledge, this is one of the first genome-wide association studies based on CNVs called by array CGH in Holstein cattle. Our results contribute substantial information about the potential CNV impacts on reproduction, health, production, and body conformation traits, which lay the foundation for incorporating CNV into the future dairy cattle breeding program.

A global analysis of CNVs in diverse yak populations using whole-genome resequencing

Background

Genomic structural variation represents a source for genetic and phenotypic variation, which may be subject to selection during the environmental adaptation and population differentiation. Here, we described a genome-wide analysis of copy number variations (CNVs) in 16 populations of yak based on genome resequencing data and CNV-based cluster analyses of these populations.

Results

In total, we identified 51,461 CNV events and defined 3174 copy number variation regions (CNVRs) that covered 163.8 Mb (6.2%) of yak genome with more “loss” events than both “gain” and “both” events, and we confirmed 31 CNVRs in 36 selected yaks using quantitative PCR. Of the total 163.8 Mb CNVR coverage, a 10.8 Mb region of high-confidence CNVRs directly overlapped with the 52.9 Mb of segmental duplications, and we confirmed their uneven distributions across chromosomes. Furthermore, functional annotation indicated that the CNVR-harbored genes have a considerable variety of molecular functions, including immune response, glucose metabolism, and sensory perception. Notably, some of the identified CNVR-harbored genes associated with adaptation to hypoxia (e.g., DCC, MRPS28, GSTCD, MOGAT2, DEXI, CIITA, and SMYD1). Additionally, cluster analysis, based on either individuals or populations, showed that the CNV clustering was divided into two origins, indicating that some yak CNVs are likely to arisen independently in different populations and contribute to population difference.

Conclusions

Collectively, the results of the present study advanced our understanding of CNV as an important type of genomic structural variation in yak, and provide a useful genomic resource to facilitate further research on yak evolution and breeding.

Electronic supplementary material

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