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

Copy number variant scan in more than four thousand Holstein cows bred in Lombardy, Italy

Copy Number Variants (CNV) are modifications affecting the genome sequence of DNA, for instance, they can be duplications or deletions of a considerable number of base pairs (i.e., greater than 1000 bp and up to millions of bp). Their impact on the variation of the phenotypic traits has been widely demonstrated. In addition, CNVs are a class of markers useful to identify the genetic biodiversity among populations related to adaptation to the environment. The aim of this study was to detect CNVs in more than four thousand Holstein cows, using information derived by a genotyping done with the GGP (GeneSeek Genomic Profiler) bovine 100K SNP chip. To detect CNV the SVS 8.9 software was used, then CNV regions (CNVRs) were detected. A total of 123,814 CNVs (4,150 non redundant) were called and aggregated into 1,397 CNVRs. The PCA results obtained using the CNVs information, showed that there is some variability among animals. For many genes annotated within the CNVRs, the role in immune response is well known, as well as their association with important and economic traits object of selection in Holstein, such as milk production and quality, udder conformation and body morphology. Comparison with reference revealed unique CNVRs of the Holstein breed, and others in common with Jersey and Brown. The information regarding CNVs represents a valuable resource to understand how this class of markers may improve the accuracy in prediction of genomic value, nowadays solely based on SNPs markers.

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.

Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds

BACKGROUND

Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated.

RESULTS

Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds.

CONCLUSION

The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.

Genomic Regions Associated with Milk Composition and Fertility Traits in Spring-Calved Dairy Cows in New Zealand

The objective of this study was to identify genomic regions and genes that are associated with the milk composition and fertility traits of spring-calved dairy cows in New Zealand. Phenotypic data from the 2014–2015 and 2021–2022 calving seasons in two Massey University dairy herds were used. We identified 73 SNPs that were significantly associated with 58 potential candidate genes for milk composition and fertility traits. Four SNPs on chromosome 14 were highly significant for both fat and protein percentages, and the associated genes were DGAT1, SLC52A2, CPSF1, and MROH1. For fertility traits, significant associations were detected for intervals from the start of mating to first service, the start of mating to conception, first service to conception, calving to first service, and 6-wk submission, 6-wk in-calf, conception to first service in the first 3 weeks of the breeding season, and not in calf and 6-wk calving rates. Gene Ontology revealed 10 candidate genes (KCNH5, HS6ST3, GLS, ENSBTAG00000051479, STAT1, STAT4, GPD2, SH3PXD2A, EVA1C, and ARMH3) that were significantly associated with fertility traits. The biological functions of these genes are related to reducing the metabolic stress of cows and increasing insulin secretion during the mating period, early embryonic development, foetal growth, and maternal lipid metabolism during the pregnancy period.

Haplotype-based association study of production and reproduction traits in multigenerational Vrindavani population

Haplotype-based association analysis promises to reveal important information regarding the effect of genetic variants on economic traits of interest. The present study aimed to evaluate the haplotype structure of Vrindavani cattle and explore the association of haplotypes with (re)production traits of economic interest. Genotyping array data of medium density (Bovine50KSNP BeadChip) on 96 randomly selected Vrindavani cows was used in the present study. Genotypes were called in GenomeStudio program while quality control was undertaken in PLINK using standard thresholds. The phenotypic traits used in the present study included age at first calving, dry days, lactation length, peak yield, total lactation milk yield, inter-calving period and service period. The haplotype structure of Vrindavani population was assessed, using a sliding window of 20 SNP with a shift of 5 SNPs at a time, in terms of the size of haplotype blocks regarding their length (in Kb) and frequency in chromosome-wise fashion. Haplotype blocks were assessed for possible association with important production and reproduction traits across three lactation cycles in Vrindavani cattle population. The first ten principal components were included in the model for haplotype-based association analysis to correct for stratification effects of assessed individuals. Multiple haplotypes were found to be associated with age at first calving, total lactation milk yield, peak yield, dry days, inter-calving period and service period. Various candidate genes were found to overlap haplotypes that were significantly associated with age at first calving (CDH18, MARCHF11, MYO10, FBXL7), total lactation milk yield (TGF, PDE1A, and COL8A1), peak yield (PPARGC1A, RCAN1, KCNE1, SMIM34 and MRPS6), dry days (CPNE4, ACAD11 and MRAS), inter-calving period (ABCG5, ABCG8 and COX7A2L) and service period (FOXL2 and PIK3CB). The putative candidate genes overlapping the significantly associated haplotypes revealed important pathways affecting the production and reproduction performance of animals. The identified genes and pathways may serve as good candidate markers to select animals for improved production and reproduction performance in future generations.